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Introduction 



An important test of the standard model (SM) of particle physics [T] is to measure and understand 
the properties of the highest momentum-transfer particle collisions, which correspond to measure- 
ments at the shortest distances. The predicted high energy behavior of the SM, however, becomes 
unphysical at an interaction energy on the order of several TeV. These phenomena beyond the SM 
may involve new elementary particles, new fundamental forces, and/or a modification of space- 
time geometry. These new phenomena are likely to show up as an anomalous production rate of 
a combination of the known fundamental particles. 

The unknown nature of possible new phenomena in the energy range accessible at the Tevatron 
is the motivation for a search strategy that does not focus on a single model of new physics, but 
presents a wide net for new phenomena. We compared SM predictions with the rates measured 
at the Tevatron with the CDF detector for final states with at least one high-PT lepton (e or n) 
and photon (7), plus other detected objects (leptons, photons, jets, and missing transverse energy, 
Pt). a priori definition of selection cuts for the search allows to test Run I anomalies, such as 
the observation of an event consistent with the production of two energetic photons, two energetic 
electrons, and large (the "6677 ^-p event"), in Run II data. Another intriguing Run I result 
that is important to test is a 2.7a excess above the Standard Model expectations in the ^7^x 
signature [2l [3]. 

The Fermilab Tevatron has the highest center-of-mass energy collisions (per nucleon) of any 
accelerator to date, and thus has the potential to discover new physics. The upgraded CDF II 
detector provides us better solid angle coverage and particle identification. The production of two 
vector gauge bosons, precisely predicted in the Standard Model, provides a set of signatures in 
which to search for the production of new particles which couple to the SM gauge sector (the top 
quark being the last new example). 

This analysis has been done with 305 pb~^ of pp collisions at ^/s= 1.96 TeV, collected with 
CDF detector at the Tevatron, Fermilab between March 23, 2002 and August 22, 2004. The main 
results of this thesis have been published in [H El [6]. Standard Model W'j and Z7 production 
CDF Run II results are published in [7]. The status of the Lepton+Photon+X search has been 
presented at the APS Conference (Philadelphia, 2003). The results have been presented at the 
SUSY 2005 conference (Durham, 2005) [8], the International School of Subnuclear Physics (Erice, 
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2005) [5], Lake Louise Winter Institute (2006) [9], and also at the CDF Collaboration Meeting 
(Sitges, 2005) and at the Exotics, Photon and Very Exotic Phenomena working group meetings. 

At the International School of Subnuclear Physics (Erice, 2005) I havereceived the "New Tal- 
ents" Award for an Original Work in Experimental Physics for the talk "Search for New Physics 
in Photon Final states" . The work has been reviewed and approved for publication in 'S] by G. 
t'Hooft, 1999 Nobel Laureate in Physics. 

One of the most important tools for a better understanding of the events that could possibly 
be New Physics candidates is a CDF Run II Event Display visualization package [IDl [H], which 
is widely used for offline analysis as well as to monitor online data taking [12]. Development and 
Support of the EventDisplay package is a responsibility of ITEP (Moscow) group at CDF. I am 
the project leader fl3\ and responsible for this task. 

The thesis consists of an introduction, 13 chapters, and conclusions. Chapter [1] presents the 
motivation for the analysis, and gives an introduction to the Run I results and Signature-Based 
searches. Chapter [2] gives a description of the CDF experiment at the Tevatron Collider. We 
describe the CDF coordinate system, and give information about the tracking, calorimetry, muon 
and luminosity systems. We introduce the trigger and data acquisition systems. 

Chapter [3] presents a detailed description of the CDF Run II Event Display (EVD) package and 
related projects. The EVD is used for online monitoring, offline analysis and for public relation 
(PR) purposes. 

Chapter H] presents the inclusive high-px electron, muon, and photon datasets from which we 
select i-y + X candidates, as well as the time intervals of data-taking, used to test the stability of 
the event yields (Section 14. ip . It also presents an overview of the kinematic selection criteria for 
the i-f + X events (Section H^D- 

The identification criteria for objects and control samples for muon candidates are described in 
detail in Chapter [5l for electron candidates in Chapter [6], and for photon candidates in Chapter [71 
Chapter [HI describes how missing transverse energy Ot) is calculated, and gives the definition 
and describes calculation of the total transverse energy (Ht). 

Chapter [9] presents the Standard Model expectations from SM physics processes that give the 
lepton-photon signature. The primary ones are production of Wj, Zj; we include estimates 
from the two-photon (3-boson) processes Wj'-f and Z'-ff. For each of these predictions we have 
used at least 2 independent Monte Carlo generators. Backgrounds from SM processes with a 'fake' 
(misidentified other object, such as a jet) photon or lepton are described in Chapter [TOl Chapter [TT] 
gives an overview of the experimental, theoretical and luminosity systematic uncertainties. 

Chapter [T2] presents the topologies of the signatures we are looking for, and gives the number of 
events observed. Section [12.41 gives the comparison of the observed event counts with expectations 
from the sum of SM physics processes and background. In conclusion we summarize the results 
and present future prospects. 



16 



Appendix Rl presents lists of Lepton-Photon- and Multi-Lepton-Photon events (Section [A. 10 
and additional plots for £7^^ and ££7 signatures (Sections IA.2I and IA.3I) . It also presents the 
stability of the event yields for W+jets and Z+jets (Section lA.41) . distributions of the isolation 
variables for different muon types (Section IA.5I1 . Finally, it presents supplementary information 
about conversion electrons (Section IA.6I) and additional checks for non-Z backgrounds for the fifij 
signature (Section lA.71) . 
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Chapter 1 
Motivation 



The goal of elementary particle physics is to find the ultimate constituents of matter and to study 
the fundamental interactions that occur among them. To address these questions we need to per- 
form measurements at the shortest distances, and therefore to study the properties of the highest 
momentum-transfer particle collisions. Particle physics seeks a classification of the elementary 
particles and a consistent theoretical description of their interactions that leads to an accurate 
description of experimental observables. 

1.1 Standard Model, Supersymmetry, or Something Else? 

The Standard Model (SM) is an effective field theory [Ij that has so far described the fundamental 
interactions of elementary particles remarkably well. All of the data from collider experiments, 
are explained and (in principle) are calculable in the framework of the SM. However, the SM 
does not include gravity and is expected to be an effective low-energy field theory [.14j. The SM 
contains no dark matter candidate(s). The SM higgs boson mass receives quadratically divergent 
loop corrections. This results in the well-known hierarchy problem [15] of the SM. 

The different approaches to solving the hierarchy problem include eliminating the Higgs scalar 
entirely from the theory (Technicolor [16]), lowering the cutoff scale (large extra dimensions |17j). 
or embedding the Higgs field in a multiplet of a symmetry group larger than the 4D Poincare 
group (supersymmetry [18], [19] ) . 

The existence of supersymmetry (SUSY) would provide solutions to the fine tuning problem [20] . 
and possibly the hierarchy problem, which we currently encounter in the SM. The experimental 
signatures of supersymmetry are complex, as all known fermions of the SM have bosons as super- 
symmetric partners while all bosons acquire fermions as superpartners. Due to the large number 
of free parameters, it is necessary to make further assumptions in the context of specific SUSY 
models [21] for specific searches. 

Part of the SM unified "electroweak" theory of the electromagnetic and weak forces is based 
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on the exchange of four particles: the photon for electromagnetic interactions, and two charged 
W particles and a neutral Z particle for weak interactions. These particles, 7, Z^, W^, are 
fundamental in the SM. 

In searching for new particles/quantum numbers, the signature of pairs of gauge bosons {W'j, 
Zj, WW, Z^Z^, '-f'-f) is natural if pairs of particles with a conserved quantum number are pro- 
duced because of flavor conservation in QCD. W'j and SM physics processes lead to inclusive 
production of events with a high-energy lepton and a high-energy photon. 

1.2 The Lepton-Photon Events 

Besides the specific theoretical models, searching for new physics with photons has several advan- 
tages. For example, the photon is one of the three SU(2) x U(l) gauge bosons and as such is likely 
to be a good probe of new interactions since it might couple to any new gauge sector. Final-state 
photons have additional distinct detection advantages over W^ or Z° bosons since they do not 
decay. Thus they do not suffer a sensitivity loss from branching ratios and momentum sharing 
between the decay products. The photon is coupled to electric charge, and thus is radiated by 
all charged particles, including the incoming states, which is important for searching for invisible 
final states. The photon is a boson and could be produced by a fermiphobic parent. For the search 
we also require a lepton: the events containing high-i?T photon and high-P^ lepton, £7 + X, are 
rare in the SM, and therefore backgrounds are low. 

There are many models [22] of new physics that could produce i'-y + X events. Gauge-mediated 
models of supersymmetry [23j, in which the lightest super-partner (LSP) is a light gravitino, 
provide a model in which each partner of a pair of supersymmetric particles produced in a pp 
interaction decays in a chain that leads to a produced gravitino, visible as If the next-to- 
lightest neutralino (NLSP) has a photino component, each chain also can result in a photon. 
Models of supersymmetry in which the symmetry breaking is due to gravity also can produce 
decay chains with photons [2lj. For example, if the NLSP is largely photino-like, and the lightest 
is largely higgsino, decays of the former to the latter will involve the emission of a photon [25j. More 
generally, pair-production of selectrons or gauginos can result in final-states with large ^rj,, two 
photons and two leptons and lead to events like the Run I 6677 candidate event (Section ll.4p . 

For example, an initial model invoked low-energy supersymmetry with a neutralino LSP as a 
possible interpretation of the CDF Run I 6677 event [26.J via the process: 

pp e+e~(+X), e X2 + e, X?7 

where e is the selectron (the bosonic partner of the electron), and Xi and X2 the lightest 
and next-to-lightest neutralinos, the fermionic partners of the neutral bosonic states formed by 
mixing the fermionic partners W^, the B, and the neutral Higgses into mass eigenstates. 

Gauge-mediated models in which a photino decays into a gravitino are also popular choices p^. 
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and have the appeahng feature that they have a natural dark matter candidate. 

Further expanding the space of parameters, a more recent SUSY interpretation of the CDF 
/i7^'P events [I9] is resonant smuon /i production with a single dominant R-parity violating cou- 
pling (Figure ll.ip . 




Figure 1.1: Resonant smuon production and subsequent decay, producing the /i7^T events 

The current interest in models of extra dimensions [T3|, which can produce events of interest to 
the search, is a good example of an innovation that was searched for before it was conceived. 

These models predict excited states of the known standard model particles. The production of a 
pair of excited electrons [27] would provide a natural source for two photons and two electrons 
(although not unless the pair were produced with some other, undetected, particle.). As in the 
case of supersymmetry, there are many parameters in such models, with a resulting broad range 
of possible topologies with multiple gauge bosons. 

However the parameter space of SUSY models is so large, and there are so many other models 
beyond SUSY, including ones that have not yet been thought of, that we have adopted the strategy 
of testing the SM predictions in promising signatures. This strategy, the Signature- Based Search, 
is nothing more than testing the SM [1]. 

1.3 Signature-Based Searches 

While it is good to be guided by theory, one should also remain open to the unexpected. Therefore 
we use a quasi-model-independent Signature- Based Searches technique, and look for significant 
deviations from the SM [28l [29l [30] . In the Run I dataset, no significant evidence for new physics 
was found, but there were some hints that the SM may be incomplete (Section II. 4p . CDF has 
preferred to highlight some potential anomalies as worth pursuing in Run II, thus setting up 
selection criteria in a priori fashion. 

We perform the search by systematically looking at events by their final state particles. The 
strategy for the Signature- Based Search is to test the SM by looking for an excess over the SM 
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prediction. The challenge also extends to the theoretical community - to look for something new 
we will need to understand the non-new, i.e. the SM predictions, at an unprecedented level of 
precision. Some amount of this can be done with control samples - it is always best to use data 
rather than Monte Carlo, but this is not always possible. 

1.4 Run I Results and Present Analysis 
1.4.1 The 6677 Candidate Event 

In 1995 the CDF experiment, measuring pp collisions at a center-of-mass energy of 1.8 TeV at 
the Fermilab Tevatron, using 85 pb~^ of data, observed an event consistent with the production 
of two energetic photons, two energetic electrons, and large missing transverse energy [311 [32| 
EH] (Figure OD- 

eeyy^TCandidate Event 




Figure 1.2: The Run I 6677 candidate event. 



This signature is predicted to be very rare in the Standard Model of particle physics, with the 
dominant contribution being from the WW77 production: 1^1^77— > (ez/) (61^)77 —> 6677 ^t, from 
which we expect 8x10"'' events. All other sources (mostly due to detector misidentification) lead 
to 5x10^'' events. Therefore, we expect (1 ± 1) x 10^^ events, which would give us one 6677 
candidate event if we have taken one million times more data than we actually had in Run I. 
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Signature (Object) 



Obs. Expected 



> 35 GeV, |A0p^_jet| > 10° 
Njet > 4, E^' > 10 GeV, < 2.0 
6-tag, > 25 GeV 
Central 7, E^^' > 25 GeV 
Central e or jj, E^"' ^ > 25 GeV 
Central r, E^ > 25 GeV 



2 



3 







2 



1 



1 0.5 ± 0.1 



0.2 ± 0.1 



0.1 ± 0.1 



0.3 ± 0.1 



1.3 ± 0.7 



1.6 ± 0.4 



Table 1.1: Number of observed and expected 77 events with additional objects in 85 pb ^ [32] 

The event raised theoretical interest, however, as the two-lepton two-photon signature is ex- 
pected in some models of physics 'beyond the Standard Model' [1] such as gauge-mediated models 
of supersymmetry [T^ [25] . 



The detection of this single event led to the development of 'signature-based' inclusive searches 
in Run I to cast a wider net: in this case one searches for two photons + X [31l [32l [28], where X 
stands for anything, with the idea that if pairs of new particles were being created these inclusive 
signatures would be sensitive to a range of decay modes or the creation and decay of different 
particle types. 

In Run I Searches for 77-l-X, all results were consistent with the SM background expectations 
with no other exceptions other than observation of 6677 candidate event (Table 11.11) [32] . 

1.4.3 From 77 to ^7: f 7 + X Search 

Another 'signature-based' inclusive search, motivated by 6677 event was for one photon plus 
one lepton -|- X [21 [3l [29] . In general data agrees with expectations, with the exception for the 
£7^^ category. We have observed 16 ^7pT events on a background of 7.6 ± 0.7 expected. The 
16 ^7pT events consist of 11 /^7pT events and 5 67 events, versus expectations of 4.2±0.5 
and 3.4±0.3 events, respectively. The SM prediction yields the observed rate of £7^t with 0.7% 
probability (which is equivalent to 2.7 standard deviations for a Gaussian distribution). One of 
the first SUSY interpretation of the CDF /i7^T events [19] was resonant smuon /i production 
with a single dominant R-parity violating coupling(Figure [1.11) . 

The Run I search was initiated by an anomaly in the data itself, and as such the 2.7 sigma excess 
above the SM expectations must be viewed taking into account the number of such channels that 
a fluctuation could have occurred in. The Run I paper concluded: "However, an excess of events 
with 0.7% likelihood (equivalent to 2.7 standard deviations for a Gaussian distribution) in one 



1.4.2 77+X Search 
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Category 




No 


P(iV > NolfisM), % 


All i-iX 


- 


77 


- 


Z-like 67 


- 


17 


- 


Two-Body I'-yX 


24.9±2.4 


33 


9.3 


Multi-Body £7^ 


20.2±1.7 


27 


10.0 


Multi-Body U^X 


5.8 ± 0.6 


5 


68.0 


Multi-Body £77^ 


0.02±0.02 


1 


1.5 


Multi-Body l-^fj^^X 


7.6 ± 0.7 


16 


0.7 



Table 1.2: Run I Photon-Lepton Results: Number of observed and expected £7 events with 
additional objects in 85 pb~^[3] 

subsample among the five studied is an interesting result, but it is not a compelling observation of 
new physics. We look forward to more data in the upcoming run of the Fermilab Tevatron." [3]. 
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Chapter 2 

The CDF Experiment at the Tevatron 
Collider 

An important part of the study of elementary particle physics is to understand experimental 
tools - the accelerators, beams and detectors by means of which particles are accelerated, their 
trajectories controlled and their collisions studied. 

The Tevatron is currently the world's highest energy particle accelerator. Protons(p) and anti- 
protons (p) are accelerated to be brought into collision with a center of mass energy of 1.96 TeV. 
Two detectors are situated at the B0 and D0 collision points, the Collider Detector at Fermilab 
(CDF) and D0. 

2.1 The Tevatron 

Fermilab uses a series of accelerators to create the world's most energetic particle beams. The 
diagram in Figure [2TT] shows the paths taken by p and p from initial acceleration to collision in the 
Tevatron. In the first stage of acceleration H~ ions are created from the ionization of the hydrogen 
gas and accelerated to a kinetic energy of 750 KeV in the Cockcroft- Walton pre-accelerator |33j . 
The H~ ions enter a linear accelerator (Linac) [3l], where they are accelerated to 400 MeV. The 
acceleration in the Linac is done by a series of "kicks" from Radio Frequency (RF) cavities. The 
oscillating electric field of the RF cavities groups the ions into bunches. Before entering the next 
stage, a carbon foil removes the electrons from the H~ ions at injection, leaving only the protons. 
The 400 MeV protons are then injected into the circular synchrotron ("Booster"). The protons 
travel around the Booster to a final energy of 8 GeV. 

Protons are then extracted from the Booster into the Main Injector [35], where they are 
accelerated from 8 GeV to 150 GeV before the injection into the Tevatron. The Main Injector 
also produces 120 GeV protons. These protons are extracted and collide with a nickel target, 
producing a wide spectrum of secondary particles, including p. In the collisions, about 20 p are 
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produced per one million protons. The p are collected, focused, and then stored in the Accumulator 
ring. Once a sufficient number of p are collected, they are sent to the Main Injector and accelerated 
to 150 GeV. 

Finally, both the p and p are injected into the Tevatron. The Tevatron, the last stage of 
Fermilab's accelerator chain, receives 150 GeV p and p from the Main Injector and accelerates 
them to 980 GeV. The p and p travel around the Tevatron in opposite directions. The beams are 
brought to collision at the center of the two detectors, CDF II and D0 II (see Figure [2TT1) . 

2.2 The CDF Detector 

A discovery will rely heavily on a thorough understanding of the detector. Two aspects are 
critical: the identification of objects that make up each signature, and the understanding of the 
calibration and resolution of the detector. The objects for which we have a good understanding 




Figure 2.1: Layout of the Fermilab accelerator complex. Protons (solid arrow) are accelerated at 
the Cockcroft-Walton, Linac, Booster, Main Injector and finally at the Tevatron. The anti-protons 
(dashed arrow) from the anti-proton source are ffist accelerated at the Main Injector and then at 
the Tevatron. 
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of the efficiencies and fake-rates are those for which tracking is essentiah electrons, muons, and 
photons (i.e. a high confidence of the absence of a track), all in the central region. Similarly, the 
energy scale and resolutions of the calorimeters are well understood in the central region, where 
the magnetic spectrometer is used to calibrate the calorimeters. 

The CDF II detector is a cylindrically-symmetric spectrometer designed to study pp collisions 
at the Fermilab Tevatron based on the same solenoidal magnet and central calorimeters as the 
CDF I detector [36] . A cross-section of one half of the detector is shown in Figure 12. 2[ 

Because the analysis described here is intended to repeat the Run I search as closely as possible, 
we note especially the differences from the CDF I detector relevant to the detection of leptons, 
photons, and prp. The tracking systems (Section 12.2.21) used to measure the momenta of charged 
particles have been replaced with a central outer tracker (COT) that has smaller drift cells |37j . 
and an enhanced system of silicon microstrip detectors [38]. The calorimeters in the regions (Sec- 
tion [2211]) with pseudorapidity jr^l > 1 have been replaced with a more compact scintillator-based 
design, retaining the projective geometry (Section [2.2.31) . The coverage in of the central upgrade 
muon detector (CMP) and central extension muon detector (CMX) systems (Section I2.2.4p has 
been extended; the central muon detector (CMU) system is unchanged. 

The main upgrades to the CDF detector from Run I to Run II, relevant to the analysis, can be 
summarized as follows: 

• Fully digital DAQ system designed for 132 ns bunch crossing times 
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• Significantly upgraded silicon detector: 

- 707,000 channels compared with 46,000 in Run I 

- Axial, stereo, and 90° strip readout 

- Full coverage over the luminous region along the beam axis 

- Radial coverage from 1.35 to 28 cm for |7y| < 2 

- Innermost silicon layer ("LOO") on the beampipe with 6 /im axial hit resolution 

• Outer drift chamber capable of 132 ns maximum drift 

- 30,240 sense wires, 44-132 cm radius, 96 dE/dx samples possible per track 

• Fast scintillator-based calorimetry out to \r]\ ~ 3 

• Expanded muon coverage 

• Improved trigger capabilities 

- Drift chamber tracks with high precision at Level- 1 

- Silicon tracks for detached vertex triggers at Level-2 

• Expanded particle identification via time-of-flight and dE/dx 

2.2.1 CDF Coordinate System 

The CDF detector uses a right-handed coordinate system. The horizontal direction pointing out 
of the ring of the Tevatron is the positive x-axis. The vertical direction pointing upwards is the 
positive y-axis. The proton beam direction, pointing to the east, is the positive 2;-axis. 

A spherical coordinate system is also used. The radius r is measured from the center of the 
beamline. The polar angle 9 is taken from the positive 2;-axis. The azimuthal angle (p is taken 
counter-clockwise from the positive x-axis. 

At a pp collider, the production of any process starts from a parton-parton interaction which 
has an unknown boost along the 2;-axis, but no significant momentum in the plane perpendicular 
to the z-axis, i.e. the transverse plane. This makes the transverse plane an important plane in app 
collision. Momentum conservation requires the vector sum of the transverse energy and momentum 
of all of the final particles to be zero. The transverse energy Et and transverse momentum px for 
a particle produced in a pp collision are defined by Et = £^ x sin 9 and px = p x sin 6*. We use 
the convention that "momentum" refers to pc and "mass" to mc^. 

Another quantity invariant under Lorentz boosts along the beamline is Rapidity, which is 
defined as |/ = ^log{^z^)i where Pl is the longitudinal momentum along the beamline and E is 
the energy. 

Pseudorapidity rj is used by high energy physicists and is defined as ?7 = — In tan |. For massless 
particles rj = y. 
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LAYER 00 SVX II INTERMEDIATE SILICON LAYERS 

Figure 2.3: CDF II tracking volume. 

Hard pp head-on collisions produce significant momentum in the transverse plane. The CDF 
detector has been optimized to measure these events. Typically, particles in a pp collision event 
tend to be more in the forward and backward regions than in the central region because there is 
usually a boost along the 2;-axis. The derivative of r] is dr] = 

A constant rj slice corresponds to variant 6 slice which is smaller in the forward and backward 
regions than in the central region. This makes the t] occupancy more uniform than 6 occupancy. 
Therefore, for example, calorimeters are constructed in rj slices instead of 6 slices. 

2.2.2 Tracking 

The CDF detector features excellent charged particle tracking and good electron and muon identifi- 
cation in the central region. The detector is built around a 3 m diameter, 5 m long superconducting 
solenoid operated at 1.4 T. The tracking volume is surrounded by the solenoid magnet and the 
endplug calorimeters as shown in Figure [231 

The CDF tracking system includes a central outer drift chamber (COT) and the silicon tracker. 
The main parameters of the CDF tracking system are summarized in Table 12. 1[ 

The Silicon Tracker 

Enhanced system of silicon microstrip detectors ^38j consists of three components: Layer 00, 
the Silicon VerteX detector II (SVX II), and the Intermediate Silicon Layers (ISL). r — view of 
the silicon tracker is shown in Figure 12.41 

A single layer rad-hard Layer 00 detector is mounted on and supported by the beam pipe. The 
Layer 00 single-sided sensors provide 6 /xm axial hit resolution. 
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COT 


Radial coverage 


44 to 132 cm 


Number of superlayers 


8 


Measurements per superlayer 


12 


Maximum drift distance 


0.88 cm 


Resolution per measurement 


180 /im 


Rapidity coverage 


\v\ < 1-0 


Number of channels 


30,240 


Layer 00 


Radial coverage 


1.35 to 1.65 cm 


Resolution per measurement 


6 /im (axial) 


Number of channels 


13,824 


SVX II 


Radial coverage 


2.4 to 10.7 cm, staggered quadrants 


Number of layers 


5 


Resolution per measurement 


12 /im (axial) 


Total length 


96.0 cm 


Rapidity coverage 


\v\ < 2.0 


Number of channels 


423,900 


ISL 


Radial coverage 


20 to 28 cm 


Number of layers 


one for < 1; two for 1 < 77 < 2 


Resolution per measurement 


16 /im (axial) 


Total length 


174 cm 


Rapidity coverage 


\v\ < 1-9 


Number of channels 


268,800 



Table 2.1: Design parameters of the CDF tracking systems 

The next five layers compose the SVX II and are double-sided detectors. The axial side of each 
layer is used for r-0 measurements. The stereo side of each layer is used for r-z measurements. 

The two outer layers compose the ISL and are double-sided detectors. This entire system allows 
charged particle track reconstruction in three dimensions. The impact parameter resolution of 
SVX II + ISL is 40 /im including 30 /im contribution from the beamline. The zq resolution of 
SVX II + ISL is 70 /im. The main parameters of the silicon tracker are summarized in Table 12.11 

The Central Outer Tracker (COT) 

The COT [37] is a multi-wire open-cell drift chamber for charged particle reconstruction, oc- 
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cupying the radial region from 44 to 132 cm and \z\ <155 cm. The COT replaced the Central 
Tracking Chamber (CTC), which, in addition to aging problems observed during Run I, would 
also suffer from degraded performance at £ > 1 x lO^^cm^^s^^. The major problem with the 
CTC would be its maximum drift time (800 ns) relative to the expected bunch crossing time in 
Run II (132 ns). 

To address this, the COT uses small drift cells (~2 cm wide, a factor of 4 smaller than the CTC) 
and a fast gas to limit drift times to less than 130 ns. Each cell consists of 12 sense wires oriented 
in a plane, tilted with respect to the radial (Figure [23]) • A group of such cells at a given radius is 
called a superlayer. There are eight alternating superlayers (Figure 12.51) of stereo (nominal angle 
of 2°, used for r-z measurement) and axial (used for r-0 measurement) wire planes. The main 
parameters of the COT are summarized in Table 12.11 

The COT is filled with a mixture of Argon:Ethane = 50:50 which determines the drift velocity. 
A charged particle travels through the gas mixture and produces ionization electrons. The elec- 
trons drift toward the sense wires in the electric field created by cathode field panels and potential 
wires of the cell. In the crossed magnetic and electric fields electrons originally at rest move in 
the plane perpendicular to the magnetic field at an angle a with respect to the electric field lines. 
The value of a depends on the magnitude of both magnetic and electric fields and the properties 
of the gas mixture. In the COT a ~ 35°. 

The optimal situation for the resolution is when the drift direction is perpendicular to that of 
the track, which is true for high px tracks because they are almost radial. To make the ionization 
electrons drift in the (j) direction all COT cells are tilted by 35° with respect to the radial. 
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When an electron gets near a sense wire the local 1/r field accelerates them causing more 
ionization and thus forms an "avalanche" producing a signal (hit) on the sense wire. By measuring 
the arrival time of "first" electrons (drift time) at sense wire relative to collision time, the distance 
of the hit is calculated. 

COT tracks above 1.5 GeV are available for triggering at Level 1 (XFT tracks). SVX layers 
0-3 are combined with XFT tracks at Level-2 (SVT tracks) (see Section 12.2.71 for the description 
of CDF Run II trigger system). 

2.2.3 Calorimetry 

The energy measurement is done by sampling calorimeters, which are absorber and sampling scin- 
tillator sandwiches with phototube readout. Outside the solenoid, Pb-scintillator electromagnetic 
(EM) and Fe-scintillator hadronic (HAD) calorimeters cover the range |?7| < 3.6. The central 
calorimeter systems have been retained from Run I, but the plug calorimeters are new detectors 
for Run II. 

Both the central {\ri\ < 1.1) and plug (1.1 < \ri\ < 3.6) electromagnetic calorimeters have 
fine grained shower profile detectors at electron shower maximum, and preshower pulse height 
detectors at approximately IXo depth. Electron identification is accomplished using E/p from 
the EM calorimeter; using HAD /EM ~ 0; and using shower shape and position matching in the 
shower max detectors. The calorimeter cell segmentation is summarized in Table 12.21 and shown 
in Figure 13.51 A comparison of the central and plug calorimeters is given in Table 12.31 



\ri\ Range 


A0 


At] 


0. - 1.1 (1.2 h) 


15° 


~ 0.1 


1.1 (1.2 h) - 1.8 


7.5° 


~ 0.1 


1.8 - 2.1 


7.5° 


~ 0.16 


2.1 - 3.64 


15° 


0.2 - 0.6 



Table 2.2: CDF II Calorimeter Segmentation 

The region 0.77 < rj < 1.0, 75 ° < < 90° is uninstrumented to allow for cryogenic utilities 
servicing the solenoid. 

Any high energy electron or photon passing through the electromagnetic calorimeters, will 
undergo pair production (7 — > e+e~) and bremsstrahlung (e^ 7^^) thus producing an electro- 
magnetic shower. The point at which the electromagnetic shower consists of the largest amount of 
particles is known as the shower maximum. At this point the average energy per particle becomes 
low enough to prevent further multiplication. After the shower maximum, the shower decays 
slowly through ionization losses for the electrons and positrons or by Compton scattering for the 
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Central 


Plug 


EM: 






Thickness 


19Xo, lA 


21X0, lA 


Sample (Pb) 


0.6X0 


0.8X0 


Sample (scint.) 


5 mm 


4.5 mm 


Stoch. res. 


14%//E^ 


16%/v^ 


Shower Max. seg. (cm) 


1.40X (1.6-2.0) Z 


0.5 X 0.5 UV 


Hadron: 






Thickness 


4.5A 


7A 


Sample (Fe) 


1 to 2 in. 


2 in. 


Sample (scint.) 


10 mm 


6 mm 



Table 2.3: Central and Plug Upgraded Calorimeter Comparison 



photons. The calorimeters measure the energy deposited by these showers, and hence the energy 
of the incident particle. Electromagnetic calorimeters are designed to fully contain showers from 
electrons and photons. 

Hadrons lose energy by nuclear interaction cascades which can have pions, protons, kaons, 
neutrons, neutrinos, muons, photons, etc. This is significantly more complicated than an electro- 
magnetic cascade and thus results in a large fluctuation in energy measurement. 

Central Calorimeters 

The central calorimeters consist of the central electromagnetic calorimeter (CEM) [39], the central 
hadronic calorimeter (CHA) [ID], and the end wall hadronic calorimeter (WHA). 

The CEM and CHA are constructed in wedges which span 15° in azimuth and extend about 
250 cm in the positive and negative z direction, shown in Figure [2l6l There are thus 24 wedges on 
both the and —z sides of the detector, for a total of 48. A wedge contains ten towers, each of 
which covers a range 0.11 in pseudorapidity. Thus each tower subtends 0.11 x 15° in 77 x 0. The 
CEM covers < jr/l < 1.1, the CHA covers < |?7| < 0.9, and the WHA covers 0.7 < < 1.3. 

The CEM uses lead sheets interspersed with polysterene scintillator as the active medium 
and employs phototube readout, approximately 19Xo in depth, and has an energy resolution of 
13.5%/v^©27c0. 

To provide more accurate information on the position of the electromagnetic shower inside the 
calorimeter, the Central Electromagnetic Shower (CES) [39] detector is embedded inside the CEM 
at the shower maximum, at a depth of approximately 6 radiation lengths. The CES detector is 

^© denotes addition in quadrature 
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a proportional strip and wire chamber situated at a radius of 184 cm from the beamhne. In the 
azimuthal direction, cathode strips are used to provide the z position and in the direction, anode 
wires are used. These wires can effectively measure the transverse shower profile to distinguish 
between a single shower from a prompt photon and two showers from a decay of a neutral meson 
to two photons, e.g. 7r° 77, with a position resolution of 2 mm at 50 GeV. 

In order to help particle identification, specifically between electromagnetic and hadronic show- 
ers the central preradiator detector (CPR) is mounted on the front of the calorimeter wedges, at a 
radius of 168 cm from the beamline, and uses the solenoid and tracking detectors as a radiator. It 
uses proportional chambers to sample the early development of the shower to measure conversions 
in the coil, helping to distinguish prompt photons and electrons from photons originating from 
tt" decays and electrons from conversions. A prompt photon has a 60% probability of converting, 
while the conversion probability of at least one photon from vr" — > 77 is about 80% [H]. 

The CHA and WHA use steel absorber interspersed with acrylic scintillator as the active 
medium. They are approximately 4.5A in depth, and have an energy resolution of 75%/-\/Et©3%, 
as measured on the test beam for single pions [IQ]. 

Plug Calorimeters 

The plug calorimeters consist of the plug electromagnetic calorimeter (PEM) [35], and the plug 
hadronic calorimeter (PHA). At approximately 6X0 in depth in PEM is the plug shower maximum 
detector (PES). Figure [2781 shows the layout of the detector and coverage in polar angle 36.8° > 
^ > 3° (1.1 < \ri\ < 3.64). Each plug wedge spans 15° in azimuth, however in the range 36.8° > 
6 > 13.8° (1.1 < < 2.11) the segmentation in azimuth is doubled and each tower spans only 
7.5°. 

The PEM is a lead-scintillator sampling calorimeter. It is approximately 21Xo in depth, and has 
an energy resolution of 16%/\/E'©l%. The PES consists of two layers of scintillating strips: U and 
V layers offset from the radial direction by +22.5° and —22.5° respectively, as shown in Figure [2791 
The position resolution of the PES is about 1 mm. The PHA is a steel-scintillator sampling 
calorimeter. It is approximately 7A in depth, and has an energy resolution of 74:%/\/E © 4%, as 
measured on the test beam for single pions [30] . 

2.2.4 Muon Systems 

The muon is a minimum ionizing particle which loses very little energy in detector materials. The 
muon's lifetime, 2.2 /is, is long enough for the muon to pass through all the detector components, 
reach the muon chambers and decay outside. 

There are four independent muon systems: the central muon detector (CMU) [13], the central 
muon upgrade (CMP) [33], the central muon extension (CMX) [45j, and the intermediate muon 



33 



Figure 2.6: CEM/CES/CHA wedge. 



Figure 2.7: CES strip and wire. 




detector (IMU). The calorimeter steel serves as a filter for muon detection in the CMU and CMX, 
over the range \ri\ < 1, pt > 1.4 GeV. Additional iron shielding, including the magnet yoke, 
provides a muon filter for the CMP in the range \r]\ < 0.6, pt > 2.2 GeV. The (non-energized) 
forward toroids from Run I provide muon filters for IMU in the range 1.0 < |?7| < 1.5 for pj- > 2 
GeV. Scintillators for triggering are included in CMP, CMX, and IMU. 

Muon identification is accomplished by matching track segments in the muon chambers with 
COT/SVX tracks; matching is available in r — (p for all detectors and in z in CMU and CMX. 
The parameters for the muon systems are summarized in Table 12.41 The IMU, which provides 
coverage in the forward region will not be discussed in detail, as it is not used for this analysis. 

The coverage for the muon systems in rj — (p space is shown in Figure 12. Ill CMU, CMP and 
CMX muon systems are also shown in Figure 13.61 
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CMU CMP 


CMX 


IMU 


Pseudo-rapidity coverage 


|r/|<0.6 |r/|<0.6 0.6 < |r/| < LO 


1.0 < \r] 


< 1.5 


Drift tubes 


Cross-section, cm 


2.68 X 6.35 2.5 x 15 


2.5 X 15 


2.5 X 


8.4 


Length, cm 


226 640 


180 


363 


Max drift time, fis 


0.8 1.4 


1.4 


0.^ 




Scintillation counters 


Thickness, cm 


2.5 


1.5 


2.5 


Width , cm 


30 


30-40 


17 


Length , cm 


320 


180 


180 


Minimum muon px, GeV 


1.4 2.2 


1.4 


1.4-2.0 



Table 2.4: Design Parameters of the CDF II Muon Detectors. 



A muon chamber contains a stacked array of drift tubes and operates with a gas mixture of 
Argon:Ethane = 50:50. The basic drift principle is the same as that of the COT, but the COT is 
a multi-wire chamber, whereas at the center of a muon drift tube there is only a single sense wire. 
The sense wire is connected to a positive high voltage (HV) while the wall of the tube is connected 
to a negative HV to produce a roughly uniform time-to-distance relationship throughout the tube. 
The drift time of a single hit gives the distance to the sense wire, and the charge division at each 
end of a sense wire can in principle be used to measure the longitudinal coordinate along the sense 
wire. The hits in the muon chamber are linked together to form a short track segment called a 
muon stub (Figure [2.101) . If a muon stub is matched to an extrapolated track in the tracking 
system (Section |2. 2. 2p . a muon is reconstructed. 

CMU and CMP 

The CMU is unchanged from Run I. It is located behind the towers of the CHA and divided into 
wedges covering 12.6° in azimuth for rj < 0.6. Only muons with a px > 1.4 GeV reach the CMU. 
Each wedge has three towers, each comprised of four layers of four drift tubes. The second and 
fourth layers are offset by 2 mm in direction from the first and third as shown in Figure [2.101 Six 
CMU wedges and their relative location with respect to CMX and CMP (outer box) subsystems 
are shown in Figure 12.121 

A 50 yum diameter stainless steel resistive sense wire is located in the center of each cell. The 
wires in the cells in the first and third (second and fourth) layers are connected in the readout. 
Each wire pair is instrumented with a time-to-digital converter (TDC) to measure the 0-position 
of the muon and an analogue-to-digital converter (ADC) on each end to measure z position via 
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Muon track 



Radial centerline 



55 mm 



^ To pp interaction vertex 

Figure 2.10: A CMU module in the r-0 plane with 4 layers of drift chambers. The drift times 
tl and t2 are used to calculate muon momentum for triggering (Section 12.2.71) . The sense wires 
connected to the readout are shown as a black circles. 



m - CMX 
-1 



3- CMP SS - CMU IMU 
1 




Figure 2.11: Location of the muon detec- 
tors in and rj. On the east side, there is 
a gap in coverage in the CMX of 30° in az- 
imuth, due to the location of the cryogenic 
utilities servicing the solenoid. 




Figure 2.12: Central muon systems: CMU, 
CMP, CMX. 



charge division. The position resolution of the detector is 250 /im in the drift direction (r-c/)) and 
1.2 mm in the sense wire direction (z). 
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Approximately 0.5% of high energy hadrons produced will pass through the CMU creating an 
irreducible fake muon background. In order to reduce this effect, an additional muon chamber 
(CMP) is installed behind 60 cm of steel. 

The CMP consists of a four-sided box placed on the outside of the CDF detector. Muons with 
Pt > 2.2 GeV can reach the CMP. The rectangular form of the CMP detector means that its rj 
varies in azimuth (Figure [2.111) . The CMP covers \ri\ < 0.6. 

The maximum drift time of the CMU is longer than the pp bunch crossing separation, which 
can cause an ambiguity in the Level 1 trigger (Section l2.2.7l) . To resolve the ambiguity scintillation 
counters are used. The scintillation counters (CSP) are installed on the outer surface of the CMP. 

Central Muon Extension (CMX) 

The CMX has eight layers and extends the r] coverage to 0.6 < \t]\ < 1.0. It consists of two 120° 
arches located at each end of the central detector, as shown in Figure 12.111 The uninstrumented 
regions have been filled by the insertion of a 30° keystone at the top, and two 90° miniskirts for 
the lower gaps. There is a gap in the coverage on the east side due to cryogenic utilities servicing 
the solenoid as shown in Figure I2.1H known as the " chimney" . 

A layer of scintillation counters (the CSX) is installed on both the inside and the outside surfaces 
of the CMX. No additional steel was added for this detector because the large angle through the 
hadron calorimeter, magnet yoke, and steel of the detector end support structure provides more 
absorber material than in the central muon detectors. 

2.2.5 Time of Flight System 

The Time of Flight detector (TOF) |16] measures the time taken by a particle to travel from the 
interaction point to the detector, and has a particle timing resolution of 100 ps. This information 
can be used to differentiate between different particle types (e.g. kaons, pions) and also to help 
tag cosmic ray events. 

The TOF is situated between the COT and the solenoid. It consists of 216 scintillator bars with 
dimensions 4x4x276 cm. At each end of the scintillator bars a photomultiplier tube is mounted 

2.2.6 Cherenkov Luminosity Counters 

Luminosity (£) is a measure of particle interaction, specifically the chance that a proton will 
collide with an antiproton. The rate of inelastic scattering in pp interactions can be used to 
determine the C . 

A gas Cherenkov Luminosity Counter (CLC) [47] measures the number of interactions per beam 
crossing to determine the luminosity of the Tevatron. There are two CLCs positioned between 
the beam-pipe and the plug calorimeters, covering the region 3.7 < \ri\ < 4.7. Each CLC consists 
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of 48 thin, conical gas-filled Cherenkov counters. They are arranged in three concentric circles, 
each consisting of 16 counters (see, for example. Figure [STT]) • 

The luminosity of a pp collider can be estimated using the equation: 



where f is the frequency of bunch crossing, /z is the average number of interactions per beam 
crossing, given by the CLC hit rate (about 5-6), and a is the inelastic cross-section of pp scattering. 
The average of the inelastic cross-sections as measured by CDF Run I and the E811 [48j is 60.7±2.3 
mb Hg. 

A total systematic uncertainty of 6% is quoted for all luminosity measurements. This includes 
a 4.4% contribution from the acceptance and operation of the luminosity monitor and 4.0% from 
the theoretical uncertainty on the calculation of the total pp cross-section [19]. 



2.2.7 Trigger and Data Acquisition 

Many interesting physics processes have cross sections which are many orders of magnitude smaller 
than the total inelastic cross section. The collision rate at the Tevatron is much higher than the 
rate at which data can be recorded. Therefore, the trigger needs to be fast and accurate to record 
as many interesting events as possible, while rejecting uninteresting events. 

To accomplish this, the CDF trigger system has a three-level architecture: Level 1 (LI), Level 
2 (L2), and Level 3 (L3). The data volume is reduced at each level, which allows more refined 
filtering at subsequent levels with minimal deadtime. Each sub-detector generates primitives 
which can be used in the trigger system to select events. The trigger system block diagram is 
shown in Figure [2.131 

At LI axial layers of the COT are used by eXtreme Fast Tracker (XFT) to reconstruct </> and px 
for the tracks. Based on the XFT tracks and a ratio of the hadronic energy to the electromagnetic 
energy of a calorimeter tower (HAD/EM ratio) electrons and photons are then reconstructed. 
Muons are reconstructed by matching XFT and muon hits. Jets are reconstructed based on a 
sum of the electromagnetic and hadronic energies for a tower, and ^ Et (a scalar sum of the 
energies of all of the calorimeter towers) are also reconstructed at LI. 

LI is a synchronous hardware trigger and it makes a decision within 4 /is. This trigger reduces 
the event rate from 7.6 MHz to 50 KHz, which is limited by the L2 processing time. Accepted 
events are then passed to the L2 hardware. 

At L2 SVX layers 0-3 are combined with XFT tracks (Figure r2. 131 ). The L2 uses jet clustering 
as well as improved momentum resolution for tracks, finer angular matching between muon stubs 
and central tracks and data from the CES for improved identification of electrons and photons. 

The L2 decision time is about 20 yus. L2 is a combination of hardware and software triggers 
and is asynchronous. If an event is accepted by LI, the front-end electronics moves the data to 
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Figure 2.13: Trigger system block diagram: flow of data through LI and L2 of trigger system. 
Silicon Vertex Tracker (SVT), based on the track impact parameter of displaced tracks, is used 
by L2. The data volume is reduced at each level which allows more refined filtering at subsequent 
levels with minimal deadtime. L3 is purely a software trigger consisting of the event builder 
running on a large PC farm. 



one of the four onboard L2 buffers. This is sufficient to process a LI 50 KHz accept rate and to 
average out the rate fluctuations. The L2 accept rate is about 300 Hz which is limited by the 
speed of the event builder in L3. 

L3 is a purely software trigger consisting of the event builder running on a large PC farm. 
Data which passes one of the specified trigger paths is reconstructed at L3 using full detector 
information and the latest calibrations. The L3 accept rate of about 75 Hz is limited by the speed 
of writing data to tape for the permanent storage. 

As soon as an event passes L3 it is delivered to the data-logger subsystem which sends the 
event out to permanent storage for offline reprocess, and to online monitors which verify that the 
entire detector and trigger systems are working properly. One of the online monitors, CDF Run 
n Event Display, is described in detail in Section [31 
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Chapter 3 

CDF Run II Event Display 



Event displays are indispensable tools for data analysis in high energy physics experiments. They 
help to understand the physics of a recorded interaction, to diagnose the apparatus, to make the 
detector geometry imaginable and to illustrate the whole matter to general audiences. 

CDF Run 11 Event Display (EVD) [llj is a major contribution to the commissioning and 
operation of CDF 11 detector. Development and support of the EVD package [lO] is one of the 
responsibilities of the ITEP (Moscow) group at CDF |13j . 

The data from the collider experiment is a stream of signals from subdetectors. The detector 
"sees" these signals as sequences of impulses, distributed over many channels of different subde- 
tectors (see Section 12^21) . The signals are analyzed by powerful pattern recognition and analysis 
programs, which create more sophisticated objects like clusters, segments, stubs and then tracks, 
muons, electrons/photons, jets, etc. 

Typically physics results are based on statistical analysis of many events. The standard forms 
of presentation are histograms, graphs and tables. However, a graphical representation of a single 
event a powerful tool for checking the validity of reconstruction or analysis algorithms. For a 
quick assessments of error conditions as well as for public presentations the visual representation 
is the most efficient way to transfer information to a human brain. 

Higher event multiplicities and higher momenta of outgoing particles are matched by detectors 
with a growing number of subunits of increasing granularity, resolution and precision. As a 
consequence pictures of detectors and events are getting more complicated and may even get 
incomprehensible. This leads to a question if the presentation of data via visual techniques is 
useful for complicated events at the Tevatron. 

3.1 Overview 

The aim of the EVD is to enable visual representation of the objects existing in the CDF Run 11 
software. EVD interacts both with the data and with the simulation and reconstruction packages. 
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For simulated data EVD visualizes the Decay Tree, which is constructed from HEPG information. 

It is natural to define three kinds of objects: Real Objects (Section I3.1.ip . Graphical Objects 
(Section I3.1.2p . and Views (Section I3.1.3p . To visualize Real Objects and to access information 
about the event, Operations (Section I3.1.4p are performed on the Graphical Objects and Views. 

3.1.1 Real Objects 

A Real Object is information from a subdetector or combined information from several subdetec- 
tors. For example, to identify electron one needs calorimeter and tracking information. Some of 
the Real Objects in CDF are listed below: 

• Tracking information (Section 12.2.20 

- axial and stereo hits from the COT 

- hits from the silicon tracker 

- tracks reconstructed with different tracking algoritms 

• Calorimeter information (Section I2.2.3P 

- Central and Plug Shower Chambers information 

- Central and Plug Preradiator information 

• Hits from TOF system (Section I2.2.5P 

• Information from the muon systems (Section I2.2.4p 

- Hits from the muon systems 

- Track segments reconstructed in the muon chambers 

• Information from East and West CLC subdetectors, which are used to monitor luminosity 
and to identify diffractive events (Section I2.2.6P 

• Information from Beam Shower Counters [50], used to identify diffractive events 

• Information from CDF Run II L1/L2/L3 Trigger System (Section I2.2.7P 

• EM Timing information [51] 

• Pre-reconstructed objects, such as Muon, Electron, Photon, Jet Candidates, 

• Full information about raw and reconstructed data 

- access to banks and collections in the event 

A major requirement has been made to keep analysis in the EVD to a minimum, as objects 
are identified differently in different analyses (i.e. loose electron for one analysis may be a jet in 
some other analysis). However, EVD has a functionality to clean complicated events by selective 
presentation of parts of the data and the detector. (Section 13.1.40 . For example, EVD can show 
or hide tracks depending on their properties (0, rj, Et, Pt, number of hits in a subsystem etc.). 
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3.1.2 Graphical Objects 

Real Objects in the EVD are mapped to their visual representation into Graphical Objects of 
different types, corresponding to different Views. The graphical objects correspond to the stored 
real objects (for example, to the list of hits and tracks) and other real objects (for example, to 
electron, muon and photon candidates) created from them. 

The properties of Real Objects are used to display Graphical Objects. For instance, information 
from the calorimeter is shown as towers with a size proportional to the deposited energy. We use 
graphics libraries available in the ROOT package [52] . 

3.1.3 Views 

View is a method of visualizing a set of graphical objects. For a user, view is generally a window 
with a defined way of displaying Graphical Objects. We define three categories of views for the 
CDF Run II Event Display: 

1. Realistic Views are obtained by either a sequence of rotations, linear scaling, and projec- 
tions of a geometry of detector /identified objects. Reahstic views are understood intuitively, 
although the pictures may become too complex (for example, see Figures 13^ and 13^ . 

Hits density as well as detector precision grow towards the interaction point. Therefore, 
an ability to hide parts of the detector obscuring the picture is crucial. Another feature of 
the EVD is to show or hide Graphical objects depending on the properties of Real Objects 
(number of hits, drift time etc.) 

2. Schematic Views can be obtained from realistic views by changing the aspect ratio or 
focal length for a subdetector. For example, the scale may be decreased with increasing 
radius, so that the inner subdetectors appear enlarged. This emphasizes the commonly used 
construction principle of detectors, namely that precision and sampling distance decrease, 
when stepping from the inner to the outer detectors. 

Schematic Views do not necessarily represent the detector in its real proportions or shapes. 
For example a box that changes color depending on error conditions may serve as a schematic 
view of any detector component. In many cases the schematic views are relatively easy to 
understand and very efficient to use (for example, see Figure 1331) . 

3. Abstract Views have little resemblance to the detector image in cartesian space. These 
views are not intuitively understood but given some training may be the most powerful. 
Examples of abstract views are angular projections, histograms, LEGO plots, etc. (for 
example, see Figure [375]) 
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3.1.4 Operations 



Operations on Graphical Objects: 

These operations only change the visual properties of the graphical objects, the real objects are 
not changed. Visual operations can be performed on the graphical objects or on the views. The 
following operations on graphical objects are supported: 

• rotation, zooming, translation, scaling 

• viewpoint changing 

• object hiding/unhiding 

• sub-view creation 

• redefinition of visual properties of graphical objects 
Operations on Real Objects: 

Operations on real objects can access information about the objects or change the state of these 
objects. In addition, visual properties of the corresponding graphics objects can also be changed 
as a result of changing real objects. The following operations on real objects are supported: 

• access to public member functions for an object from the EVD 

• obtaining the detailed information on an object properties (for example, for tracks one can 
access information on number/type of hits in COT/Silicon detectors, dE/dx, x^, do, zq etc.) 

• access to event record information (data banks and collections) 

• application of identification cuts 

• interface to a histogramming package 



3.2 CDF Run II Application Framework 

The design of the EVD is based on the belief that both requirements and graphics software abilities 
will be very broad at any time and will constantly evolve. The EVD is designed to accommodate 
that diversity and change. This can be accomplished only by sufficient flexibility and modularity 
of the core control structure [53] . 

An application framework, in the context of a high energy physics experiment, is a "system" 
that allows physicist-developed code to be combined with code developed by other people and 
to be used in both the online and offiine environments. Either real or simulated data can be 
used as the input and the output can include (modified) copy of the input as well as additional 
reconstructed quantities. This output then forms the input in the next stage of a multi-stage data 
reduction environment. CDF Run II Application Framework (AC++) is based on the ROOT 
object oriented analysis framework [52], which is incorporated in CDF Run II C++ software. 
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The AC++ provides a unified framework for event reconstruction, post-production analysis and 
online monitoring, triggering and calibration. The goals of the framework is to provide a simple 
and straightforward means of combining any number of independent classes, called modules, into 
a single executable and to provide a flexible system for specifying (either interactively or in a 
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(b) The EVD control panel with the Loop tab (c) Data Menu to select which data to read, 
selected. 

Figure 3.1: CDF Run II Event Display Graphical User Interface. The start of the EVD is signalled 
by the appearance of the manager window (a). Clicking the Panel button brings up the EVD 
control panel (b). The control panel allows you to configure the behavior of the EVD. For example, 
user may select which data to read using the DataMenu (c). 
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batch mode) how these modules are run. 

Therefore, EVD has the flexibihty to enable/disable different modules, so that one can work 
with unprocessed data for the immediate feedback. Alternatively, one can also run reconstruction 
modules inside the EVD with user-defined parameters for debugging purposes. 



3.3 Graphical User Interface 



The start of the EVD is signalled by the appearance of the manager window (Figure [3?T| a). From 
the manager window user select one of the event displays, and control automatic looping through 
events, as is done in the CDF Run II control room at B0. The name of the input file is displayed 
at the bottom. 

Clicking the Panel button brings up the EVD control panel (Figure ISTTl b). The control panel 
allows user to configure the behavior of the EVD. User may select which data to read using the 
DataMenu (Figure \3A] c) 

Complicated events might be cleaned by selective presentation of parts of the data and the 
detector. User can specify cuts on physical quantities of the displayed objects, such as the r], 0, 
Pt of tracks or the energies deposited in the calorimeter towers. This helps to clean up the event 
by removing low-px tracks and 1ow-Et towers, or to debug reconstruction problems by requiring 
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Figure 3.2: Cut Manager. Complicated events might be cleaned by selective presentation of the 
data. User can specify cuts on physical quantities of the displayed objects, such as the rj, 0, px of 
tracks or the energies deposited in the calorimeter towers. 
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EVD to show objects which pass some cuts (e.g. number of COT or Sihcon hits, or tracks with 
large do etc.). Figure [3^ shows sample Cut Manager session. 



3.4 Displays 

3.4.1 r — (f) and r — z Views 

In layered projections, each geometry object acquires a 2-dimensional shape which can be different 
in each projection. For example, the drift chamber outline is a circle in the r — cf) view and a 
rectangle in the r — z view. 

The r — view (COT Display, Figure [373]) shows hits in the subdetectors (COT, Silicon trackers, 
Muon Chambers, TOF, XTRP etc.), energies in the central EM and HAD calorimeter towers 
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Figure 3.3: COT Display (r — view). 
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(summed over rj), missing Et information, information about CDF reconstructed objects, such as 
e, fi, 7, jet candidates. It also gives the information about run/event number, as well as trigger 
information. 

Figure [33] is an example of a schematic view (Section 13. 1.3p . and it is obtained from a realistic 
r — (f) projection by applying a fish-eye transformation to a COT volume. 

A fish-eye view introduces a nonlinear transformation of radius in the layered r — (p projection, 
with the aim to enable simultaneous inspection of tracking chambers, calorimetry and muon system 
within the same picture (Figure 13731) . 
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Figure 3.4: RZ Display {r ~ z view) 



The r — z view (Figure 13. 4p is designed to show the same information as r — view, but in 
r — z projection. The only exception are COT hits, which do not have z coordinate and therefore 
are displayed in r — projection only. Energies in the central EM and HAD calorimeter towers 
are summed over (p. Figure [3741 is an example of a realistic view (Section 13. 1.3p . 
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User can "slice" the detector by selectin pair of opposite wedges ("Select Slice" option at the 
bottom of the RZ display menu, Figure 13.41) to reduce the amount of information displayed. In 
this case one will have only those wedges' tracks, hits, bits, and calorimeter towers. There are 24 
15° wedges, which are displayed in 12 opposite pairs. The default slice value is -1, which folds all 
upper wedges onto the top and all lower wedges onto the bottom, which is the way the RZ display 
used to work in CDF Run I. 

3.4.2 Lego Displays 

These are generic LEGO plot windows (Figure [3.51) showing a variable (E, Et, ADC counts) as 
a function of eta-phi. Figure [375] is an example of an abstract view (Section I3.1.3p . rj — (p grid 
corresponds to 77 — segmentation of CDF calorimeter system (Section [223]); each bin on the lego 
display represents a tower. Towers corresponding to CDF Run II particle candidates (e, /x, 7, jets) 
have been added to LEGO views to improve the display (see Figure [331) . User can interactively 
rotate LEGO display/change default settings to obtain a better view of an event. 

There also many other LEGO-based views, used to monitor/debug specific subdetectors, such 
as the PLUG/PPR Views to visualize Plug Calorimeter information, and the CES View to show 
the sums of energies for CES Strips and Wires. In addition to the views itself there is an interface 
to a histogramming package to obtain specific distributions. For example, one can see distributions 
of measurements in strips/wires of the CES. Views are designed to give a general idea about an 
event, and histograms allow to see more detailed picture for a part of a subdetector. 

3.4.3 3D Displays 

Perspective 3-dimensional views with hidden lines and hidden surface removal are very useful for 
understanding detector geometry, and provide attractive pictures for public relation (PR) purposes 
(Figure [XB]) . Analyzing the event itself is often less successful in this mode, since the complicated 
geometry tends to obscure the tracks and hits. 

Figure [376[ is an example of a realistic view (Section l3.1.3l) . There are several available 3D views 
which one can use separately for specific needs or combine to obtain a complicated view. 

• CDF 3D display 

- 3-dimensional view of the CDF detector with tracks/silicon hits/muon hits and stubs) 

• 3D calorimeter display 

- central and plug calorimeter towers together with tracks 

• SVX 3D display 

- dedicated silicon detector display, which shows the silicon hits/strips together with tracks). 
SVX 3D display is designed to obtain detailed information for all silicon hits associated with 
a given track, which is helpful for debugging purposes. 
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Three-dimensional Views with hidden hnes and surface removal are possible through the special 
OpenGL viewer which is integrated in view. Figure [3l6] has been obtained using OpenGL. In this 
mode live rotations are possible given a suitable hardware acceleration for the instantaneous 
response. 



3.4.4 Other Displays 

There many other views (see Figure 13771) . used to monitor different CDF Run II subsystems and 
to analyze real and simulated data: 

• Wedge Display 

- CES/CPR information together with central calorimeter information 
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Figure 3.5: The Lego Display shows the variable Et as a function of eta-phi. 
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Figure 3.6: 3-dimensional OpenGL view of detector geometry with hidden hnes removal. 

• Trigger Display 

- Level 1, Level 2, Level 3 trigger bits and corresponding names 

• BSC Display 

- shows ADC counts from beam shower counter subdetectors 

• CLC Display 

- shows ADC counts from Cherenkov Luminosity Counters 

• MC Decay Tree Display 

- decay tree constructed from HEPG information for the simulated data 

3.5 Live Events 

The Live Events page has been designed to provide attractive pictures for PR purposes [12]. The 
Views displayed on Figure 13.71 from the top to the bottom are as follows: 

• 1st row: 

COT Display. COT Display zoomed to access central outer tracker information. COT 
Display zoomed to access silicon tracker information 
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Figure 3.7: CDF "Live Events" public page. 



• 2nd row: 

RZ Display. SVX 3D Display. Calorimeter LEGO Display 

• 3rd row: 

East Plug LEGO Display. Calorimeter 3D Display. West Plug LEGO Display 
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• 4th row: 

East CLC Display. BSC Display. West CLC Display 

3.6 Conclusions 

Prom the very beginning of Run II data taking EVD is in continuous use for the online moni- 
toring and for the analyses. This answers the question whether a presentation of data via visual 
techniques is possible for complicated events at the Tevatron. More importantly, features of the 
EVD make it to be one of the most important tools for a better understanding of the events which 
could possibly be New Physics candidates in a CDP Run II data. 
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Chapter 4 

+ X Selection 



4.1 Datasets 

The data presented in the thesis was taken between March 21, 2002, and August 22, 2004 and 
represent 305p6~^ for which the sihcon detector (Section 12.2. 21) [38] , and all three central muon 
systems (CMP, CMU and CMX), described in Section 12.2.41 were operational. 

The /i7 candidates are taken from a logical 'OR' of the inclusive high-px muon sample and the 
inclusive high-Ex photon sample; this was done to ensure a high and stable trigger efficiency for 
the muons. For consistency, 67 candidates are also obtained from a logical 'OR' of the inclusive 
high-Ex electron sample and the inclusive high-Ex photon sample. Each of the samplei^ was 
ntupled using the UC flat ntuple [M] . 

To accept an event from the inclusive high-px lepton sample we require the event to have a loose 
lepton and a photon, or two leptons (either tight or loose), or a tight lepton and ^x > ISGe^ 
(see Tables I5J1IO lO and rO) . 

To accept an event from the inclusive high-px photon sample we require the event to have a 
tight photon (see Table El]) and a loose lepton. The muon selection criteria are listed in Table [5?T| 
the electron selection criteria are listed in Tables 16.11 and 16. 2t 

We check data integrity during the run by plotting the stability of the event yields for the 
control samples. We use the 8 time intervals defined in Table 14.11 [55]. The boundaries of the 
intervals have been chosen to correspond to shutdowns or to major changes in the trigger table. 
The luminosity in each bin is plotted in Figure 14.11 

To summarize, in the resulting inclusive electron and muon samples every event contains either 
a /7 candidate, or a candidate for the W and Z° control samples, described in detail below in 
Chapters [5] and [61 

^We used bhelOd as high-Ex electron sample; bhmuOd as high-px muon sample; and cphlOd as high-Ex photon 
sample 
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Run 


Table 


Date 


£, pb-^ 


Comment 


141544 
152949 


PHYSICS_1_01 4_275 
PHYSICS_1_02 7_175_323 


2002.03.23 
2002.10.16 


15.8 ± 0.9 


Start 
TrigTable 


152593 
156487 


PHYSICS_1_03 1_185_325 
PHYSICS_1_03 2_194_329 


2002.10.16 
2003.01.12 


36.8 ± 2.2 


TrigTable 
Shutdown 


159603 
163113 


PHYSICS_1_04 4_255_357 
PHYSICS_1_04 9_288_373 


2003.02.28 
2003.05.19 


45.9 ± 2.8 


Startup 
TrigTable 


163130 
166325 


PHYSICS_1_05 1_290_375 
PHYSICS_1_05 3_298_382 


2003.05.19 
2003.07.18 


30.1 ± 1.8 


TrigTable 
TrigTable 


166328 
168889 


PHYSICS_1_05 5_319_391 
PHYSICS_1_05 8_345_402 


2003.07.18 
2003.09.06 


39.0 ± 2.3 


TrigTable 
Shutdown 


175066 
179056 


PHYSICS_1_05 8_345_404 
PHYSICS_2_01 4_416_424 


2003.11.26 
2004.02.13 


42.4 ± 2.5 


Startup 
COT bad 


182843 
184835 


PHYSICS_2_05 1_475_455 
PHYSICS_2_05 11_508_473 


2004.05.19 
2004.07.06 


41.1 ± 2.5 


COT good 
TrigTable 


184868 
186598 


PHYSICS_2_05 11_508_473 
PHYSICS_2_05 17_531_484 


2004.07.07 
2004.08.22 


50.0 ± 3.0 


TrigTable 
Shutdown 



Table 4.1: The intervals in run number used to check the stability of W and Z event yields versus 
time. The boundaries of the intervals have been chosen to correspond to shutdowns or to major 
changes in the trigger table. The offline luminosity is shown; the luminosity scale factor of 1.019 
has not been applied. The systematic error of 6% (Chapter [TTl) in the C is shown. 

4.2 Selection Overview 

Events with a high transverse momentum (px) lepton or photon are selected by a three-level 
trigger [13] that requires an event to have either a lepton with px > 18 GeV or a photon with 
Ex > 25 GeV within the central region, \ri\ < 1.0. The trigger system selects photon and electron 
candidates from clusters of energy in the central electromagnetic calorimeter. Electrons are further 
distinguished from photons by requiring the presence of a COT track pointing at the cluster. The 
muon trigger requires a COT track that extrapolates to a reconstructed track segment ("stub") 
in the muon drift chambers. 

We use the same kinematic event selection as in the Run I analysis: inclusive £7 events are 
selected by requiring a central photon candidate with E^ > 25 GeV, a central lepton candidate 
(e or /i) with E^ > 25 GeV passing the "tight" criteria listed below, and a point of origin along 
the beam-line not more than 60 cm from the center of the detector. 

A muon candidate (Chapter [5]) passing the "tight" cuts has the following properties: a) a 
well-measured track in the COT; b) energies deposited in the electromagnetic and hadron com- 
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Figure 4.1: The luminosity in each of the 8 time intervals defined in Table 14.11 used to check 
the stability of control event yields during the run. The boundaries of the intervals have been 
chosen to correspond to shutdowns or to major changes in the trigger table. The total luminosity 
is 305 pb~^. The 6% luminosity uncertainty (Chapter [TTj) is not shown. 

partments of the calorimeter consistent with expectations; c) a muon "stub" track in the CMX 
detector or in both the CMU and CMP detectors [13] consistent with the extrapolated position 
of the COT track; and d) COT timing measurements consistent with a track from a pp collision 
and not from a cosmic ray. 

An electron candidate (Chapter [6]) passing the "tight" selection has the following properties: a) 
a high-quality track with px of at least half the shower energy, unless the Et > 100 GeV, in which 
case the pt threshold is set to 25 GeV; b) a transverse shower profile consistent with an electron 
shower shape and that matches the extrapolated track position; c) a lateral sharing of energy in 
the two calorimeter towers containing the electron shower consistent with that expected; and d) 
minimal leakage into the hadron calorimeter. 

Photon candidates (Chapter [7j) are required to have no track with px > 1 GeV, and at most 
one track with px < 1 GeV, pointing at the calorimeter cluster; good profiles in both transverse 
dimensions at shower maximum; and minimal leakage into the hadron calorimeter. 

To reduce background from photons or leptons from the decays of hadrons produced in jets, 
both the photon and the lepton in each event are required to be "isolated" . The Et deposited in 
the calorimeter towers in a cone in rj — (p space of radius R = OA around the photon or lepton 
position is summed, and the Ex due to the photon or lepton is subtracted. The remaining Et in 
the cone is required to be less than 2.0 GeV + 0.02 x (Ex — 20 GeV) for a photon, or less than 
10% of the Et for electrons or pt for muons. In addition, for photons the sum of the Pt of all 
COT tracks in the cone must be less than 2.0 GeV + 0.005 x Et. 
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Missing transverse energy (Section 18.11) is calculated from the calorimeter tower energies in 
the region \ri\ < 3.6. Corrections are then made to the for non-uniform calorimeter response [56] 
for jets with uncorrected Et > 15 GeV and rj < 2.0, and for muons with px > 20 GeV. 

We use and Z° production as control samples (see Section [572] for the details for the muon 
channel and Section [672] for the electron channel) to ensure that the efficiencies for high-px electrons 
and muons, as well as for ^rp, are well understood. The photon control sample is constructed from 
the events in which one of the electrons radiates a high-Ex photon, with an additional requirement 
that the 67 invariant mass be within 10 GeV of the Z° mass. 

The first search we perform is in the ^7^t + subsample, defined by requiring that an event 
contain > 25 GeV (Section 18. ip in addition to the photon and "tight" lepton. 

A second search, for the i£^+X signature, is constructed by requiring another muon (Chapter [5]) 
or electron (Chapter [6]) in addition to the "tight" lepton and the photon. The additional muons 
are required to have px > 20 GeV and to satisfy at least one of two different sets of criteria: 
the same as those above for "tight" muons but with fewer hits required on the track, or a more 
stringent cut on track quality but no requirement that there be a matching "stub" in the muon 
systems. Additional central electrons are required to have Ex > 20 GeV and to satisfy the same 
criteria as tight central electrons but with a track requirement of only px > 10 GeV (rather 
than 0.5 X Ex), and no requirement on a shower maximum measurement or lateral energy sharing 
between calorimeter towers. Electrons in the end-plug calorimeters (Section l6.1.3p . 1.2 < \ri\ < 2.0, 
are required to have Ex > 15 GeV, minimal leakage into the hadron calorimeter, a "track" 
containing at least 3 hits in the silicon tracking system, and a shower transverse shape consistent 
with that expected, with a centroid close to the extrapolated position of the track [57] . 

The analysis includes a search for efij events, for which the estimated SM expectation is of 
order of 0.2 events. We also search for ^77 events by requiring another photon with Ex >25 GeV 
in addition to the "tight" lepton and the photon. The additional photons are required to pass 
standard photon cuts, described in Chapter [71 
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Chapter 5 



Muon Identification and Control 
Samples 

This chapter describes the selection of muon objects that are used both in the searches and for 
the control samples. We require at least one 'tight central muon' in an event for it to be classified 
as a /i7 event. In both 67 and /i7 events we search for additional muons using a definition of 
'loose central muon'. In this chapter we describe these two sets of cuts and the numbers of muon 
objects passing each cut below. As this is a chapter on object identification, the tables show the 
number of objects passing each cut. 

The summary on the number of events in the Muon Sample is shown in Table [531 The counting 
experiments based on event counts rather than object counts for the /i7 candidates are described 
in Chapter [121 This chapter also describes the control samples of — > fi^u and Z° f^~^fJ'~ 
decays used to check the temporal stability, and also the product of acceptance and efficiency. 

5.1 Muon Selection Criteria 

The muon selection cuts are similar to the standard CDF Run II cuts [58]. We describe the selec- 
tion of the tight muons in Section [5. 1.11 The selection of loose muons is described in Sections 15. 1.21 
and I5.1.3[ 

5.1.1 Tight Central CMUP and CMX Muons 

The muon selection criteria for a tight central muon are listed in Tables ISTTl and are described below. 
Tight central muons are identified by extrapolating tracks in the COT through the calorimeters, 
and the extrapolation is required to match to a stub either in both the CMU and CMP muon 
detectors ('CMUP' muon) or in the CMX system ('CMX' muon), see Table [5^ . Tight central 
muons are required to have a track-stub matching distance less than 3 cm for CMU, less than 5 
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Variable 


Tight 


Loose 


Stubless 


Track Pt 


> 25 Gey 


> 20GeV 


> 20GeV 


Track quality cuts 


3x3SLx5 hits 


3x2SLx5 hits 


3x3SLx5 hits 


Track \zO\ 


< 60 cm 


< 60 cm 


< 60 cm 


Calorimeter Energy (Em) 


< 2 + sliding 


< 2 + sliding 


< 2 + sliding 


Calorimeter Energy (Had) 


< 6 + sliding 


< 6 + sliding 


< 6 + sliding 


Fractional Calorimeter Iso- 
lation Et 


< 0.1 


< 0.1 


< 0.1 


Cosmic 


False 


False 


False 


Clii2/(N of COT hits-5) 


- 


- 


< 3 


Cal.Energy (EM+Had) 






> 0.1 


CMUP muons cuts 


yes 


yes 


no 


CMX muons cuts 


yes 


yes 


no 



Table 5.1: Muon identification and isolation cuts for Tight, Loose, and Stubless muons. Tight 
and Loose are further subdivided into CMUP and CMX categories. CMUP muon cuts are: 
\AX{CMU)\ < 3 cm, \AX{CMP)\ < 5 cm. CMX muon cuts are: \AX{CMX)\ < 6 cm, 
COT exit radius of the muon track p(COT) > 140 cm. 

cm for CMP, and less than 6 cm for CMX. For a CMX muon we also require COT exit radius of 
the muon track p(COT) to be greater than 140 cm to ensure that the track is well-measured. 

The CMUP and CMX muon identifications require a muon object with the requisite muon 
stubs. There are 355105 such objects in the 370679 events of the muon sample. These are then 
divided into CMUP muon candidates and CMX muon candidates. Stubless muon candidates are 
treated separately in Section [5. 1.3t there are 55346 stubless muon objects in the 370679 events. 

The impact parameter calculation uses the default muon track rather than the parent COT 
track, and a tighter impact parameter cut is applied if the track does in fact contain silicon 
hits. Instead, we have tabulated this do cut for reference but we do not use it to select tight or 
loose muons. The muon tracks used in the initial selection for this analysis are beam-constrained 
COT-only [58j. For default muon tracks that contain silicon we link backwards to the COT-only 
parent track and use that track for all subsequent analysis. This technique, while losing valuable 
information from the silicon at this stage, puts all prompt COT tracks on the same footing. 

The Tj — (f) distributions are shown in Figure 15.11 for the muons that pass the loose muon 
identification cuts (Table 15.11) . Muon candidates which have stubs reconstructed from hits in 
either the 'bluebeam', 'miniskirt' or 'keystone' regions [58] of the detector are rejected. These 
sections of the detector were not fully operational for the entire data sample. 

All central muons are required to have l^ol < 60 cm so that the collision is well-contained within 
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Variable 


Cut 


Cumulative 


This Cut 


II Candidates 




355105 


355105 


Track Pt 


> 25 GeV 


309679 


309679 


Track quality cuts(loose) 


3x2SLx5 hits 


292208 


314742 


Track quality cuts(tight) 


3x3SLx5 hits 


284759 


304392 


Track zq 


< 60 cm 


280969 


348881 


Calorimeter Energy (Em) 


< 2+max (0, 0.0115X (p-100)) GeV 


278614 


341205 


Calorimeter Energy (Had) 


< 6+max (0, 0.028 x (p-100)) GeV 


277925 


346820 


Fractional Calorimeter 


< 0.1 


277508 


321987 


Isolation Et 








Cosmic 


FALSE 


194038 


248861 


CMU stub 


TRUE 


113369 


201714 


|AX CMU\ 


< 3 cm 


113057 


347140 


CMP stub 


TRUE 


111498 


203653 


|AX CMP\ 


< 5 cm 


111407 


347649 


Region is OK 


TRUE 


111407 


349193 


Track do 


< 0.02 cm (si hits) OR 

< 0.2 cm (no si hits) 


101933 


211425 


CMX stub 


TRUE 


76066 


89334 


|AX CMX\ 


< 6 cm 


75989 


351303 


Region is OK 


TRUE 


73367 


349193 


Track o?o 


< 0.02 cm (si hits) OR 

< 0.2 cm (no si hits) 


60554 


211425 



Table 5.2: Tight CMUP and CMX muon identification. The cumulative totals of tight central 
CMUP and CMX muons, showing the behavior as the cuts are applied. The initial entries in the 
table, before the rows in which the stub requirement is applied, start with the total number of 
muon objects with muon stubs in either the CMUP or CMX systems. Each entry corresponds to 
a muon in the CDF Muon Collection in an event. The "Region is OK" cut for the CMUP muons 
requires \AX{CMU)\ < 3 cm, \AX{CMP)\ < 5 cm. The "Region is OK" cut for the CMX 
muons requires |AX(CMX)| < 6 cm, COT exit radius p(COT) > 140 cm. The (Iq cut is not 
applied to select muons, but is tabulated for reference and to see its effect. The column labeled 
'Cumulative' gives the effect of each successive cut on the number of muon candidates in the 
370679 events in the muon subsample; The heading 'This Cut' represents the effect of applying 
only the cut listed. 
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DetEta Vs Phi for mu 
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(a) 1] — (f) distribution for the muons which pass 
loose muon identification cuts 



Figure 5.1: (a) The rj — 
in Table O 



distribution for the muons which pass muon identification cuts, hsted 



the CDF detector. In order to be well-measured, the muon track is required to have minimum of 
3 axial and 3 stereo superlayers with at least 5 hits in each superlayer. 

High energy muons are typically isolated 'minimum-ionizing' particles that have limited calorime- 
ter energy. A muon traversing the central electromagnetic calorimeter(CEM) deposits an average 
energy of ~ 0.3 GeV. Therefore we require muon candidates to deposit less than 2 GeV total in 
the CEM towers (we take into account two towers in the CEM) the muon track intersects. Simi- 
larly, muons transversing the central hadronic calorimeter (CH A) deposit an average energy of ~ 2 
GeV; we consequently require muon candidates to deposit a total energy less than ~ 6 GeV, also 
increasing with muon momentum, in the CHA towers intersected by the track extrapolation. To 
take into account the (slow) growth of energy loss with momentum, for very high energy muons 
{p > 100 GeV) we require the measured CEM energy to be less than 2.0 + 0.0115 x (p — 100) GeV 
and CHA energy to be less than 6.0 + 0.028 * {p - 100) GeV. 

To suppress hadrons and decay muons created from hadrons in jets we require the total trans- 
verse energy deposited in the calorimeters in a cone of R=0.4 around the muon track direc- 
tion(known as the fractional calorimeter isolation Et) to be less than 0.1 of the muon track pt. 

The COT cosmic finder by itself is essentially fully efficient. Therefore, to suppress cosmic 
rays we use the COT-based cosmic rejection [59j and reject events which it tagged as cosmic ray 
muons. 
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Variable 


Cut 


Cumulative 


This Cut 


fi Candidates 




74888 


74888 


Track Pt 


> 20 GeV 


42701 


42701 


Track quality cuts(loose) 


3x2SLx5 hits 


25328 


39120 


Track quality cuts(tight) 


3x3SLx5 hits 


22362 


34366 


Track \zo\ 


< 60 cm 


21041 


72372 


Calorimeter Energy (Em) 


< 2+max (0, O.OllSx (p-100)) GeV 


18634 


62363 


Calorimeter Energy (Had) 


< 6+max (0, 0.028 x (p-100)) GeV 


18337 


67502 


Fractional Calorimeter 


< 0.1 


18115 


46783 


Isolation Et 








Cosmicu 


FALSE 


6427 


55172 


Calorimeter Energy 


> 0.1 GeV 


6175 


56551 


(Em+Had) 








Track quality cuts 


xV(N of COT hits-5) < 3 


6052 


33472 


Region is OK 


TRUE 


6052 


74888 


Track \do\ 


< 0.02 cm (si hits) OR 

< 0.2 cm (no si hits) 


6024 


35160 



Table 5.3: Stubless muon identification and isolation cuts. Each entry corresponds to a 'Stubless' 
muon in the CDF Muon Collection. The do cut is not apphed to select muons, but is tabulated for 
reference and to see its effect. The column labeled 'Cumulative' gives the effect of each successive 
cut on the number of stubless muon candidates in the 370679 events in the muon subsample. The 
heading 'This Cut' represents the effect of applying only the cut listed. 

5.1.2 Loose Central CMUP and CMX Muons 

While each fi'j event has to contain at least one tight CMUP or CMX muon, both ej and fi'j events 
are searched for additional high-pT muons that could come from the decays of heavy particles. 
There are two types of secondary muons we accept: 'Loose' CMUP and CMX muons, described 
here, and stubless muons (see Section [5.1.3p . 

Loose muons are muon objects with either CMUP or CMX stubs, but with looser COT cuts 
than the tight CMUP or CMX muons (see Table [5TT]) . We require 3 axial and 2 stereo COT super 
layers with at least 5 hits each for loose CMUP and CMX muons. 

5.1.3 Loose Central Stubless Muons 

The cuts for the Stubless muons, described in Table 15. 3[ are looser than the tight cuts, and in 
particular do not require a stub in the muon chambers. There are three types of 'Stubless' muons: 
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• CMU muons (muon track matches the CMU stub only) 

• CMP muons (muon track matches a stub in the CMP only) 

• CMIO muons (muon track does not match a stub in CMU, CMP or CMX) 

To identify stubless muons, we require at least some energy in the calorimeter towers that the 
muon extrapolates to. Calorimeter Energy (Em+Had) > 0.1 GeV, and a good fit to the COT track, 
•y^ I (Number of COT hits-5)<3 [58]. These two cuts are used to reject charged kaon decays in flight 
in which a low-momentum kaon (~ 5 GeV, typically) decays inside the COT with the kaon and 
decay-muon tracks forming a 'seagull' pattern which is reconstructed as a single high- momentum 
track. 

The pattern-finding algorithm often removes a complete stereo layer in order to get a good 
fit, and so these tracks are badly mis-reconstructed in polar angle. They consequently often are 
recorded to leave zero energy in the extrapolated traversed calorimeter towers [60] . 

5.2 Muon Control Samples 

The W and provide control samples for the ^7 samples. We use data triggered by the high-px 
muon trigger (MU0N_CMUP_18 or MU0N_CMX_18) for both Z° and W candidates, where tight 
muons are the muons that pass the cuts in Table 15.21 

For comparisons with data we used the Z° /i+yU^ Monte Carlo sample [6T|. We applied the 
trigger efficiencies and scale factors [58] , listed in Table 15.41 



Trigger efficiencies 


trigger 
^CMUP 


(90.78 ± 0.47)% 


trigger 
^CMX 


(96.49 ± 0.40)% 


Scale factors 


tight CMUP 


0.8738 ± 0.0086 


tight CMX 


0.9889 ± 0.0063 


loose CMUP 


0.8921 ± 0.0088 


loose CMX 


0.9990 ± 0.0060 


Stubless 


0.9760 ± 0.0026 



Table 5.4: Scale factors and trigger efficiencies for the muons, applied to TP ji'^ ji MC sample. 
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5.2.1 The Z" ji^ jT Control Sample 

The selection criteria for the TP —>■ fj^^fj^^ control sample are listed in Table 15.51 We require 
two muons in the event. One must pass the tight cuts (Table [52]), and another must pass the 
loose cuts (Table \5A\i . We further require the two muons to have opposite charge, and require 
the difference of the Zq beam-line coordinates of the muon tracks to be less than 4 cm. The last 
requirement for counting events is that the invariant mass of the muon pair of the candidate 
should be between 66 GeV and 116 GeV. We find 9175 Z^ events. We also find 5 same-sign 
events, indicative of the maximum level of track reconstruction problems. The selection variable 
distributions for the Z° fx^fi^ control sample are shown in Figure 15.21 



Variable 


Cut 


Tight Muon 


Cuts are listed in Table |5.2| 


Loose Muon 


Cuts are listed in Tables 15.21 with looser COT cuts. Table 15.31 


Muon tracks must be 
of opposite charge 




Delta Z cut 


1 7^^ 7^^ ^ A nm 
l^TRACK ' ^TRACK ^ ^ 


Mass Window cut 


M > 66 GeV and < 116 GeV 


Trigger 


MU0N_CMUP_18 or MU0N_CMX_18 



Table 5.5: The selection cuts for the Z° — > /i+yU control sample. The superscripts /il and /i2 
stand for the 2 muons in the event. 

To normalize the Z° /^^/^^ MC sample, we used the measured a{7P —>■ /i+yU^)xBR [62j. The 
comparison of data vs MC for different muon types is shown in Table 15.71 The Z° — ^ /^"'"A*" event 
yields are shown in Figure 15.31 

5.2.2 The W± ^ fi^iy Control Sample 

The selection criteria for the W^ —>■ fi^u control sample are listed in Table 15.61 We require one 
tight muon (px > 25 GeV, see Table and > 25 GeV (Section EH]). We require the 
transverse mass of the W candidate to be in the mass window 20-140 GeV. 

To reject Z° fi~^fi~ events in which one muon is not identified, events with a second track 
with pt > 10 GeV and associated EM and HAD calorimeter energies less than 3 and 9 GeV, 
respectively, are rejected. We find 118321 W"*^ — > /x^z/ events, 59387 positive and 58934 negative 
VT's. The W^ fi^u event yields are shown in Figure E31 
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I Z Invariant Mass[GeV] ] 



Entries 


9175 
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Entries 
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Figure 5.2: The Z° — > control sample 'sanity-check' plots. The distributions for Invariant 

Mass of Z° Pt of Z° linear plots(a, c), log plots(b, d). The histogram is the 

prediction from the Z° — yU+yU^ MC sample; the points are the Z° fi^fi~ candidates from the 
data. Background estimates are not included. 



Variable 


Cut 


Tight Muon 


Cuts are listed in Table |5.2| 


Missing Cut 


Et > 25 GeV 


Z-Rejection 


No 2nd Track with Pt > 10 GeV, Eem < 3 GeV, Ehad < 9 GeV 


Mass Window Cut 


Mt > 20 GeV and < 140 GeV 


Trigger 


MU0N_CMUP_18 or MU0N_CMX_18 



Table 5.6: The selection cuts for the n^u control sample. 
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5.2.3 Summary of the Muon Control Sample Event Counts. Stability 
Plots 



We use the control samples of jji^u and Z° decays to check temporal stability of 

the event yields (Figure [5l3|) . The summary on the number of events in the Muon Sample is shown 
in Table ESI 



Entries vs Runs for stability 




Entries vs Runs for stability I Entries 

' — H Overflow 

/ ndf 
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Figure 5.3: Stability plots of: (a) the Z° —>■ fJ^^fi", and (b) —>■ fi^u control sample cross 
sections versus run number. The bins are those of Table 14.11 Shown are the statistical errors, 
the luminosity systematic error of 6% (Chapter [TT!) is not included. The trends in Z*^ — > fi~^fi~ 



(W^ — ^ fi^jy) are similar to e^e 
effects (luminosity, trigger, COT). 



(W^ 



see Figure 16.41 We attribute this to common 





DATA 


MC 








ZO 


— i> 11^ jjT 




ZO- 






ALL Tight 


3465 


3403 


ZO- 






CMUP-CMUP Tight 


1462 


1428 


ZO- 






CMUP-CMX Tight 


1533 


1529 


ZO- 




V- 


CMX-CMX Tight 


470 


446 


ZO- 




V- 


ALL 


9175 


8784 


ZO- 






CMUP-STUBLESS 


3594 


3464 


ZO- 






CMX-STUBLESS 


1929 


1868 



Table 5.7: Z'^ fi^jj, summary: data vs Z° ji'^ ji MC. Background estimates are not included. 
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Category 


N 


Events with at least one Tight Muon 


180340 


Tight Muons 


183982 


— > /i^i/(triggered+Mvi/ window cut) 


118321 


Tight + Loose Muons 


9777 


Z° — > (triggered+M^ window cut) 


9175 



Table 5.8: The numbers of events for the muon control samples. The muon selection cuts are given 
in Table ISTTj the Z° fi~^fi~ selection cuts in Table 1531 and the fi^u cuts in Table [531 
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Chapter 



6 



Electron 



Identification and Control 



Samples 



We require at least one 'tight central electron' in an event for it to be classified as an 67 event. 
In both 67 and /i7 events we search for additional 'loose' electrons in the central and end-plug 
electromagnetic calorimeters. We describe the tight central and loose central and plug cuts below. 
The counting experiments based on event counts rather than object counts for the 67 candidates 
are described in Chapter [121 This chapter also describes the control samples of W^ — e^z/ and 
Z° — > e^e~ decays used to check the temporal stability, and also the product of acceptance and 
efficiency. 

6.1 Electron Selection Criteria 

The electron selection cuts are similar to the standard CDF Run II cuts [63]. We describe the 
selection of the tight central electrons in Section 16.1.11 The selection of loose central electrons 
is described in Section I6.1.2t the selection of electrons in end-plug is presented in Section 16.1.31 
The selection cuts are standard [63] with the exception that the CES fiducial requirement (see 
Section 16.1.11) and the conversion cut (see Sections IA.6I and 110. 2p are not applied. 

6.1.1 Tight Central Electrons 

The selection criteria for tight central electrons are listed in Table 16.11 and are described below. 

Electrons are identified in the CEM by matching high momentum tracks to high-energy CEM 
clusters. The electron track is the highest momentum track which intersects one of two towers in 
the CEM cluster. The electron tracks that we use in this analysis are beam-constrained COT-only. 
We apply the same corrections to the electron tracks as we do to the muon tracks. 

Fiduciality is a variable, which can have following values: 
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-1 : error (null strip/wire cluster) 

: not fiducial in central or plug 

1 : fiducial in central or plug 

2 : fiducial, but in CEM Tower 9 

3 : fiducial, but in Chimney wedge tower 7 

4 : fiducial in CEM using max pt track extrapolated to plane of CES 

5 : fiducial in PEM using PES 

6 : fiducial in PEM using max px track extrapolated to plane of PES 

An electron candidate is required to have tracking momentum (P) which exceeds half of its 
calorimeter energy (E). The electron track is required to have a minimum of 3 axial and 2 stereo 
SL segments containing at least 5 hits each. In order that the momentum resolution docs not 
make for inefficiencies for very high-energy electrons, for Et > 100 GeV the E/P cut is not 
applied (leaving only the px > 25 GeV cut as the requirement on the track). The electrons are 
required to have the track extrapolate to the beam line within l^ol < 60 cm so that CDF detector 
contains the collision well. 

The position of the track extrapolated to the CES radius must satisfy the following require- 
ments: it must fall within charge-signed CES shower position of the cluster in the r-phi view 
-3.0 cm < Qtrk X AX < 1.5 cm and it must fall within 3 cm of the CES shower position in the 
Z-direction( AZ) . 

The CEM shower characteristics should be consistent with that of a single charged particle. 
We require the ratio of the total energy of the CHA towers located behind the CEM towers in 



Variable 


Tight 


TightlOO 


Loose 


Et, GeV 


> 25 


> 100 


> 20 


Track Pt, GeV 


> 10 


> 25 


> 10 


Track \zO , cm 


< 60 


< 60 


< 60 


Had/Em 


<0.055+0.00045xE 


<0.055+0.00045xE 


<0.055+0.00045xE 


E/P 


< 2.0 






Lshr 


< 0.2 






Chi2 Strips 


< 10 






AX, cm 


-3.0 < Qtrk^AX <1.5 


AX <3.0 




AZ , cm 


< 3.0 


< 5.0 




Fractional Calorime- 
ter Isolation Et 


< 0.1 


< 0.1 


< 0.1 


Track quahty cuts 


3x2 SL x5 hits 


3x2 SL x5 hits 


3x2 SL x5 hits 



Table 6.1: Central electron identification and isolation cuts. 
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Variable 


Tight 


Phoenix Tight 


Et, GeV 


> 15 


> 15 


Had/Em 


< 0.05 


< 0.05 


Fractional Calorimeter Isolation Et 


< 0.1 


< 0.1 


Chi2 Strips 


< 10 


< 10 


AR^y, cm 


< 3.0 


< 3.0 


PES 5by9 U and V 


> 0.65 


> 0.65 


PEM \r]\ 


2.0 < \r]\ < 1.2 


2.0 < \r]\ < 1.2 


PhxMatch 




TRUE 


Number of Silicon Hits 




> 3 


Z(Phoenix) , cm 




< 60 



Table 6.2: Identification and isolation cuts for the electrons in end-plug. We are using the "Phoenix 
Tight" selection [6] 



the electron cluster to that of the electron itself to be less than than 0.055+0.00045 xE GeV. A 
comparison of the lateral shower sharing with neighboring towers in the CEM cluster with test- 
beam data is parameterized by a dimensionless quantity, Lshr [39], which must have a value less 
than 0.2. 

We require the for the profile of energy deposited in the CES strips compared to that 
expected from test beam data [32] to be less than 10. No cut is made on the profile in the CES 
wires as bremsstrahlung will separate from the electron in the r0 view. 

As an additional isolation requirement, the total transverse energy deposited in the calorimeter 
in a cone R=0.4 around the electron track, must be less than 0.1 of the E^ of the electron. The 
isolation is corrected via the standard algorithm [65] , for leakage, but not the number of vertices. 

We do not apply 'Conversion Flag' and 'Fiducial' cuts to select electrons, they are tabulated 
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Figure 6.1: Plots for electrons which pass cuts all but CES fiducial requirement. 
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Variable 


Cut 


Cumulative 


This Cut 


e Candidates 




470379 


470379 


Et, GeV 


> 25 


323506 


323506 


Track Pt, GeV 


> 10; 

Et > 100 GeV: > 25 


316943 


324739 


Track \zo\, cm 


< 60 


316930 


468731 


Had/Em 


< 0.055 + 0.00045xE 


315752 


420940 


E/P (for < 100 GeV only) 


< 2.0 


308932 


365229 


Lshr 


< 0.2 (for Et < 100 GeV only) 


307955 


362783 


Clii2 Strips 


< 10 (for Et < 100 GeV only) 


306883 


433201 


AX, cm 


-3.0 < Qtrk X AX < 1.5; 
Et > 100 GeV: |AX| < 3.0 


306768 


362872 


AZ , cm 


< 3.0; 

Et > 100 GeV: < 5.0 


306660 


361208 


Fractional Calorimeter 
Isolation Ej- 


< 0.1 


269624 


289848 


Track quality cuts 


3 Axial, 2 Stereo SL x 5 hits 


269621 


398424 


Conversion Flag 




252900 


439187 


Fiducial 


based on shower max 


222844 


367024 



Table 6.3: Tight central electron identification and isolation cuts. Each entry corresponds to an 
CDF Central EM object in an event. The column labeled 'Cumulative' gives the effect of each 
successive cut on the number of electron candidates in the 490355 events in the electron subsample. 
The heading 'This Cut' represents the effect of applying only the cut listed. We do not apply 
'Conversion Flag' and 'Fiducial' cuts to select electrons, they are tabulated for reference and to 
see their effects. 

for reference and to see their effects. 

The acceptance gain by removing the fiduciality requirement is approximately 14%. Figure EH] 
shows the distributions for the electrons which pass all cuts but fiducial requirement. 

6.1.2 Loose Central Electrons 

While each 67 event has to contain at least one tight electron, both 67 and fij events are searched 
for additional high-px electrons that could come from the decays of heavy particles. The cuts for 
these additional electrons are described in Table 16. 1[ These cuts are looser than the tight cuts, 
and in particular do not require any of the CES variables, i.e. no track-cluster match in AX or 
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AZ and no cut on strip x^, and also no cut on LsHt [39]. 
6.1.3 Plug Electrons 

Additional isolated electrons in the plug calorimeter with Et > 15 GeV are identified for measured 
PEM rapidities of 1.2 < \ri\ < 2.0. The cuts used for plug electron identification are given in 
Table [631 We require minimal leakage or activity in the hadron calorimeter, Had/Em < 0.05, a 
fractional isolation (isolation energy over the electron energy) less than 0.1, and the shower shape 
to satisfy the PEM 3x3 and PES 5by9 5-strip to 9 strip ratio cuts. These cuts are similar to 
standard cuts [63] . 

We apply face corrections to the PEM energy of the plug electron candidate, add the PPR 
energy and scale resulting number by 1.0315 [66j, as shown in Equation 16.11 

Eplug electron {ElZ^ + E,,r) X 1.0315 (6.1) 



6.2 Electron Control Samples 

As in the muon case, the W and provide control samples for the 67 samples. We require at 
least one tight electron pass the high-Ex electron trigger (ELECTR0N_CENTRAL_18) for both 
W and candidates, where tight electrons are the electrons that pass the cuts in Table 16.11 



Variable 


Cut 


Cumulative 


This Cut 


e Candidates 




54247 


54247 


Ej' 


> 15 GeV 


15341 


15341 


Had/Em 


< 0.05 


11971 


37048 


Fractional Calorimeter 
Isolation E^ 


< 0.1 


10274 


21053 


Chi2 Strips 


< 10 


9931 


21708 


Delta R 


< 3.0 cm 


9886 


47213 


PES 5by9 U and V 


> 0.65 


9478 


25977 


PEM \r]\ 


1.2 < 7] < 2.0 


9057 


33376 


PhxMatch 


TRUE 


9028 


13654 


Number of Silicon Hits 


> 3 


8836 


11647 


Z (Phoenix) 


< 60 cm 


8836 


11733 



Table 6.4: Identification and isolation cuts for additional plug electrons. 
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Trigger efficiencies 


trigger 
^CEM 


l-2.784*exp[-1.749*(ET - 17.86)] 


trigger 
^Track 


^0.96-0.98 


Scale factors 


tight CEM 


0.999 ± 0.006 


loose CEM 


1.005 ± 0.005 


tight PHX 


0.95 ± 0.01 



Table 6.5: Scale factors and trigger efficiencies for the electrons, applied to the TP — > e^e MC 
samples. 

For comparisons with data we used the Z° e+e~ MC sample [61]. We applied trigger effi- 
ciencies [67j and scale factors [63], which are listed in table [631 The CEM Trigger efficiency for 
high-^^T electrons is l-2.784*exp(-1.749*(ET - 17.86)) [67]. The total Track Trigger efficiency is 
~ 96-98%, depending on the run number and silicon/non-silicon list of good runs [67j.We apply 
this trigger efficiency to our run-dependent Z° — > e^e~ MC sample. 

6.2.1 The e^e~ Central- Central Control Sample 

The selection criteria for the e~^e~ central-central control sample are listed in Table [6l6l For 

this sample a Z° candidate is required to have two central electrons, one passing the tight cuts, 
and the other the loose cuts (Table [O]) . The mass of the Z candidate is required to be within 
the window 66 GeV to 116 GeV. The difference in the zq coordinates of the two electron tracks 
must be less than 4 cm. We find 10128 opposite-sign events and 199 same-sign central-central 
events satisfying these criteria (the large number of same-sign events in the electron sample but 
not in the muon sample is largely due to photon conversions - see Chapter [TOj) . 

The distributions in mass, px of the pair, and A(f) are shown in Figure 16. 2[ To normalize the 
Z° — > e~^e~ MC sample, we used the measured cr(Z° — > yU+yU~)xBR [62] . 

6.2.2 The e+e" Central-Plug Control Sample 

We also form a central-plug dielectron Z^ control sample to monitor the identification performance 
for the electrons in end-plug calorimeters. We require a tight central electron, and a plug electron 
passing the cuts of Table 16. 4[ The mass of the Z candidate is required to be within the window 
66 GeV to 116 GeV. We find 3996 (4004) central-plug Z° events satisfying these criteria with the 
plug electron in the East (West) calorimeter. 

The distributions in mass, px of the pair, and A0 are shown in Figure 16.31 The comparison of 
data vs MC is shown in Table 16. 8[ The Z° — > event yields are shown in Figure 16.41 
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Variable 


Cut 


Tight Electron 


Cuts are listed in Table |6.3| 


Loose Electron 


Cuts are listed in Table |6.1| 


Electron tracks must 
be of opposite charge 


Q-i + Q^^ = 


Delta Z cut 


l^TRACK ~ ^TRACkI < ^ Cm 


Mass Window cut 


M > 66 GeV and < 116 GeV 


Trigger 


ELECTR0N_CENTRAL_18 



Table 6.6: The selection cuts for the CC and CP 
el and e2 stand for electrons in the event. 



e^e control samples. The superscripts 
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Figure 6.2: The distributions for Invariant Mass of Z e e~ , Pt of Z e e~, linear plots(a, 
c), log plots(b, d). The histogram is the prediction from the Z° e^e^ MC sample (ztop2i), the 
points are Z° — >• e~^e~ candidates in the data. Background estimates are not included. 
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Figure 6.3: Central-plug dielectron Z° e^e^ control sample to monitor the plug electron iden- 
tification performance. We require a tight central electronfTable . and a plug electron passing 
the cuts of Table [631 The distributions for Invariant Mass of Z° e+e^, Pt of Z° e~^e~ , linear 
plots(a, c), log plots(b, d). We find 3996 (4004) central-plug events satisfying these criteria with 
the plug electron in the East (West) calorimeter. The points are Central-East Plug Z*- 



e^e 



candidates; the histogram is Central- West Plug 



e ' e 



candidates. 



6.2.3 The W± e^v Control Sample 

The selection criteria for the W^ e^v control sample are listed in Table 16.71 We require one 
tight central electron (i.e. Et > 25 GeV), and corrected (Section 18.11) greater than 25 GeV. 
We require the transverse mass of the W candidate to be in the mass window 20-140 GeV. We find 
184805 W=^ e^u events, 92670 positive and 92135 negative W's. The W^ ^ e^u event yields 
are shown in Figure 16. 4[ 
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Variable 


Cut 


Tight Electron 


Cuts are listed in Table |6.3| 


Missing Et Cut 


Et > 25 GeV 


Mass Window Cut 


Mr > 20 GeV and < 140 GeV 


Trigger 


ELECTR0N_CENTRAL_18 



Table 6.7: ^ e^u selection cuts. 

6.2.4 Summary of the Electron Control Sample Event Counts and 
Stability Plots 

We use the control samples of W^ e^u and Z° e^e^ decays to check temporal stability of 
the event yields (Figure [63!) • Table [63] presents summary of the electron control sample event 
counts. 





DATA 


MC 










z 


^ e+e- 




ZO - 


^e^ 


"e" 


~cc 




10128 


9955 


ZO- 




"e" 


CC Same— Sign 




199 


127 


ZO- 


^e^ 


"e" 


~CE 




3996 


4257 


ZO- 


-^e^ 


"e" 


~CW 




4004 


4131 



Table 6.8: Z° e"^e summary: data vs Z° e'^e MC. The material is underestimated in MC, 
so we estimate e — 7 fake rate from data. 
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Entries vs Runs for stability 




Entries vs Runs for stability I Entries 
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Figure 6.4: Stability plots for Z° ^ e+e-(a), W± e±z/(b), Z° ^ e+e" Central-East (c), 
Z° e+e~ Central-West(d). The bins are those of Table 14. 1[ Shown are the statistical er- 
rors, the luminosity systematic error of 6% (Chapter [TTj) is not included. The trends in 
Z° — i> e"'"e~(W^ e^z/) are similar to Z° ji^ ji~ (W^ f^"^^), see Figure 15. 3[ We attribute 
this to some common effects (luminosity, trigger, COT). 
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Category 


N 


Events with at least one Tight Electron 


254664 


Tight Electrons 


262555 


e'^i/(triggered+Mv(/) 


184805 


Tight + Loose Central Electron 


11091 


Tight + East Plug Electron 


4181 


Tight + West Plug Electron 


4150 


CC ZO ^ e+e-(triggered+Mz) 


10128 


CE Plug Z° e+e-(triggered+Mz) 


3996 


CW Plug Z° e+e-(triggered+M^) 


4004 



Table 6.9: The numbers of events for the electron control samples. The electron selection cuts are 
given in Table 16.31 the —>■ e~^e~ selection cuts in Table 16.61 and the —>■ e^v selection cuts 
in Table 16. 7[ CC refers to Central- Central events; CE refers to Central-East Plug, CW refers to 
Central- West Plug. 
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Chapter 7 

Photon Identification 



The photon selection criteria are identical for photons in both the muon and electron samples; the 
photon cuts, and the number of events passing in each sample, are enumerated in Table 17.21 and 
are described below. 

A sample of photons from Z*^ — > e~^e~ events used to measure the probability of an electron 
radiating an energetic photon is also introduced. 

7.1 Photon Selection Criteria 

A photon candidate is required to have corrected transverse energy greater than 25 GeV. For 
photons or electrons the CES shower position is determined by the energy-weighted centroid of 
the highest energy clusters of those strips and wires in the CES which correspond to the seed 
tower. The direction of the photon is determined by the line connecting the primary event vertex 



Variable 


Cut 


rpcorr 

Ulrp 


> 25 GeV 


Had/Em 


< 0.125 or < 0.055 + 0.00045 xE'^"^'' 


(Strips+Wires)/2.0 


< 20 


N Tracks 


< 1 


Track Pt 


< 1+0.005 xE^°'-'- GeV 


Cone 0.4 IsoEf?^^ 


< 2.0+0.02 x(£;^°^'^ - 20) GeV 


Cone 0.4 Tracklso 


< 2.0+0.005x^1?^'- GeV 


2nd CES Cluster (Strip and Wire) 


< 2.4+0.01 xEfP'''' GeV 


Fiducial 


CeslXl < 21 cm, 9 cm < Ces \Z\ < 230 cm 



Table 7.1: Photon identification and isolation cuts. 
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Subsample 


Muon 


Electron 


Variable 


Cut 


Cumulative 


This 


Cumulative 


This 


7 Candidates 




75026 


75026 


524626 


524626 


GeV 


> 25 


4567 


4567 


333650 


333650 


HAD/EM 


<0.125|| 

0.055+0.00045xE 


4169 


68819 


330675 


499963 


CES 


< 20 


2978 


39468 


279701 


389644 


N Tracks 


< 1 


2500 


58325 


210366 


371793 


Track Pt, GeV 


< 1 +0.005xEt 


1787 


43186 


4073 


96410 


Cone 0.4 IsoE™^^ GeV 


< 2.0+0.02x(Et -20) 


1647 


19068 


3818 


302862 


Track Isolation, GeV 


< 2.0+0.005xEt 


1610 


36352 


3771 


64776 


E (2nd CES Cluster) 
(Strip and Wire), GeV 


< 2.4+O.OIxEt 


1604 


67048 


3762 


471879 


Fiducial: 

Ces X , Ces Z , cm 


Ces|X| < 21, 

9 < Ces|Z| < 230 


1598 


40879 


3735 


409371 



Table 7.2: Photon identification and isolation cuts. Each entry corresponds to one CDF EM 
Object, Central or Plug. The column labeled 'Cumulative' gives the effect of each successive cut 
on the number of photon candidates. The heading 'This' represents the effect of applying only 
the one cut listed. After a final cut requiring a Tight muon to be in the event with the photon 
we find a total of 66 yU7 candidate events. After a final cut requiring a Tight electron to be in the 
event with the photon we find a total of 508 67 candidate events. 

to the shower position in the CES. To ensure that events are well-measured the shower position 
of the photon is required to fall within the fiducial region of the CES so that the shower is fully 
contained in the active region. 

Photon candidates are required to have characteristics consistent with those of a neutral 
electromagnetically-interacting particle. No COT track with px > 1 GeV may point at the photon 
cluster. One track with px < 1 GeV may point at the cluster. 

The variable 'Isoii^™'^''' is the Run I cone Rr|-^p = 0.4 isolation energy with the Run I correction 
to isolation energy due to phi-crack leakage |65]. The tracking isolation variable 'Track Isolation' 
is the sum of the px of tracks in a in a cone in 77 — 93 space of radius R = 0.4 surrounding the 
photon, measured in GeV. 

Table 17.21 summarizes the selection of photons in the muon and electron subsamples (Sec- 
tion SI]). 

For the muon subsample we require the event to be triggered by either a high-Ex muon trig- 
ger (which is MU0N_CMUP_18 or MU0N_CMX_18) or by a high-Ex photon trigger (PHO- 
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TON_25JSO). After requiring a Tight muon (Table IS.ip to be in the event with the photon 
(Table 17. ip we find a total of 66 /i7 candidate events. 

For the electron subsample we require event to be triggered by either high-Ex electron trigger 
(ELECTR0N_CENTRAL_18) or by high-Ex photon trigger (PHOTON_25_ISO). After requiring 
a tight electron(Table \QA^ to be in the event with the photon(Table mi) we find a total of 508 ej 
candidate events. The disparity between this number and the 66 fij events is due to a number of 
causes, in particular hard photon bremsstrahlung of an electron, as discussed later in Section 17.21 

7.2 Introducing the Photon Control Sample 

Unlike for the electron or muon, finding a pure sample of high-Pt photons is difficult; the 'Compton' 
sample of 7-jet events has QCD fake backgrounds, for example. We describe here a control sample 
of high-Pt photons derived from the 

In looking for photons in the electron sample, one has to take into account that the dominant 
source of fake background for 67 events is —>■ e~^e~ production, wherein one of the electrons 
undergoes hard photon bremsstrahlung in the detector material, or the COT fails to detect one 
of the electron tracks, and that electron subsequently passes all of the photon cuts. There are 
approximately 7890 tight central electron pairs in the CDF data, so a photon fake rate as low 
as 1% will give rise to 158 67 background events, which would be unacceptably high for finding 
sources of new physics comparable to W/Z^ + j production. The Run II detector has significantly 
more material inside the outer tracking chamber than the Run I detector had, and so the number 
of '67' events from Z° e^e" production will be significantly higher. We measure this fake rate 
directly from the data (the material in simulated data is underestimated, see Table [6^ . 

A control sample of Z°-like events, e + '7', is selected from the 397 two-body 67 candidates (see 
Chapter [T2l) by requiring that the invariant mass of the 67 pair, Mg^, be within 10 GeV of the Z^ 
mass (91 GeV). There are 209 such events in the CDF data. It is observed that the shapes of the 
distributions of the two samples (Z° e^e~ and e -|- '7') are similar to each other (Figure 17. ip . 
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Figure 7.1: The distributions for: (a) Invariant Mass of e + '7', (b) pt of e + '7', (c) A0 of e + '7', 
(d) AR of e + '7'. The points are the Z^-like 67 sample; the shaded histogram represents electron- 
electron events from data with the same kinematic cuts, normalized to the number of events in 
the control sample. Invariant mass for e + '7' is slightly shifted to the left, because energy of the 
radiated photon is less than energy of the original electron. 
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Chapter 8 



Calculating the Missing Transverse 
Energy and Hj^ 

This chapter describes how missing transverse energy Ox) is calculated, and gives the definition 
and describes calculation of the total transverse energy. 

8.1 Calculating the f^^^ 

Missing transverse energy (^t) is associated with particles that escape detection. For example, 
is the signature of weakly interacting neutral particles such as neutrinos, or possible new 
particles such as the gravitino or LSP. It also can come from mismeasurement of the true Et of 
objects, or from backgrounds such as cosmic rays or beam halo. 

Missing Et (-^t) is defined by = — Xli E^^t, where i is the calorimeter tower number for 
\7]\ < 3.6, and rij is a unit vector perpendicular to the beam axis and pointing at the i*'^ tower. 
We define the magnitude = \^t\- 

Corrections are made to the for non-uniform calorimeter response [56] for jets with uncor- 
rected Et > 15 GeV and t] < 2.0, and for muons with pt > 20 GeV: 

• Muons: 

- correct for Ex — Pt, where E^ is a transverse energy deposited in electromagnetic and 
hadron calorimeters, and Pt is a transverse momentum of a muon track. We correct ^t fo'^ 
all muons with Et > 20 GeV. 

• Jets: 

- correct for Et — Et^'^, where Et is a transverse energy of an uncorrected jet, and Et ""^ is 
a transverse energy of a jet, corrected for non-uniform calorimeter response. We correct for 
jets with E^'-"- > 15 GeV. 
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When identifying jets we check that jet object does not have any of the objects identified in 
the current anafysis close to it (within AR < 0.5). 

For the £7^t anafysis we set the cut on to be > 25 GeV. 

8.2 Calculating the Total Transverse Energy 

Total transverse energy Ht is defined for each event as the sum of the transverse energies of the 
leptons, photons, jets, and that pass the analysis selection criteria. To calculate we use 
Tight and Loose Central Electrons (Table WAh . Tight Phoenix Electrons (Table [6l2|) . Tight and 
Loose CMUP and CMX muons, Stubless muons (Table ISTTl) . ^rp, and jets in the event with |?7| < 2 
and E^r'"''' > 15. 
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Chapter 9 

Standard Model Predictions 



The dominant source of £7 events at the Tevatron is electroweak diboson production (Figure [9?T]) . 
in which an electroweak boson {W or decays leptonically {iu or i£) and a photon is radiated 
from either an initial state quark, a charged electroweak boson (W), or a final state lepton. The 
number of £7 events from electroweak diboson production is estimated using several leading- 
order (LO) Monte Carlo (MC) event generator programs. These programs are MadGraph |68j . 
CompHep p], and Baur [TOlin]- 




A) u-channel B) l-channel C) s-channel(W channel only) D) Inner Bremsstrahlung 



Figure 9.1: Tree-level diagrams for Z7 and W'y production. 

These programs output 4-vectors and helicities of particles emanating from a diboson produc- 
tion event in an ASCII format. In addition the information on how the particles are produced 
( "mother" and "daughter" ) is recorded, including the energy scale and other parameters used for 
the matrix element calculation. 

These files are then fed into the LesHouchesModule [72], which runs Pythia [73] to add parton 
fragmentation and final-state radiation and initial-state radiation (both QED and QCD) , and 
then writes out the events in CDF HEPG format. These files are then used as input to the CDF 
detector simulation program. This program outputs simulated data in a format identical to that 
of an actual CDF Run II event. Simulated £7 event rates are then estimated in a manner identical 
to that of CDF data. 
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9.1 W-i, Z-i, W-f-i and Z77 MC Sets 



The details on the generator level settings of the Baur [7T] MC 14^7, Z'-j datasets can be found in 
Ref. [71]. 

The MadGraph 14^7, Z7, W'-yy, Z'j'j datasets created for this analysis are listed in Table 19.31 
Details on the MadGraph and CompHep MC samples can be found in Ref. [75l [76] . 

The kinematic cuts used for the generation of W'j, Z7, W'yy and Z77 MC Sets are listed in 
Table O 



Object 


Cuts 


MadGraph and CompHep Samples 




Et, GeV 


V 


AR 


'First' Lepton 


6.0 


4.0 


0.2 


Additional Leptons 


6.0 


4.0 


0.2 


Neutrinos 


1.0 


10.0 





'First' Photon 


6.0 


4.0 


0.2 


Additional Photons 


6.0 


4.0 


0.2 


Baur Samples 




Et, GeV 


V 


AR 


'First' Lepton 


0.0 


10.0 


0.2 


Additional Leptons 


0.0 


10.0 


0.2 


Neutrinos 


0.0 


10.0 





'First' Photon 


5.0 


10.0 


0.2 


minimum m{ii), GeV 




>20 




minimum mliij), GeV 




>20 





Table 9.1: The cuts used at generator level to produce the Z7, W'y, Z77 and 1^77 samples for 
CompHep, MadGraph and Baur datasets. 

To account for NLO corrections to the Wj, Z7, W'-yy, Z'j'j processes we use the E'r-dependent 
K-factor=^^ obtained using Baur's NLO calculation programs [77] . 

We apply these corrections to our LO MC(MadGraph and Baur). For W'yy and Z77 we used 
the same K-factor formulas as for W'j and Z7, respectively. 

Since the Baur NLO program only considers the s, t and u channel contributions and not the 
bremstrahlung off the lepton lines(in this case on the generator level Mw < 76.0, Mz < 86.0) we 
apply the inclusive W cross-section K-factor of 1.36 to the W+photon processes. 

The K-factor applied to Wj and 1^77 MadGraph MC samples is shown in Equation 19.11 
The K-factor applied to Zj and Z'j'j MadGraph MC samples is shown in Equation 19. 2[ In the 
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DataSet Name 


Events 


Cross section (pb) 




Z(ee)7 


395482 


6.855990 




395482 


6.855990 




199047 


6.423524 


W^(ez/)7 


199831 


27.2 




199850 


27.2 




199891 


24.0 


Z77, 1^77 


Z(ee)77 


198830 


0.089137 


Z(/i/i)77 


198830 


0.089137 


Z(rr)77 


198700 


0.078612 


l^(ez/)77 


199351 


0.126 




199351 


0.126 




198910 


0.0939 



Table 9.2: The 1^7, Z7, 1^77 and Z77 MadGraph datasets. 



DataSet Name 


Events 


Cross section (pb) 


Z7, W-f 


Z(ee)7 


429979 


8.62 


Z(/i/i)7 


438468 


8.61 


W^(ez/)7 


140130 


31.9 


Vr(/ii/)7 . 


164732 


31.9 



Table 9.3: The W'^ and Z7 Baur datasets. 



Equations l9.1l and l9.2I MT4/ (Mz) is the mass of the generated W(Z) system , and is a (generated) 
photon transverse energy. 

Mw < 76.0 ^ K-f actor = 1.36 

Mw > 76.0 K-f actor = 1.62 + 0.00001 x - 0.386 x exp(-0.100 x P^) (9.1) 



Mz < 86.0 K-f actor = 1.36 

Mz > 86.0 K-f actor = 1.46 - 0.000728 x P^ - 0.125 x exp(-0.0615 x P}) (9.2) 
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Figure 9.2: The integral cross-sections in fb from MadGraph and COMPHEP at 1.96 GeV versus 
the gamma Et for (a) ez/7 production, (b) e~^e~'j production [75| [76j. 



Every MC event is weighted with the appropriate K-factor. 



9.2 Checks 

To be highly confident in the SM predictions, we compared the predictions from the three inde- 
pendent LO matrix-element generators at generator and HEPG level. 

There is excellent agreement (within 10% or within statistics) between MadGraph and Com- 
pHep in all channels. As the two generators are really different in technique (a helicity amplitude 
calculation in MadGraph, as opposed to the symbolic evaluation of squared matrix element in 
CompHep), this gives us confidence in the predictions. 

Figure 19.21 shows the integrated cross section versus Ex of the photon for the MadGraph and 
CompHep W'j and Z7 samples [75], [76] . Figure 19.31 shows the integrated cross section versus Et 
of the photon for the MadGraph and CompHep W'yy and Z77 samples [751 ES] ■ 

We have compared MadGraph, CompHep and MadGraph samples at GENERATOR [75] and 
HEPG [78]. The more detailed study for MadGraph and Baur Z7 and W'j samples have been 
performed at Ref. [79l [80] . For example, the distributions for the muon channel is shown in 
Figure [931 
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Figure 9.3: The integrated cross section in attobarns (10^'^ fb or 10^^ pb) versus the gamma Et 
for for (a) 14^77 production and (b) Z77 production. The cross section for the highest-Ex photon 
to be above the threshold is in blue; and the 2nd photon is in red [75l [76] . 



9.3 The SM Diboson and Processes as Sources of 
Lepton-Photon Events 

The 1^7 channel was the main SM contributor to the ^7^t signature with the Run I cuts |29j . 
In the Run I analysis, W--^ was expected to contribute 1.93 ± 0.26 events to the 67 channel, 
out of a total of 3.41 ± 0.34, and 1.99 ± 0.27, out of a total of 4.23 ± 0.46 for the /i7^T channel. 
Having a reliable prediction of this signature is crucial. 

The photon can be radiated from the incoming quarks, from the outgoing electron, or the 
intermediate W (Figure [9TTil . More detail on the kinematic distributions and the cross sections is 
available in Ref [75]. Initial-state radiation is simulated by the pythia MC program [73] tuned to 
reproduce the underlying event. The generated particles are then passed through a full detector 
simulation, and these events are then reconstructed with the same code used for the data. 

The predicted numbers of detected 67 and /i7 events in 305j96'^ from SM W-^ production 
satisfying the analysis cuts from MadGraph are given in Table 19.41 

The uncertainties on the SM contributions include those from parton distribution functions 
(5%), factorization scale (2%), K-factor (3%), a comparison of different MC generators (~ 5%), 
and the luminosity (6%) (Chapter [TT]). 

We have studied predictions from MadGraph, CompHep and Baur generators for data and MC 
reconstructed in Ref. [7^, and these studies showed good agreement in predicted rates. 
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Figure 9.4: Comparison of distributions for MadGraph (black) vs Baur(red dots) for Z7 {fi 

channel) after fragmentation (HEPG level). See |79j for details. 



The process pp—^Z^/j* + 7 is also one of the major SM backgrounds for the lepton+ photon+X 
searches. It is the largest contributor to the inclusive multi-body category in the Run I search, 
with expected contributions in 86 pb~^ of 5.01 ± 0.54 events in the ej mode and 4.60±0.54 in the 
/X7 mode [2]. It is also significant in the ^7^t channel, especially for muons, as one muon can 
be missed inducing ^rp. The expected contributions for the electron and muon channels in Run I 
were 0.32 ± 0.5 and 0.96 ± 0.15 events, respectively, smaller than the W'y contributions, but still 
significant, with a total (e + /i) of 1.27 ± 0.17 events out of the 7.64 ± 0.71 events expected in 86 

The photon can be radiated from the incoming quarks or from the outgoing electron (Fig- 
ure [9Tl]). More detail on the kinematic distributions and the cross sections is available in Ref |75j . 
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1 


:S 
O 


V 






o 
o 

a, 

— V 

^= 


Wj, electron channel 


MadGraph 


15.3±0.99 


1.61±0.32 


13.6±0.93 


0±0.042 


0±0.042 


Baur 


15.9±1.31 


2.05±0.47 


13.8±1.21 


0±0.069 


0±0.069 


Average 


15.60±0.82(stat) 


1.83ih0.28(stat) 


13.70±0.76(.tat) 


0.0±0.040(stat) 


0.0±0.040(stat) 




lhl.70(sys) 


±0.48(sy.) 


±1.41{sys) 


ibO.O(sys) 


ihO.O(sys) 


W'y, muon channel 


MadGraph 


10.4±0.77 


1.97±0.34 


8.44±0.70 


0±0.042 


0±0.042 


Baur 


10.8±0.95 


1.56±0.36 


9.25±0.88 


0±0.059 


0±0.059 


Average 


10.60±0.61(stat) 


1.77±0.25(stat) 


8.84±0.56(Btat) 


0.0±0.036(stat) 


0.0±0.036(stat) 




lhl.l8(sys) 


±0.45{Bys) 


lhl.23{sys) 


ihO.O(sys) 


ihO.O(sya) 



Table 9.4: The SM contributions from the W'y channel to the analysis categories. The LO 
calculations from Baur and MadGraph have been corrected by the K-factor and CDF efficiencies. 
The difference between the two LO generators has been included as a systematic uncertainty (see 
text); the uncertainty on the average includes both the statistical uncertainties and this systematic 
uncertainty. 



The predicted numbers of detected (inclusive) 67 and /i7 events in 305pb^^ from SM Z'j pro- 
duction satisfying the analysis cuts from MadGraph are given in Table 19.51 The uncertainties 
on the SM contributions include those from parton distribution functions (7%), a comparison of 
different MC generators (~ 5%), and the luminosity (6%) (Chapter [TT]). 

9.4 The SM Triboson VI/77 and Z77 Processes as Sources 
of Lepton-Photon Events 

While small, the W'-yy, Z77 processes are the largest true SM sources of a signature of a high-px 
lepton plus two photons. The observation of several such events has motivated a careful study of 
the cross-sections for these sources [75j . In this study we have used both MadGraph and Comphep; 
we get excellent agreement between the two generators (Figure [973]) . giving us confidence in the 
predictions. 

The final state of one lepton and two photons, £1/77, is produced in the SM through an inter- 
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1 


:S 





V 









a, 
— V 

^= 


Z7, electron channel 


MadGraph 


21±0.39 


7.77±0.24 


0.97±0.084 


12.1±0.29 


0.12±0.030 


Baur 


22.7±0.43 


8.31±0.26 


1.34±0.11 


12.9±0.33 


0.17±0.038 


Average 


21.85±0.29(stat) 


8.04±0.18(stat) 


1.16±0.069(stat) 


12.50±0.22(stat) 


0.15±0.024(stat) 




lh2.80(sys) 


±0.98(sys) 


±0.39{sys) 


±1.51{sys) 


±0.052(sy.) 


Z7, muon channel 


MadGraph 


15.7±0.32 


3.46±0.15 


4.28±0.17 


7.96±0.22 


0.026±0.013 


Baur 


16.6±0.35 


4.18±0.18 


4.7±0.19 


7.67±0.23 


0.022±0.013 


Average 


16.15±0.24(.tat) 


3.82ih0.12(stat) 


4.49±0.13(stat) 


7.81±0.16(stat) 


0.024±0.01(.tat) 




±1.92(sys) 


lh0.82{sys) 


±0.63{sys) 


lh0.87{sys) 


±0.0047{Bys) 



Table 9.5: The SM contributions from the Z7 electron channel to the analysis categories. The LO 
calculations from Baur and MadGraph have been corrected by the K-factor and CDF efficiencies. 
The difference between the two LO generators has been included as a systematic uncertainty (see 
text); the uncertainty on the average includes both the statistical uncertainties and this systematic 
uncertainty. 



mediate W, with radiation off of any of the charged lines in the diagrams. We denote the final 
state of ^1/77 as W77' for convenience, although the W is virtual and the kinematics of the final 
state are more complicated than the name would suggest. Note that this is process has three 
spin-one bosons in the final state; MadGraph treats the helicities correctly and writes them into 
the output file as input to the next steps. 

The predicted numbers of detected 67 and /i7 events in 305p6~^ from SM 1^77 production 
satisfying the analysis cuts from MadGraph are given in Table I9l6l The agreement on the generator 
level for 1^77 between the two monte carlos, CompHep and MadGraph is good (Figure [9731) . and 
we will use MadGraph 1^77 to get predicted rates. More details, including kinematic distributions, 
can be found in Ref. [75] . 

The final state of two leptons and two photons, ££77, is generated through the intermediate 
photon and Z'^ states, with radiation off of any of the charged lines in the diagrams. We denote 
this as '^'77' for convenience, but the two amplitudes modify the mass spectra and angular 
distributions, and so are both important. 
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C3 








or 





V 




so 

or 

So 

1 




So 

1 


electron channel 


MadGraph 


0.14± 


0.019± 


0.11± 


0.0052± 


0.0067± 




0.0064(stat) 


0.0023(Btat) 


0.0057(stat) 


0.0012(Btat) 


0.0014(stat) 




±0.017(sys) 


±0.0023(sys) 


±0.013(sys) 


±6.2e-04{sy.) 


±8.0e-04(sys) 


M^77, muon channel 


MadGraph 


0.069± 


O.OlOOi 


0.055± 


0.0± 


0.0037± 




0.0042(stat) 


0.0016(stat) 


0.0037(stat) 


1.9e-04(stat) 


9.5e-04{stat) 




±0.0083(sy.) 


±0.0012(sys) 


±0.0066(sy=) 


±0.0(Bys) 


±4.4e-04(sys) 


Z77, electron channel 


MadGraph 


0.54± 


0.26± 


0.029± 


0.24± 


0.015± 




O.OlOO(stat) 


0.0070(stat) 


0.0023(stat) 


0.0067(stat) 


0.0017(stat) 




±0.065(sys) 


±0.031(sys) 


±0.0035(sy.) 


±0.029(syB) 


±0.0018(sys) 


Z77, muon channel 


MadGraph 


0.37± 


0.12± 


0.12± 


0.12± 


O.Olli 




0.0080(stat) 


0.0046(stat) 


0.0046(stat) 


0.0044(stat) 


0.0014(stat) 




±0.044(sys) 


±0.014(sys) 


±0.014(sys) 


±0.014(sys) 


±0.0013(ByB) 



Table 9.6: The predicted number of £7 + X events in 305p6"^ from SM 1^77 and ^'77 produc- 
tion satisfying the analysis cuts from MadGraph. The uncertainty includes both the statistical 
uncertainties and this systematic uncertainty. 



This process is one of the SM mechanisms that could produce the /i/i77j j event [28J (although 
with an extra two jets), and is of interest in the dilepton-diphoton searches as a background. As 
one can see below the SM cross sections are small, typically one femtobarn or less. 

The predicted numbers of detected 67 and /i7 events in 305pb~^ from SM Z77 production 
satisfying the analysis cuts from MadGraph are given in Table l9T6l The agreement on the generator 
level for W'j'j between the two monte carlos, CompHep and MadGraph is good (Figure 19731) . and 
we will use MadGraph 1^77 to get predicted rates. More details, including kinematic distributions, 
can be found in Ref. [75] . 
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9.5 The Sum of Contributions for SM ^7, Vl^77, ^77 
Processes 



Table 19.71 gives the sum of the expected contributions to the 67 and fij channels from SM Wj, 
Z7, W'yy, and ^'77 processes. We have multiplied the average LO predictions for each channel 
by the K-factors listed in Equations 19.11 and 19.21 As we require the event to be triggered either 
by the lepton or by the photon trigger, this combination of triggers is fully efficient. 
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38. 13ih0.87(stat) 


10.15±0.33(stat) 


15.00±0.76(stat) 


12.75±0.22(Btat) 


0.17±0.047(stat) 




±4.58(.ys) 


±1.49(sys) 


lbl.82{sys) 


±1.54(sys) 


±0.055{ByB) 


/i7 


27.19±0.66(stat) 


5.72±0.28(.tat) 


13.51±0.57(stat) 


7.93±0.16(stat) 


0.039±0.037(stat) 




±3.15(sys) 


±1.29(sy.) 


lbl.88{sys) 


lh0.88(sys) 


±0.0064(sys) 


£7 


65.32±1.09(.tat) 


15.87±0.44(stat) 


28.51±0.95(Btat) 


20.68±0.28(Btat) 


0.21±0.060(stat) 




±7.74(=y=) 


lh2.78(sys) 


lh3.70{sys) 


±2.42(sys) 


±0.061{sys) 



Table 9.7: The sum of the expected contributions to the 67 and /X7 channels from SM 14^7, ^7, 
W'yy, and ^'77 processes. The average of the two LO predictions for each channel has been 
multiplied by the K-factors listed in Equations 19.11 and 19.21 



9.6 and Followed by W ^ t u or Z ^ r^r and 

T evv or iivv 

The last SM direct contribution (as opposed to misidentification) we consider is and Z7 
production followed by the boson leptonic decay in the r channel("r7 background"). The tau can 
then decay into an electron or muon. These events are not fakes, in the sense that the electron or 
muon is real, although not a direct product of the vector boson decay. 

Table [23] gives a summary of the contributions to the 67 and /i7 channels from r7 events (yV^ 
and Z7 decaying to taus). Shown are numbers of expected tau events from different processes 
(PF7, Z7) making a signature of 67, /i7 and Z7 for the different categories defined in the analysis. 
The final lepton is either an electron or muon. 
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67 


1.13±0.20(.tat) 


0.42±0.12(=tat) 


0.71±0.16(stat) 


0±0.041(stat) 


0±0.041(.tat) 




±0.14(Bys) 


±0.051(sy.) 


±0.085(.ys) 


±0{sys) 


±0{sys) 




0.66±0.12(stat) 


0.38±0.11(.tat) 


0.26±0.072(stat) 


0.013±0.042(.tat) 


0±0.041(stat) 




±0.079(sys) 


±0.045(sys) 


±0.031(.y.) 


±0.0016(sys) 


±0(sys) 


£7 


1.79±0.24(stat) 


0.80±0.16(.tat) 


0.97±0.18(stat) 


0.013±0.059(.tat) 


0±0.059(stat) 




i0.21(sys) 


±0.096(sys) 


i0.12(sys) 


±0.0016(sys) 


±0(sys) 



Table 9.8: SM contributions from r7 events {W'-) and Z7 decaying to taus). 



Background from W — > tv, where r ^ pu [^1/4 of the r branching fraction), then p — > 
7r7r°, which could mimic a single track + photon signature that looks like £7px ^ P^^t of tj 
background estimated. 



94 



Chapter 10 
Backgrounds: Fakes 



In addition to the expectations from real SM processes that produce real £7 events described in 
Chapter [9l there are backgrounds due to misidentified leptons and photons, and also incorrectly 
calculated ^-p. We generically call these misidentifications 'fakes'. In this chapter we first treat 
backgrounds from fake photons, then from fake leptons, including backgrounds to the W samples 
due to events with a fake lepton and false ^rp. 

10.1 Fake Photons 

We consider two sources of fake photons: QCD jets in which a neutral hadron or photon from 
hadron decay mimics a direct photon, and electron bremsstrahlung, in which an energetic photon 
is radiated off of an electron which is then much lower energy and curls away from the photon. 

10.1.1 Fake Photons from Jets 

High px photons are copiously created from hadron decays in jets initiated by a scattered quark 
or gluon. In particular, mesons such as the tt" or rj decay to photons which may satisfy the photon 
selection criteria. 

The number of lepton- plus-misidentified-jet events in the ^7^t ^^7 samples is determined 
by measuring the jet Et spectrum in £^rp+jet and ii+jet samples, respectively, and then mul- 
tiplying by the probability of a jet being misidentified as a photon, P;J'^*(Et). The uncertainty 
on the number of such events is calculated by again using the measured jet spectrum and the 
upper and lower bounds on the Ex-dependent misidentification rate. An overview of the fake rate 
method is given in Ref. [¥T] . 

Any photon that is due the decay of a meson (tt^ 77, rj 77) is classified as FAKE and 
any photon that is created in the hard scattering process or radiated off a quark is classified as 
a TRUE photon. The strategy is first to measure the RAW fake rate. The fraction of jets which 
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E-r (GeV) E-^ (GeV) 



Figure 10.1: The dN/dE:^^* vs E:^^* distribution for jets in W sample (black points) and in jet 
samples. 2nd high-Exjet in the jet samples is shown in blue squares and the 345th jet in red open 
triangles. 



are matched to a photon candidate 



AT jVfTRUE , jvrFAKE 
D n?jf^i\ ■y- candidate -"7 7 /in 1 

r^RAWKI^T ) - ^7 - ^7 l-LU-iJ 



We estimate TRUE fake as 



j\jFAKE 

PtRUe{^T^) = PrAw(E^t^) X Fqcd, Fqcd = ^TRUE p^FAKE {10.2) 

To distinguish 7 from 7r° or other hadrons, the following variables have been used in the jet 
fake rate studies: 

• CES x' 

- 7r° — i> 77 typically have higher than prompt photons 

• Isolation in a cone in r/ — (/? space of radius R=0.4 around the 7 candidate 

- the background is usually produced as a part of a jet =^ Et is higher than for 7 

• Hit rate in the CPR (Central Preradiator Detector, see Section [2.2.31) 

- CPR is between solenoid and calorimeter. 7 converts in the coil and therefore we measure 
charge. Photons and fakes have different conversion probabilities. 

Technical details on the studies of jets faking photons in CDF II are available in Ref. [81j- The 
most recent numbers are available in Ref. [H] and resulting distribution for Fqcd and Ptrue are 
shown in Figure [10. 2[ We follow this study closely in our estimates. We use jets from the inclusive 
high-px muon and high-Ex electron lepton samples (Section 14.11) to evaluate the background from 
jets faking photons. 



96 




The above studies show that the (fake) photon carries about 94% of the jet Et, and the 
resolution on the photon Et is about 5%. We consequently scale and smear the jet Et by these 
numbers to get the Et of the fake photon. We then select jets with \ri\ < 1.1 and Et > 25 GeV 
after the scaling and smearing. The jets are then weighted with by the jet fake rate [81]. 



1 

or 





V 






c? 


so 



/ 


67 


3.76±3.76(tot) 


0.69 ±0.69{tot) 


2.8 ib2.8(tot) 


0.3 ±0.3(tot) 


0.0003 ±0.0003{tot) 




1.88ihl.88(tot) 


0.14 ±0.14{tot) 


1.6 ±1.6(tot) 


0.2 ±0.2(tot) 


0.00008 ±0.00008{tot) 


£7 


5.64±5.64(tot) 


0.83 ±0.83{tot) 


4.40 ±4.40(tot) 


0.50 ±0.50(tot) 


0.00038 ±0.00038{tot) 



Table 10.1: The predicted backgrounds from jets faking photons in the analysis subcategories. 
We estimate systematic uncertainty on the predicted number to be 100%. Statistical errors are 
negligible. The uncertainties in the electron and muon samples are assumed to be 100% correlated. 



The number of fake events versus the Et of the fake photon are shown in Figure 110.31 for the 
£7^^ and ££7 samples in both the electron and muon channels. We count lepton-'fake photon' 
candidates in the same way as we do for real photons to calculate background for each subcategory. 
Table 110.11 summarizes the number of events in the 67 and /i7 samples for the sub-categories used 
in the analysis. 
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I I Jet Fakes ^ 
;...] Jet Fakes: Upper bound ; 
Jet Fakes: Lower bound - 
Systematic Error is 99% - 




30 40 50 60 70 80 90 100 
Fake Photon 



(a) i'7jErp (Electron Channel) 



I I Jet Fakes 

...] Jet Fakes: Upper bound _: 
'] Jet Fakes: Lower bound 
Systematic Error is 89% 




30 40 50 60 70 80 90 100 
Fake Photon 



(b) (Electron Channel) 



I I Jet Fakes t 
Jet Fakes: Upper bound J 
Jet Fakes: Lower bound : 
Systematic Error is 89% J 




90 100 
Fake Photon 



(c) ^7pT (Muon Channel) 



I I Jet Fakes 

Jet Fakes: Upper bound 
Jet Fakes: Lower bound 
Systematic Error is 92% 




90 100 
Fake Photon 



(d) (Muon Channel) 



Figure 10.3: The Et spectrum of events expected in 305pb with a fake photon from a jet versus 
the Et of the fake photon. The four plots show the number of events per 2 GeV bin expected in 
the •^7pT ^^7 signatures in both the electron and muon channels. The upper and lower error 
bands from Ref. [H] are shown. 



10.1.2 Fake Photons from Electron Bremsstrahlung 

We can measure the probability that a high-Ex electron 'brems' an energetic photon in the material 
before the COT tracking volume in such a way that it fakes a photon, by using the as a source 
of 'tagged' electrons. We look for a back-to-back e+e^ pair close to the mass (A0ee > 150° 
and 81 GeV < M^e < 101 GeV) - this is the sub-category labeled 'Z-like'. For the ej sample the 
requirements are an exactly 1 electron and 1 photon, with A0e7 > 150° and 81 GeV < rrief < 
101 GeV. From the number of these events and the number of Z° e+e~ events we can measure 
the probability per 'leg' of a Z^ that an electron is misidentified as a photon. We require both 
electrons in Z° e+e" to be central and pass tight cuts, so that the kinematic requirements will 
be the same as for the ej sample (i.e. two tight central electromagnetic objects). 

The events in the 2nd row (labeled as Z° —>■ e'^e") of the Table [TO^ have the same signature as 
i'j events (see Figure [T2?T1) . but instead of a photon we require a tight central electron(Table [673]) . 
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For instance, to estimate contribution from electron faking photon to ^7^t category we count 
events with ee+^rp (both electrons are tight central, Table [673]) . For Multi-Lepton + Photon we 
count events with three electrons, of which two electrons should be tight central, and the third one 
can be tight or loose central (Table [6?T1) . or phoenix tight (Table 16^ . Finally, for Multi-Photon 
+ Lepton category we use ee7 events, both electrons are Tight Central. 

Using the numbers in Table 110.21 we calculate the background from electrons misidentified 
as photons. For example, for 67 subcategory the estimated number of electron-fake-photon 
events("e 7 fakes") is calculated as follows: 

/\re7 ]\rdiboson ]\rjets 



jye ^ 7 fakes _ jyZ° 6+6 ^ ZO-like Z^-like ZO-like Qg 

67 ~ ^^i^T /vZ° 6+6" 

ZO~Hke 



N -^T^ ^ ^ = 66 events in the Z° 6+6~ sample 
^zo-iike = ^7 ^°-like events 

^zo-iike = Z^-like events expected from diboson events(iy7, Z'j) 

^zolri, = e+jet faking photon Z°-like events expected from misidentified jets 



' ZO-like 

^zo-J^ke^^^~ = ZO ^ 6+6- events (81 GeV < < 101 GeV, A^^e > 150) 



We take all numbers from data, with the exception that we get number of — like events 
expected from diboson events, W'-f and Z7, from MC (we take into account contribution from 
Z77 and W'-fj, although it's ~ 1% of that from Z7 and W'j). 

Finally, we estimate the number of electron-fake-photon events in ^7^^ to be 

„^ ; . M.. , ,6 , ± 14.45 - 9.03 ± 0.23 - 0.97 ± 0.97 ^ 

ei^T 6174 ^ ' 

Table 110.21 summarizes the calculated number of events in each analysis subcategory for elec- 
trons faking photons by catastrophic bremsstrahlung. The upper three rows are the input numbers 
used in the calculation, which is given by Equation 110. 3[ The last row gives the estimated number 
of events detected with fake photons from electron bremsstrahlung. 



10.2 QCD ('Non-W/Z') Backgrounds 

To measure the QCD backgrounds from fake leptons and or fake ^rp, we form a 'non-W/Z' sample 
we expect to have very little real lepton content [82] by selecting on loose leptons, rejecting events 
from the W or Z. 

To estimate Non-W/Z background we use four samples: 

• Non-W/Z sample: 1 i + jet(s), no W or Z candidates 
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On 




1? 
O 


>-/ 

V 

rC 


4^ 


43 


o 

1 


TP ^ e+e- 


6174 


723 


637 


76 


7 


6 


Diboson 


9.03 ±0.23 


38.13±0.87 


10.15±0.33 


15.00±0.76 


12.75±0.22 


0.17 ±0.047 


Jet fakes 


0.97 ±0.97 


3.76 ± 3.76 


0.69 ±0.69 


2.8 ± 2.8 


0.3 ±0.3 


0.0003 ± 0.0003 


e — > 7 fakes 


199.00± 14.67 


23.30±1.90 


20.53±1.70 


2.45 ±0.33 


0.23 ±0.09 


0.193 ± 0.080 



Table 10.2: Bottom row: The calculated number of events in each analysis subcategory for elec- 
trons faking photons by catastrophic bremsstrahlung. The upper three rows are the input numbers 
used in the calculation, which is given by Equation 110.31 Only statistical errors are quoted. Sys- 
tematic errors estimated by varying Z mass window are found to be negligible. 



• Signal Sample: ^7pT ^^7 

• G olden- Lepton: tight e's, CC Z° —>■ e~^e~] tight /i's —>■ fi~^fi 

• Jet Faking Photon: ij or £ij 

In these samples we define three track isolation regions: 

• Track- Isolated: GeV < Track Isolation < 2 GeV 

• N on- Track- Isolated: Track Isolation >4 GeV 

• Intermediate: 2 GeV < Track Isolation <4 GeV 



The procedure we use is to: 



• Estimate fraction of golden leptons with bad track isolation 

• Estimate fraction of non-W/Z leptons with good track isolation 

• Estimate QCD( Jet faking lepton and ^rpj background 

• Vary track isolation regions to get systematics 

We describe the selection and the procedure in detail below. We define 'Non-W/Z' sample in 
Section [10.2.1[ We describe the basic method in Section [10. 2. 2[ and then modify it to avoid double 
counting and to include systematics in Section [10. 2. 4[ We estimate QCD (non-W/Z) background 
for W^ — > e^v and W^ — >■ /z^z/ in Section [10.2.51 to make sure we obtain results consistent with 
Ref. [Sg. 
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10.2.1 Non-W/Z Sample 



To select an event for the Non-W/Z sample we require the event to have no W or Z candidates. 
Therefore, we require the event to have exactly one tight lepton and no additional loose leptons 
or > 10 GeV. We also require the event to have at least one hadronic jet. For the electron 
Non-W/Z sample the jet is required to have EM Fraction < 0.8, Nn-acks in the jet > 2, E^* > 20 
GeV. For the muon Non-W/Z sample the jet is required to have 0.1 < EM Fraction < 0.9, 

•acks in the jet 

> 2, E^r^' > 20 GeV. 

To reject Z° — > fi~^fi~ events in which one muon is not identified, events with a second track 
with px > 10 GeV and associated EM and HAD calorimeter energies less than 3 and 9 GeV, 
respectively, are rejected (we used the same requirements to select W^ js^u control sample in 
Section EXSl). 
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Figure 10.4: A comparison of distributions from fake 'electrons' from the Non-WZ sample(black 
histogram) and from electrons from tight Z's (red dots). The top set of plots has no (calorimeter) 
IsoEt cut applied; the IsoEt cut is applied in the bottom plots. The plot labeled 'Conversion' 
shows the value of the conversion flag, where the meaning is: - not a conversion, 1 - conversion 
electron, -2 - trident. 



10.2.2 Track Isolation Method 

To calculate the backgrounds from fake W's and fake Z's in which the lepton comes from a jet, 
we use the track isolation of the lepton and the samples of good electrons from Z's and QCD 
background from the non-WZ sample described above. The procedure is as follows: 



101 



Track ISO of mu[cm] | 



Entries 


338 


Overflow 





Entries 


2572 


Overflow 


156 



Muon IsoEt I 




10 12 14 16 18 20 

Track (so of mu[[:m] 



(a) Track Isolation 



Track Iso of mu[cm] | 


Entries 322 
Overflow 




Entries 390 
Overflow 






(d) Track Isolation 



Entries 


338 


Overflow 


8 


Entries 


2572 


Overflow 


2182 



MomenlumTof mu[GeV] | 




(b) Iso Et 



Muon IsoEt I 



Entries 


322 


Overflow 





Entries 


390 


Overflow 







(e) Iso Et 




Entries 


338 


Overflow 





Entries 


2572 


Overflow 






20 30 40 50 60 70 



(c) Muon Pt 

MomentumT of mu[GeV] | 




Entries 


322 


Overflow 





Entries 


390 


Overflow 






20 30 40 50 



70 80 90 100 



Figure 10.5: A comparison of distributions from fake 'muons' 
histogram) and from muons from tight Z's (red dots). The top 
IsoEt cut apphed; the IsoEt cut is apphed in the bottom plots. 
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• Assume W, Z, 14^7, Z7 all have the same underlying event structure, including jets (good 
assumption to first order in the SM, see for example Figure fTO.61) . 

• Define 3 samples: 

— £7pT ^^7 (signal region) 

— golden-lepton (tight central-central electrons TP e+e^, tight muons Z° yU^yU ) 

— Non-W/Z QCD background 

• [0] Ntot'- Number of events in a signal sample 

Count Ntot events of £7^rp(//7) 

• [1] fs'- fraction of golden leptons with bad track isolation 

From the golden-lepton sample find the fraction fs of golden leptons that have bad track 
isolation (tiso>4). 

• [2] Ng'^°^^: number of signal in a tiso > 4 region 

fs X Ntot represents the number of real electrons we will lose by subtracting off electrons 
with tiso>4 in the i'j'^r^^U'y) sample. This has an uncertainty {6fs x Ntot) =^ 

j^UsO, ^j^^ ^^^^ ^ ^ ^^^^ 
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Figure 10.6: Check of the same underlying event structure assumption: muons from fifij (red 
points) vs "golden" Z (black histogram). 



[3] Rb'- ratio of background with good track isolation to the bad track isolation 

From the non-WZ sample find the ratio Rb of background with good track tiso {tiso < 2) 
to those with bad track tiso {tiso > 4). We use this to estimate how many of the tiso<2 
candidates are really QCD background. 



R 



_ Non-WZ^"°<^ 
^ Non-WZ^^''°>'^ 



[4] N^'^"^'^: the number of background with tiso > 4 

From the £7^'p(//7) sample, the number of background with tiso>4 should be the number 
of candidates with tiso>4 minus the expected number of real electrons with tiso>4: 

[5] A'^^'^"^^: the number of background with tiso < 2 
The number of background with tiso < 2 is 

[6] ]\['^<^^^"<'^: the number of background with 2 < tiso < 4 

From number of background in tiso > 4 ]\[^^^°>^ and number of events in Non-WZ sample 
with 2 < tiso < 4 and tiso > 4 we estimate number of background with 2 < tiso < 4. The 
number of background with 2 < tiso < 4 is 

]\T2<tiso<4: _ (]\Ttiso>i\ y Non-WZ^<*'''°<^ 
^^B — V^B I ^ Non-WZi"°>'^ 

[7] Nqcd- QCD background 

Resulting Number of QCD background 

Nqcd = Nj^'°<^ + Ar2<*i^°<4 + p^tiso>4 
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10.2.3 Non-W/Z Results 

We have used track isolation to make estimates of the non-W backgrounds in the —>■ l^u 
control sample (Section 110.2.51) . 

Table 110.31 summarizes the numbers of track-isolated and non-track-isolated events in the good 
electron, non-WZ, and signal samples. In addition, we add numbers of events from W+jet and 
Z+jet categories, used to estimate jet faking photon background. 

The good electrons are heavily track-isolated (0 < tiso < 2); the QCD background 'electrons' 
are predominantly track- non-isolated (4 < tiso). 



GeV < tiso < 2 GeV I 2 GeV < tiso < 4 GeV I tiso > 4 GeV I total 



Electrons: 



^tight 


11204 


610 


534 


12348 


Non-WZ 


1179 


870 


3709 


5758 




21 


2 


2 


25 


ee7 


35 


1 


2 


38 




1.98655 


0.164143 


0.62847 


2.77916 


eej 


0.46559 


0.0241221 


0.0371426 


0.526854 



Muons: 



Ztight 


5686 


190 


36 


5912 


Non-WZ 


173 


78 


93 


344 


/^7pT 


16 


1 





17 




22 


2 





24 




1.41881 


0.0601055 


0.0618934 


1.54081 


m 


0.264149 


0.00932185 


0.0318951 


0.305366 



Table 10.3: The numbers of events that are track-isolated (0 < tiso < 2) versus non-track-isolated 
{tiso > 4) for good electrons, QCD background 'electrons', and the signal samples. The good 
electrons are heavily track-isolated; the QCD background 'electrons' are predominantly track- 
non-isolated. 



Table 110.41 shows step-by-step calculation of predicted QCD (Non-WZ) backgrounds. 



10.2.4 Modified Track Isolation Method 

From the jet fake rate we obtain (fake gamma -|- real W) and (fake gamma -|- fake W) backgrounds. 
From Track Isolation Method we obtain (real gamma -|- fake W) and possibly (fake gamma + fake 
W). 

To summarize, there are three sources of faking W'y with either fake gamma and/or fake W: 
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Sample 


ftiso>4 

Js 


]ytiso>4 


Rb 


]\Ttiso>4 


]\Ttiso<2 
^^B 


■\T2<tiso<4 
^^B 




67 

667 

/X7^T 


0.043±0.0023 
0.043±0.0023 
0.0061±0.0011 
0.0061±0.0011 


1.1±0.057 
1.6±0.086 

0.1±0.019 
0.15±0.026 


0.32±0.014 
0.32±0.014 
1.9±0.33 
1.9±0.33 


0.92±0.057 
0.36±0.086 
0±0.019 
0±0.026 


0.29±0.031 
0.11±0.032 

0d=0.035 
0±0.049 


0.22±0.024 
0.084±0.024 
0±0.016 
0±0.022 


1.4±0.11 
0.55±0.14 
0±0.069 
0±0.097 


Table 10.4: Predicted QCD backgrounds: step-by-step calculation. 


Sample 


ftiso>4: 

Js 


T\jtiso>4 
^^S 


Rb 


i\jtiso>4 
^^B 


'\Ttiso<2 
^^B 


]\.T2<tiso<4 
^^B 


^tot 


67 

667 


0.043±0.0023 
0.043±0.0023 
0.0061±0.0011 
0.0061±0.0011 


0.96±0.05 
1.6±0.085 
0.094±0.017 
0.14±0.026 


0.32±0.014 
0.32±0.014 
1.9±0.33 
1.9±0.33 


0.41±0.05 
0.34±0.085 
0±0.017 
0±0.026 


0.13±0.022 
0.11±0.032 
0±0.031 
0±0.048 


0.096±0.017 
0.08±0.024 
0±0.014 
0±0.022 


0.64±0.089 
0.53±0.14 
0±0.063 
0±0.096 


Table 10.5: Predicted QCD backgrounds: with jet fakes subtracted. Step-by-step calculation. 


Sample 


ftiso>4: 

Js 


]\ltiso>4 


Rb 


j\^tiso>4 


j\qtiso<2 


jY-2<tiso<4 


^tot 


67 
667 


0.043±0.0023 
0.043±0.0023 
0.0061±0.0011 
0.0061±0.0011 


0.84±0.044 
1.6±0.084 
0.085±0.015 
0.14±0.026 


0.32±0.014 
0.32±0.014 
1.9±0.33 
1.9±0.33 


0±0.044 
0.33±0.084 
0±0.015 
0±0.026 


0±0.014 
0.1±0.031 
0±0.028 
0±0.048 


OiO.Ol 
0.077±0.023 
0±0.013 
0±0.021 


0±0.068 
0.51±0.14 
0±0.056 
0±0.095 



Table 10.6: Predicted QCD backgrounds: with jet fakes double-subtracted. Step-by-step calculation. 



1. fake gamma + real W (from jet fake rate) 

2. fake gamma + fake W (also part of jet fake rate and possibly part of non-W/Z background) 

3. real gamma + fake W (from non-W/Z background) 

If fake W is a part of non-W/Z background it should contribute to track- non-isolated regions. 
Real W's from W+jet should have the same track isolation distribution as W+7 (we have per- 
formed the checks for all the samples we use, see for example Figure [T0.6p . 

To avoid double counting, we modify Track Isolation Method as follows: we subtract W+jet 
faking photon contribution from ^7^x trackiso regions, we subtract Z+jet faking photon contri- 
bution from ££7 trackiso regions. Then we repeat procedure, documented in Section 110.2.21 with 
these modified £7pT ££7 numbers (Table [T03|) . 

To take into account systematic error of 100% for jet faking photon rate (Table [TO. 51) we repeat 
these procedure, subtracting W+jet and Z+jet contribution multiphed by 2. Therefore, to avoid 
double-counting we use background estimates from "Predicted QCD Backgrounds with Jet Fakes 
Subtracted" (Table [103]). 



Region, GeV 


e7^T 


667 






0-2; 2-00 


1.4+0.11 


0+0.17 


0.58+0.12 


2.1+0.23 


0-2; 4-cx) 


0.64+0.09 


0.53+0.14 


0+0.06 


0+0.10 


0-2; 5-00 


1.2+0.09 


1.2+0.14 


0+0.07 


0+0.1 


0-2.5; 4.5-00 


0.93+0.09 


0.9+0.14 


0+0.07 


0+0.1 


0-3; 3-00 


0.03+0.08 


0+0.14 


0+0.08 


1.6+0.2 


0-3; 5-00 


1.2+0.09 


1.2+0.14 


0+0.07 


0+0.1 


Final 


0.7 


0.6 


0.3 


1.0 




+ 0.1(stat) 


+ 0.1(stat) 


+ 0.1(stat) 


+ 0.1(stat) 




±0.7(sys) 


+ 0.6(sys) 


+ 0.3(sys) 


+ 1.0(sys) 



Table 10.7: QCD studies table: estimating the QCD background faking a £7^t ^^7 event. The 
method uses track isolation and two regions, one at low track-iso and one at high. Leptons from 
the Z° —>■ sample are used to estimate the number of leptons in the high-track-iso region; 
'leptons' from the non-WZ sample are used to estimate the number of leptons in the low-track-iso 
region. The first column of the table gives the ranges in GeV for the low and high track-iso regions, 
respectively. The last line gives the final estimates. 

Finally, Non-W/Z backgrounds for ^7^x ^'^'^ ^^7 signatures are summarized in Table fT0.2.4[ 
Systematic errors are obtained by varying track isolation regions. Further checks for Non-WZ 
background for W^ — > 6^z/ and W^ — > fi'^u are described below in Section [10.2.5[ 
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10.2.5 Non-W/Z background for ^ e^u and fi^u 

We have used track isolation to make estimates of the non-WZ backgrounds in the — > /^z/ 
control sample and iij and £7^t signal samples. Table fTO^ shows the estimated 'non-W back- 
ground in the W samples for 5 different track-iso regions. This is for a check- the values are 
consistent with known QCD backgrounds for W's [62]. The method uses track isolation and two 
regions, one at low track-iso and one at high. Leptons from the golden Z sample are used to 
estimate the number of leptons in the high-track-iso region; 'leptons' from the non-WZ sample 
are used to estimate the number of leptons in the low-track-iso region. The first column of the 
table gives the ranges in GeV for the low and high track-iso regions, respectively. 



Re^ 


^ion, GeV 


W± ^ e^u 


W± ^ fi^u 


0-2 


4- 


oo 


6.5e+02±6.6e+02 


1.8e+03±7.3e+02 


0-3 


5- 


oo 


l.le+03±6.5e+02 


1.6e+03±7.2e+02 


0-2 


5- 


oo 


l.le+03±6.5e+02 


1.6e+03±7.3e+02 


0-2 


2- 


oo 


56±8.3e+02 


5.7e+02±7.7e+02 


0-3 


3- 


oo 


4.5e+02±7.1e+02 


l.le+03±7e+02 



Table 10.8: QCD backgrounds for W's The method uses track isolation and two regions, one at 
low track-iso and one at high. Leptons from the golden Z sample are used to estimate the number 
of leptons in the high-track-iso region; 'leptons' from the non-WZ sample are used to estimate the 
number of leptons in the low-track-iso region. The first column of the table gives the ranges in 
GeV for the low and high track-iso regions, respectively. 
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Chapter 11 

Systematic Uncertainties 



In this chapter we summarize estimates of the systematic uncertainties on the SM predicted rates 
and on the measured event counts. 

The errors are categorised as theoretical (Section lll.ip . luminosity (Section 111.21) and experi- 
mental (Section 111.31 The contributing effects for the SM predictions we have considered are: 

• 7% error is on the total theoretical prediction, including the NLO uncertainties. 

• Luminosity: 6% 

• Trigger Efficiencies: 2% for muons and 1% for electrons for lepton triggers only. We OR'ing 
lepton trigger with photon trigger, and therefore this combination of triggers is fully efficient. 

• |z_vert| < 60: 1% 

• Muon ID Efficiencies: 2% 

• Electron ID Efficiencies: 1% 

• Photons ID Efficiencies: 4% 

The systematic uncertainties on the backgrounds are included in the background estimates, 
discussed in Chapter |9] and Chapter [TOl For the SM predictions the total systematic uncertainty 
is 10.2% for W'y and Z7 for electrons, and 10.5% for W'y and Z7 for muons (Chapter [9]). 

11.1 Theoretical Systematic Uncertainties 

Limitations in the theoretical precision of the calculation, result in an uncertainty on the cross- 
section prediction. The effect of these errors on the cross-section is studied in [71 ST], [H3] and is 
summarized in Table [TLT] 

11.1.1 Factorization Scale 

The factorization scale is the minimum value calculated [CT [83] for photon emission in the 
ZGAMMA and WGAMMA programs [70]. This value will affect the maximum value for post 
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Source 


% 


Factorization Scale 


2 


PDF 


6 


K-factor 


3 


Total 


7 



Table 11.1: Systematic errors on the ^7, W'f, ^77 and W'y'y generation 

generation Pythia fragmentation. The default factorization scale was s, the square of the collision 
energy of the event. The cross-section and acceptance were measured using four other values, 2 
s, 3/2 s, 2/3 sand 1/2 s. The greatest variation in the cross-section from the default value of Q2 
was 2%. 

11.1.2 Parton Density Function Choice 

Protons and anti-protons are composite particles. Therefore, any interactions between them must 
be described using parton density functions (PDF). The PDF describes the energy distributions 
of the valence quarks, gluons and sea quarks inside the proton/anti-proton. 

The PDF chosen for use with ZGAMMA [70] was the CTEQ5L PDF. In order to determine 
the systematic error from this choice, the LO cross-section is compared to the corresponding 
predictions calculated from the MRST 72 - 76 PDFs. The MRST cross-sections range between 

I. 604 and 1.625 pb^^ whereas the cross-section using CTEQ5L is 1.72 pb~^. The difference between 
the two was taken to be the systematic error [HI [83], of 6%. 

11.1.3 K-factor 

The calculated K-factor only takes into account 0{as) corrections. To take into account higher 
order corrections, the Q scale in the NLO calculation was varied by factors of 2 and 1/2. A 3% 
variation in the cross-section calculation was observed, and taken to be a systematic error. 

II. 2 Luminosity Systematic Uncertainties 

The luminosity error is estimated to be 6%, which includes a 4.4% contribution from the accep- 
tance and operation of the luminosity monitor and 4.0% from the theoretical uncertainty on the 
calculation of the total pp cross-section [19] . 



109 



11.3 Experimental Systematic Uncertainties 



The sources of experimental systematic errors [3 HUES] for the + X analysis subcategories are 
summarized in Table lll.2[ Jet Fake systematic error is discussed in Section llO.l.ll 



Source 


% 


Central 


Plug 


CMUP 


CMX 


Jet Fake 


^100 


X 


X 


X 


X 


Zq cut efF 


1.0 


X 


X 


X 


X 


photon cut eff 


2.0 


X 


X 


X 


X 


energy scale (7) 


3.0 


X 


X 


X 


X 


conversion rate uncertainty 


1.5 


X 


X 


X 


X 


momentum scale (a) 


2.0 






X 


X 


acceptance (e) 


1.0 


X 


X 






acceptance (/i) 


2.0 






X 


X 


central e ID 


1.0 


X 








central e trigger 


1.0 


X 








energy scale (e) 


1.0 


X 


X 






plue e ID 


2.5 




X 






plug trig eff 


1.0 




X 






plug e vertex eff 


1.0 




X 






plug e track eff 


1.5 




X 






cosmic 


0.01 






X 


X 


Cot track reconstruction 


0.4 


X 




X 


X 


CMUP ID 


0.7 






X 




CMUP reconstruction 


0.6 






X 




CMUP trigger 


0.7 






X 




CMX ID 


0.8 








X 


CMX reconstruction 


0.3 








X 


CMX trigger 


0.6 








X 



Table 11.2: Systematic errors summary for £7. 'x' means that channel needs to take into account 
its systematic uncertainty. Jet Fake systematic error is discussed in Section llO.l.ll 

Systematic uncertainty on jet faking photon rate is one of the dominating errors. The uncer- 
tainty is limited by the statistics for the high-Ex photons, so we expect it to significantly improve 
with more data. At that point we'll be dominated by the systematic uncertainty. Therefore, the 
biggest contribution will be from the SM estimates on Wj, Zj, W'-ff and ^'77 production. 
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Chapter 12 

The £7 + X Search 



This chapter presents the results of the £7 + X search for the three signatures of interest - ^7pT? 
U-i and ££7. 

Section [12.11 describes the 'analysis subcategories' established in the Run I analysis [2], and 
used again here so as to be a priori. Section 112.21 presents the number of events in each analysis 
subcategory for the 67 + X, /i7 + X, and ^7 + X (the sum of e + /i) samples. Section [12.31 discusses 
the stability of the observed numbers during the course of the run. The predicted and observed 
totals for the ^7^x5 ^^7 ^^"^ ^77 ^"^^ comparison of observed kinematic distributions to the SM 
predicted shapes is done in Section [1 2. 4[ 

12.1 Defining the Event Categories by Topology 

Categories of photon-lepton events were defined a priori in a way that characterized the different 
possibilities for new physics. For each category, the inclusive event total and basic kinematic 
distributions can be compared with standard model expectations. The decay products of massive 
particles are typically isolated from other particles, and possess large transverse momentum and 
low rapidity. 

Therefore, inclusive £7 events are selected by requiring a central tight photon with Ei^ > 25 
GeV and a central e or fi with > 25 GeV. Both signal and control samples are drawn from 
this sample (Figures [HJl [12ll] and MM- 

Considering the control samples first, from the £7 sample we select back-to-back events with 
exactly one photon and one lepton (i.e. < 25 GeV); this is the dominant contribution to the £7 
sample, and has a large Drell-Yan component. A subset of this sample is the 'Z-like' sample, which 
provides the calibration for the probability that an electron radiates and is detected as a photon, 
as discussed in Section [l0.1.2[ The remaining back-to-back events are called the Two-Body Events 
and were not used in this analysis. 

All events which either have more than one lepton or photon, or in which the lepton and 
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photon are not back-to-back (and hence the event cannot be a Two-Body event), are classified 
as 'Inclusive Multi-Body i'j + X\ These are further subdivided into three categories: ^7^t 
(Section 112.51) ('Multi-Body ^7^t Events'), for which the ^-p (Section 18.11) is greater than 25 
GeV , i£-f (Section [nSD and ^77 (Section [I22D ('Muhi-Photon and Multi-Lepton Events'), and 
events with exactly one lepton and exactly one photon, which are not back-to-back. The events 
with exactly one lepton and exactly one photon, which are not back-to-back were not used in the 
analysis. 

12.2 The Number of Events Observed 

Figure fT2. II shows the results of this classification for the inclusive electron data sample. We find 
508 67 events, of which 111 are in the Inclusive Multi-Body category. Of these, 25 are classified 
as 67 events and and 19 as Multi-Photon and Multi-Electron events respectively. 

Figure [T2.2I shows the results for the inclusive muon sample. We find 66 n'j events, of which 41 
are in the Inclusive Multi-Body category. Of these, 17 are classified as fi'j^T events and and 12 
as Multi-Photon and Multi-Muon events. 

Fi gure [T2 . 3 1 shows the sum of the electron and muon entries in the analysis subcategories. There 
are 42 /7^t events, Multi-Photon events, and 31 Multi-Lepton events. It is these categories, 
shown in red in the figures, that are of particular interest due to the Run I results. 
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Electron-Photon Sample 
1 Electron and 1 Photon 
Et > 25 GeV 
508 Events 



Exactly 1 Electron 
Exactly 1 Photon 

Acpe-y > 150 

< 25 
397 Events 



Inclusive Multi-Body Events 
(All Other Electron-Photon) 
111 Events 



Z-Like Electron-Photon 
81 Gev < Me^ < 101 Gev 
(Background Calibration) 
209 Events 



Exactly 1 Electron 
Exactly 1 Photon 

A0e7 < 150 

< 25 GeV 
67 Events 





Muhi-Body ej'^T^ 


Multi-Photon and 


Two-Body Events 


Events 


Multi-Electron Events 


188 Events 


;Et > 25 GeV 


and 19 Events, 




25 Events 


resp. 



Figure 12.1: Electron-photon sample: the subsets of inclusive 67 events analyzed. 
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Muon-Photon Sample 
1 Muon and 1 Photon 
Et > 25 GeV 
66 Events 



Exactly 1 Muon 
Exactly 1 Photon 
A0^^ > 150 

< 25 
25 Events 



Inclusive Multi-Body Events 
(All Other Muon-Photon) 
41 Events 



Z-Like Muon-Photon 
81 Gev < M^^ < 101 Gev 
(Background Calibration) 
9 Events 



Exactly 1 Muon 
Exactly 1 Photon 
A0^.^ < 150 
< 25 GeV 
12 Events 





Multi-Body fi-f'^^ 


Multi-Photon and 


Two-Body Events 


Events 


Multi-Muon Events 


16 Events 


> 25 GeV 


and 12 Events, 




17 Events 


resp. 



ure 12.2: Muon-photon sample: the subsets of inclusive /X7 events analyzed. 



114 



Lepton-Photon Sample 
1 Lepton and 1 Photon 
Et > 25 GeV 
574 Events 



Exactly 1 Lepton 
Exactly 1 Photon 

A^i^ > 150 
< 25 

422 Events 



Inclusive Multi-Body Events 
(All Other Lepton-Photon) 
152 Events 



Z-Like Lepton-Photon 
81 Gev < Me^ < 101 Gev 
(Background Calibration) 
218 Events 



Exactly 1 Lepton 
Exactly 1 Photon 
A(f)i^ < 150 

< 25 GeV 
79 Events 





Multi-Body i-f'^r^ 


Multi-Photon and 


Two-Body Events 


Events 


Multi-Lepton Events 


204 Events 


> 25 GeV 


and 31 Events, 




42 Events 


resp. 



ure 12.3: Lepton-photon sample: the subsets of inclusive £7 events analyzed. 
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12.3 Stability of the Event Rates versus Run Number 



This sub-section looks at the rate for the analysis subcategories as a function of run number. We 
use the same eight luminosity bins used to check the stabihty of the control samples, described in 
Section 14.11 We see no obvious problems, although the statistics are low. Figure 112.41 shows the 
rates in events/p6~^ for the ££7 and ^7^t signal subcategories in each run segment in the muon 
channel; Figure 112.51 does the same for the electron channel. 



0.4- 
0.3^ 

0.2- 
0.1- 



1450 1500 1550 1600 1650 1700 1750 1800 1350 



1450 1500 1550 1600 1650 1700 1750 1800 1850 



xlO^ 



(a) U-i 



(b) l^^^ 



Figure 12.4: Stability plots for the muon channels of the rate in events per pb ^ for the 8 run 
segments (see Section WA\ for: a) ££7, and b) ^7^t- 



Entries vs Runs for stability Entries 
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/ ndf 
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Entries vs Runs for stability Entries 



Overflow 

z'/ndf 



4+- 



1450 1500 1550 1600 1650 1700 1750 1800 1850 



3x10^ 



(a) «7 



(b) ^7^, 



Figure 12.5: Stability plots for the electron channels of the rate in events per pb ^ for the 8 run 
segments (see Section WA\\ for: a) ££7, and b) ^7^t- 
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12.4 Results 



In this section we summarize the predicted and observed totals for the ^7^x5 ^^7 ^11 searches. 
We compare observed kinematic distributions to the SM predicted shapes. 

12.5 £7Pt Search 

The predicted and observed totals for the ^7^t search are shown in Table 112.11 We observe 42 
£7^rp events compared to the expectation of 37.3 ±5.4 events. 



Lepton+Photon+^rp Events 



SM Source 


e7^T 




(e + /i)7;pT 




13.70±1.89 


8.84±1.35 


22.54±2.80 


Z0/7* + 7 


1.16±0.40 


4.49±0.64 


5.65±1.03 


W^±77,Z°/7*+77 


0.14±0.02 


0.18±0.02 


0.32±0.03 


iy±7,Z°/7* +7-^x7 


0.71±0.18 


0.26±0.08 


0.97±0.22 


W^+Jet faking 7 


2.8±2.8 


1.6±1.6 


4.4±4.4 


Z^/y* — >e^e~, e^j 


2.45±0.33 




2.45±0.33 


Jets faking £ + '^j, 


0.7±0.7 


0.3±0.3 


1.0±0.8 


Total SM 








Prediction 


21.7±3.4 


15.7±2.2 


37.3±5.4 


Observed 








in Data 


25 


17 


42 



Table 12.1: A comparison of the numbers of events predicted by the standard model(SM) and 
the observations for the ^7pT search. The SM predictions for the search are dominated by Wj 
production, respectively [68l [711 EH]. Other contributions come from Z7 production, from the 
tri-boson processes W'yy and ^'77, leptonic r decays, and misidentified leptons, photons, or ^rp. 

There is no significant excess in the ^7^t signature. Figure [12.61 shows the observed distri- 
butions summed over the 67 and /i7^T events in a) the Et of the photon; b) the Et of the 
lepton; c) the missing transverse energy, ^rp; and d) the transverse mass of the ^7^t system, 
where Mt = [(E^^ + + ^tY - (^t + ^? + ^r)T^'- 

The predicted and observed kinematic distributions for /i7^x compared in Figure [TZTl 
The distributions for 67 signature are compared in Figure [T2T81 
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Figure 12.6: Distributions for the events in the ^7^^ sample (points) in a) the Et of the photon; b) 
the Et of the lepton; c) the missing transverse energy, prp; and d) the transverse mass of the ^7^t 
system. The histograms show the expected SM contributions, including estimated backgrounds 
from misidentified photons and leptons. 

The additional plots of the identification variables for 67^^ and /i7^T available in Sec- 
tion ID 



118 





14 = 




12: 


> 


10: 


(5 




o 


8: 






w 




c 


6- 






UJ 


4: 




2 




of- 



I I ' I I I ' ' I ' I ' I ' ' ' I ' ' ' I 
CDF Run II (a) 

• Data(\i), 305 pb'^ 
I I IVy 



l///ef, VI, QCD, Zyy, Wyy 




± 



40 60 80 100 120 140 
Photon Et (GeV) 




20 40 60 80 100 120 140 160 180 200 
Lepton Et (GeV) 




80 100 120 140 
^^(GeV) 



100 150 200 
(GeV) 



300 



Figure 12.7: Distributions for the events in the Ai7^T sample (points) in a) the Et of the photon; b) 
the Et of the muon; c) the missing transverse energy, ^t! d) the transverse mass of the •^7^t 
system. The histograms show the expected SM contributions, including estimated backgrounds 
from misidentified photons and leptons. 
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Figure 12.8: Distributions for the events in the 67 sample (points) in a) the Et of the photon; 
b) the Et of the electron; c) the missing transverse energy, ^t! d) the transverse mass 
of the ^7pT system. The histograms show the expected SM contributions, including estimated 
backgrounds from misidentified photons and leptons. 
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12.6 U-i Search 



The predicted and observed totals for the ££7 search are shown in Table [1221 We observe 31 ££7 
events compared to the expectation of 23.0 ± 2.7 events. 



Multi-Lepton + Photon Events 


SM Source 


ee7 


/i/i7 


//7 


Z0/7* + 7 


12.50±1.53 


7.81±0.88 


20.31±2.40 


^°/7* + 77 


0.24±0.03 


0.12±0.02 


0.36±0.04 


Z°/7*+Jet faking 7 


0.3±0.3 


0.2±0.2 


0.5±0.5 


Z°/7*^e+e-,e^7 


0.23±0.09 




0.23±0.09 


Jets faking £ + 


0.6±0.6 


l.Oil.O 


1.6±1.2 


Total SM 








Prediction 


13.9±1.7 


9.1±1.4 


23.0±2.7 


Observed 








in Data 


19 


12 


31 



Table 12.2: A comparison of the numbers of events predicted by the standard model(SM) and 
the observations for the ££7 search. The SM predictions for the search are dominated by Z7 pro- 
duction [681 ITU [69] . Other contributions come from the tri-boson process 2^77, and misidentified 
leptons or photons. 

The ££7 search criteria select 31 events (19 ee7 and 12 /i/i7) of the 574 ^7 events. No 6/17 events 
are observed. Figure 112.91 shows the observed distributions in a) the Et of the photon; b) the 
Et of the leptons; c) the 2-body mass of the dilepton system; and d) the 3-body mass mu^. For 
the Z7 process occurring via initial state radiation, the dilepton invariant mass irtu distribution 
is peaked around the Z^-pole. For the final state radiation, the three body invariant mass mu-y 
distribution is peaked about the Z°-pole. 

The predicted and observed kinematic distributions for |UyU7 are compared in Figure [T2. 101 The 
distributions for ee7 signature are compared in Figure 112. 11[ The dominant contribution for the 
ee7 and /i/i7 signatures is from the SM Z7 production. 

We do not expect events with large the ££7 sample, based on the SM backgrounds; the 

Run I ee77 event was of special interest in the context of supersymmetry [TH [19] due to the 
large value of (55 ± 7 GeV). Figure [r2.12l shows the distributions in for the /i/i7 and ee7 
subsamples of the ££7 sample. No events are observed with > 25 GeV . 

The additional plots of the identification variables for ee7 and /i/i7 are available in Section IA.31 
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Figure 12.9: Distributions for the events in the ii'j sample (points) in a) the Ex of the photon; 
b) the Et of the leptons (two entries per event); c) the 2-body mass of the dilepton system; and 
d) the 3-body mass mu^. The histograms show the expected SM contributions. 
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Figure 12.10: Distributions for the events in the sample (points) in a) the Et of the photon; 
b) the Et of the muons (two entries per event); c) the 2-body mass of the dimuon system; and d) 
the 3-body mass m^^-y. The histograms show the expected SM contributions. 
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Figure 12.11: Distributions for the events in the ee^ sample (points) in a) the Ex of the photon; 
b) the Et of the electrons (two entries per event); c) the 2-body mass of the dielectron system; 
and d) the 3-body mass rriee^. The histograms show the expected SM contributions. 




Figure 12.12: Distributions in missing transverse energy observed in the inclusive search for 
a) /i/i'j events and b) ecy events. The histograms show the expected SM contributions. 
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12.7 ^77 Search 



The predicted and observed totals for the £77 search are shown in Table 112. 31 We do not observe 
any £77 candidate events compared to the expectation of 0.23 ± 0.080 events. 



Multi-Photon + Lepton Predicted Events 


SM Source 


677 


/i77 


hi 




0.0067 ±0.0014 


0.0037 ±0.00095 


0.010 ±0.0017 




0.015 ±0.0017 


0.011 ±0.0014 


0.026 ± 0.0022 


Z°7, e^7 


0.193 ±0.080 




0.193 ±0.080 


Total SM 








Prediction 


0.22 ±0.080 


0.015 ±0.0022 


0.23 ±0.080 


Observed 








in Data 












Table 12.3: The predicted number of multi-body events with additional photons in 305p6~^ from 
Standard Model sources, and the numbers of events observed. The Z77 and 1^77 predictions are 
from MadGraph. The "2^°7, e 7" prediction is taken from the data. 
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Conclusions 



In conclusion, we have searched for physics beyond the Standard Model in a channel in which 
new phenomena due to supcrsymmetry or extra dimensions, for example, could well appear. In 
particular, we have repeated the Run I CDF search for inclusive lepton + photon production, a 
final state in which both a very rare event appeared and also in which there seemed to be an excess 
over SM predictions. The new analysis, the subject of this thesis, was done specifically with the 
same kinematic requirements as the Run I search, but with a significantly larger data sample and 
a higher collision energy. 

We conclude that the excess in the ^7^t signature in Run I of 2.7 sigma was at least largely 
a statistical fluctuation. If the Run I ratio of the number of ^7jEx events observed to the number 
expected, 16/7.6, had held up, the "2.7 o excess" in this channel observed in Run 1 would have 
resulted in an observation of 78±11 events when 37.3 ± 5.4 are expected in the Run II repeat of 
the analysis, versus the 42 events observed. 

We find that the numbers of events in the i?7^T ££7 subsamples of the £7 + X sample 
agree with the SM predictions. We find no events like the 6677 candidate event of Run I, and, 
even more generally, observe no il^ events with anomalous large or with multiple photons. 
We have no explanation for the Run I 6677 event, and nothing we have measured leads us to 
believe that is background. 

However, we still find that in the ^7^t signature in the Run II data the number of observed 
events is higher than predicted, although the excess is now shght. It is possible that the leading- 
order + K-factor theoretical calculations of the diboson 1^7 and Z7 channels, which contain many 
diagrams including initial state radiation, are not precise enough for the precision we have now 
reached. 

Alternatively, or perhaps in addition, this analysis has observed a small number of events on the 
'tails' of the kinematic distributions, in regions we expect few SM events. These events contribute 
to the observation of more events than expected in the •^7pT signature, much as the 6677 
event contributed to the excess in the Run I search. Whether these are very rare backgrounds or 
something new will require yet more data. 

The Fermilab plan is to have a factor of 10-20 more data than presented here by the end of Run 
II of the Tevatron. The increased statistics will require an improved understanding of backgrounds 
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as well as better SM predictions. In particular the estimate of the rate for a jet to be misidentified 
as a photon is limited now by the statistics for high-Ex photons; we expect the estimate will 
significantly improve with more data. 

In summary, while it would have been very exciting to find physics beyond the Standard Model, 
we found no more 6677 events in a much larger sample than in Run I, and the Run I excess in 
£7^^ became less significant rather than more. However, we have conclusively settled a question 
that generated much interest in the theoretical community. The channels we have investigated 
will remain interesting, and the techniques we have developed and the knowledge gained will be 
useful for similar searches at the LHC. 
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Appendix A 
Appendices 



A.l List Of Lepton-Photon Events 





7 




^7 




run/event 


Ex (p Tj 


type PT r] T 






179043/12444345 


+29.13 -0.88 0.07 


CMUP 27.16 -2.18 -0.20 1 


1.30 34.60 


27.67 1.36 


178602/5535619 


25.97 -1.29 1.04 


CMUP 33.03 0.29 -0.24 1 


1.58 58.37 


30.98 2.96 


178816/108292 


+31.28 1.78 0.77 


CMUP 77.49 -1.35 -0.20 1 


3.12 110.7 


34.17 3.06 


178855/422245 


+31.27 0.48 0.57 


CMX 42.14 0.05 0.91 1 


0.43 17.13 


31.94 2.39 


151843/1584392 


+25.40 -1.55 0.10 


CMUP 41.13 3.12 -0.21 1 


1.61 47.87 


48.05 0.76 


151870/1255798 


+79.18 -0.24 -0.28 


CMX 62.62 2.32 0.80 1 


2.56 158.8 


27.78 -1.77 


153739/175671 


+28.20 -1.41 0.78 


CMUP 41.18 -1.26 0.26 1 


0.15 16.61 


32.21 -0.42 


155895/G33C800 


+55.41 2.30 -0.29 


CAIUP 52.70 -O.OG 0.24 1 


2.35 104.5 


31.93 -1.76 


160230/4222557 


+52.97 1.22 1.02 


CMX 49.38 0.83 0.66 1 


0.39 25.13 


43.43 -2.12 


162686/2327952 


+50.31 0.10 -0.23 


CMUP 79.31 0.86 0.33 1 


0.75 58.83 


44.12 -2.52 


166406/10446136 


+119.8 0.86 0.05 


CMX 25.17 -2.62 0.78 1 


2.81 114.3 


110.2 -2.26 


166653/3270001 


+32.07 -0.03 0.43 


CMUP 26.39 0.63 0.51 1 


0.66 18.81 


48.83 -2.89 


164274/2876183 


+35.10 -1.17 0.19 


CMX 38.78 -2.49 0.70 1 


1.31 49.92 


45.23 1.37 


166008/3466824 


42.82 -0.13 0.90 


CMUP 184.1 2.05 -0.47 1 


2.18 201.6 


35.74 2.54 


164386/2366320 


25.30 0.84 0.45 


CMUP 35.43 0.34 0.39 1 


0.49 14.67 


39.06 -2.23 


186145/10202075 


+27.36 -2.31 -0.33 


CMUP 26.58 2.81 0.29 1 


1.16 32.93 


47.40 0.25 


185332/13796632 


+28.88 2.72 0.44 


CMX 38.52 -2.61 0.92 1 


0.96 33.46 


57.47 0.14 



Table A.l: List of muon + photon + events. Et, Ptj ^"^^ ^'^^ i'^ GeV . Column "T" 
shows if an event has been triggered by the high-Ex niuon trigger. "+" in front of Ei^ value means 
that an event has been triggered by the high-Ex photon trigger. 
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/: 




)5-r 


run/event 




type PT (f) ?7 T 


A(f) m£j 




177371/273794 


25.24 -1.91 -0.39 


CMP 40.37 2.10 -0.20 - 
CMUP 34.51 -0.49 -0.39 1 


2.27 58.09 
1.42 38.43 


7.84 1.33 


179043/9452199 


-F28.91 -0.03 0.29 


CMX 29.72 -2.47 0.79 1 
CMX 20.95 1.30 -0.64 1 


2.44 57.61 
1.34 37.71 


5.55 2.79 


152507/1700045 


34.62 0.99 -0.96 


CMX 65.31 -3.00 -0.73 1 
CMUP 27.43 0.06 -0.47 1 


2.30 87.55 
0.93 31.31 


7.18 -0.82 


153345/1348276 


+35.05 -0.13 -0.65 


CMX 25.75 -3.13 0.78 
CMUP 21.63 0.27 -0.37 1 


3.00 75.63 
0.40 12.46 


1.85 -2.04 


154175/360663 


-K32.17 -1.23 0.12 


CMP 29.45 0.76 0.43 - 
CMX 28.42 2.99 0.75 1 


1.99 52.65 
2.06 55.60 


3.67 -2.65 


155895/2377214 


+26.58 -0.21 0.70 


CMUP 29.76 -1.33 0.33 1 
CMIO 23.10 0.00 -0.54 - 


1.12 31.49 
0.22 39.49 


16.3 -3.08 


156089/1783191 


-1-27.40 -2.75 -0.84 


CMX 39.87 0.50 -0.73 1 

CMUP 25.97 -2.71 -0.38 1 


3.03 66.16 

0.04 13.48 


4.24 -0.63 


161330/3293805 


-F27.42 -0.22 0.32 


CMUP 42.07 2.56 0.11 1 
CMIO 21.70 -0.70 0.63 - 


2.77 66.85 
0.48 13.49 


8.00 2.89 


162479/90695 


-1-26.73 -0.29 -U.DO 


CMX 46.39 1.95 -U.91 1 
CMIO 41.69 -1.37 -1.00 - 


2.24 64.14 
1.08 36.36 


3.76 2.27 


162238/185410 


+28.27 -2.27 -0.24 


CMUP 59.06 -0.12 -0.26 
CMP 40.10 1.83 0.47 - 


2.14 71.82 
2.19 66.24 


0.49 0.15 


164261/192043 


+25.70 -0.68 -0.20 


CMUP 37.11 2.82 0.47 1 
CMUP 20.54 -0.13 -0.25 1 


2.79 65.35 
0.55 12.42 


2.71 2.21 


167955/33039 


+26.68 -2.66 0.58 


CMUP 55.97 -0.04 0.31 1 

CMIO 31.75 2.61 0.60 - 


2.63 75.52 

1.01 29.11 


4.43 0.32 



Table A. 2: List of multi-muon + photon events. Et, Pt, Pt "^^7 GeV. Column "T" 

shows if an event has been triggered by the high-ET muon trigger. "+" in front of Ei^ value means 
that an event has been triggered by the high-ET photon trigger. 
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7 


e 


£7 




run/event 


Ex (t> f] 


type pT 4> 7/ T 


A(p nii^ 




178537/410195 


+54.36 -3.09 0.97 


TCEM 31.41 -1.11 0.53 1 


1.98 71.53 


67.39 0.60 


178677/5534996 


+38.69 -0.13 0.39 


TCEM 37.70 1.27 0.07 


1.40 50.56 


49.60 -2.45 


178758/3669012 


+53.70 -2.28 -0.76 


TCEM 27.63 -0.46 -0.86 1 


1.82 60.93 


40.55 1.38 


178785/12963076 


+58.01 -3.09 -0.62 


TCEM 28.12 -2.63 0.77 1 


0.46 63.80 


56.56 0.30 


151515/2273934 


+34.33 -3.12 -0.18 


TCEM 51.43 0.86 0.67 1 


2.30 85.51 


37.75 -1.54 


152602/728988 


+27.60 3.00 -0.19 


TCEM 57.76 -0.67 -0.80 1 


2.61 81.15 


34.44 1.37 


156083/1091219 


+35.04 2.13 -0.19 


TCEM 33.24 1.07 0.37 1 


1.06 39.69 


65.11 -1.57 


161170/427847 


+26.24 -2.29 -0.27 


TCEM 25.04 0.44 -0.15 1 


2.73 50.25 


34.92 -2.44 


163064/10108920 


+81.04 -1.15 -0.12 


TCEM 70.13 1.64 1.07 1 


2.80 178.7 


29.38 1.58 


163431/1462399 


+31.00 -0.54 -0.40 


TCEM 27.76 0.03 -0.32 1 


0.57 16.82 


57.92 2.90 


163526/14750 


+25.41 -2.06 -0.46 


TCEM 29.98 1.34 -0.54 


2.88 54.76 


29.75 -1.00 


167299/670904 


+27.78 2.29 -0.77 


TCEM 26.16 1.92 -0.43 1 


0.37 13.47 


43.32 -1.04 


167849/2063706 


+48.13 3.12 0.89 


TCEM 41.48 1.28 -0.14 1 


1.84 84.77 


42.78 -0.70 


168599/3868597 


+32.85 0.89 0.72 


TCEM 53.58 -1.71 -0.66 1 


2.60 101.9 


28.62 2.19 


183752/4059116 


+46.67 -0.50 -0.28 


TCEM 45.66 2.00 0.76 1 


2.51 101.0 


27.34 -2.60 


184762/1221041 


+30.45 2.64 0.47 


TCEM 59.61 2.95 0.75 1 


0.31 18.43 


101.3 -0.32 


183965/5394458 


+36.76 0.45 0.68 


TCEM 33.65 -1.72 0.80 1 


2.17 62.32 


28.56 2.48 


184519/1108274 


+33.70 -2.64 -0.50 


TCEM 41.63 -0.84 0.18 1 


1.81 64.29 


35.83 2.10 


184453/1736470 


+28.81 1.89 0.35 


TCEM 53.71 2.93 0.49 1 


1.04 39.35 


36.29 -1.22 


184778/6449604 


+28.14 -0.85 -0.71 


HCEM 176.9 -0.09 -0.59 1 


0.77 53.58 


26.61 -3.07 


184067/336957 


+142.1 2.20 0.38 


TCEM 68.80 -3.11 -0.45 1 


0.97 125.0 


32.88 -0.44 


1848(38/4710858 


+37.C3 -0.32 -0.55 


TCEM 47.4(3 -3.05 0.31 1 


2.73 90.83 


33.97 1.18 


185176/11940 


+25.94 1.24 -0.91 


TCEM 44.78 2.44 0.44 1 


1.19 61.29 


41.11 -1.32 


184778/345250 


+28.03 2.82 0.13 


TCEM 33.77 1.54 -0.57 1 


1.28 42.98 


26.55 -1.21 


185848/2995941 


+43.42 -2.57 0.50 


TCEM 43.14 0.65 1.07 1 


3.06 90.03 


28.09 0.61 



Table A. 3: List of electron + photon + events. Et, Pt- f^T GeV. TCEM 

stands for Tight CEM electron. Column "T" shows if an event has been triggered by the high-Ex 
electron trigger. "+" in front of Ei^ value means that an event has been triggered by the high-Ex 
photon trigger. 
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7 


£ 






run /event 


Et ([> 1] 


type pt (t> rj T 


A(p m^y 


</> 


178602/3018939 


+25.72 2.74 -0.89 


TCEM 60.24 -0.68 -0.74 1 
TCEM 32.77 2.23 0.65 1 


2.85 78.14 
0.52 51.37 


5.95 -1.56 


178852/3194148 


+94.25 -2.93 0.64 


TCEM 101.7 0.58 0.65 1 
LCEM 31.52 -1.14 0.65 1 


2.77 192.5 
1.79 84.90 


6.18 -1.46 


153074/1339595 


+25.54 -0.33 -0.18 


TCEM 76.89 2.87 -1.02 1 
TCEM 38.95 -1.34 -0.99 1 


3.08 95.67 
1.01 39.36 


16.3 1.31 


160346/1528176 


+45.61 -0.82 -0.87 


TCEM 50.34 -3.12 -1.00 1 

LCEM 24.53 0.81 -0.94 1 


2.30 87.71 
1.63 48.68 


19.7 1.84 


155394/2469758 


+27.80 -1.77 -0.28 


WEST 52.56 1.12 1.82 1 
TCEM 29.65 -2.30 0.80 1 


2.89 123.4 
0.53 35.45 


1.62 1.45 


155996/1191192 


+37.75 1.21 -0.64 


TCEM 66.64 2.81 -1.06 
EAST 19.33 -0.27 -1.37 


1.60 75.27 
1.49 41.49 


7.42 -1.72 


161830/69435 


+135.0 -0.19 -0.57 


HCEM 138.5 3.07 0.82 1 
LCEM 26.95 1.91 -0.16 1 


3.02 340.2 
2.11 107.5 


12. 2.50 


162396/1323030 


+27.61 0.55 -0.99 


TCEM 39.96 -2.81 -0.95 1 

LCEM 20.18 0.15 -0.39 1 


2.92 66.01 

0.40 15.82 


2.94 1.13 


164844/6642760 


+26.85 -0.66 -0.07 


TCEM 57.45 -2.92 0.88 1 
TCEM 41.01 0.54 0.98 1 


2.26 80.95 
1.20 52.39 


9.50 1.53 


166038/6509526 


+31.01 2.27 0.36 


TCEM 27.90 -1.89 0.96 1 
LCEM 22.41 -0.13 -0.15 1 


2.13 55.19 
2.40 50.98 


4.32 1.68 


167623/4691216 


+27.49 1.92 0.70 


TCEM 28.78 1.70 0.21 1 
LCEM 23.13 -1.49 -0.72 1 


0.21 13.31 
2.88 62.58 


10.4 -1.32 


167866/443088 


+67.15 -1.26 -0.29 


TCEM 148.1 1.73 0.05 1 
TCEM 99.93 -1.29 0.54 1 


2.99 201.8 
0.03 70.28 


17.1 1.47 


183913/878106 


+30.10 -0.75 -0.92 


TCEM 42.48 2.21 0.61 1 
EAST 17.65 -1.26 -1.57 1 


2.96 91.99 
0.50 20.09 


4.56 2.77 


184519/3570367 


+27.58 2.43 0.44 


TCEM 35.08 -0.80 -0.18 1 

WEST 1C.02 0.97 1.28 1 


3.05 65.14 

1.4C 33.13 


9.03 -2.57 


185037/1000257 


+33.05 -1.15 -0.22 


EAST 51.01 1.45 -1.35 1 
TCEM 34.93 -2.46 0.93 1 


2.60 92.97 
1.30 59.26 


2.53 -1.22 


185075/1540099 


+64.57 2.85 0.17 


TCEM 59.86 -0.71 0.48 1 
WEST 18.92 1.35 1.78 1 


2.71 123.0 
1.51 78.68 


12.2 0.48 


185281/13145621 


+32.43 -1.69 -0.95 


TCEM 44.60 0.95 -0.40 1 
TCEM 33.45 2.93 -0.52 1 


2.64 77.33 
1.66 51.20 


5.00 -0.63 


185634/4852157 


+37.10 1.98 0.83 


TCEM 28.27 -1.85 0.56 1 


2.45 61.58 


1.85 -0.98 


Continued on the next page 
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Table A. 4 — continued 







TCEM 25.95 -0.12 0.36 1 


2.10 55.73 




186145/10988173 


+31.88 1.18 0.14 


TCEM 27.08 -2.73 0.58 1 
LCEM 22.91 -1.20 -0.68 1 


2.37 56.00 

2.38 55.12 


1.76 -0.85 



Tabic A. 4: List of multi-electron + photon events. Ex, Pt; 
;Et and niiy are in GeV. TCEM stands for Tight CEM Elec- 
tron. LCEM stands for Loose CEM Electron. HCEM stands 
for TightlOO CEM Electron. EAST stands for Phoenix East 
Plug Electron, WEST stands for Phoenix West Plug Elec- 
tron. Column "T" shows if an event has been triggered by 
the high-Ex electron trigger. in front of E^p value means 
that an event has been triggered by the high-Ex photon trig- 
ger. 
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A.2 Additional ^7^^ Plots 



In this section we present additional plots of the identification variables for /i7^T ^I'^t 
signatures. 



A. 2.1 Additional ^^7^^ Plots 

CDF Run II Preliminary Jl dt = 305 pb"^ 
8F 



Entries 


17 


Overflow 







0.5 1 1.5 2 2.5 3 3.5 4 4.5 

li^jj AR(nY) 



CDF Run II Preliminary | L dt = 305 pb~ 

14 

12: 
10- 

U) 

> 

LU 5 



Entries 


17 


Overflow 






. Data 
□Wy 

■ Wjet 

"Zyy, Wyy 

■ Non-WZ 



1 2 3 4 5 6 7 
Muon Type (3:CMUP, 4:CMX) 

(c) Muon Detector Type 



CDF Run II Preliminary Jl dt = 305 pb"^ 



Entries 


17 


Overflow 







100 200 300 400 

H^tY HT(GeV) 



500 



CDF Run II Preliminary Jl dt = 305 pb"^ 




Entries 


17 


Overflow 






. Data 
□Wy 
□Zy 

■ Wjet 
• Zyy, Wyy 
|Ty 

■ Non-WZ 



-0.5 0.5 

Muonri 

(d) Muon 77 



Figure A.l: Muon + photon + distributions: AR{fi'j), Ht, Detector Type(/i), ri{fj,). 
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Figure A. 2: Muon + photon + distributions: do{^). 
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A.2.2 Additional e7^T Plots 
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Figure A. 3: Electron + photon + distributions: Ai?(e7), Ht, fiduciality (see Chapter [6]), 
electron t]. There are 3 electron candidates non-fiducial in central or plug, which is in the agreement 
with the expectation. 
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Figure A. 4: Electron + photon + distributions: CES X, CES Z, Acot{9) and Axy (see 
Section rA.6p . There is only one trident in the 67 sample. 
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A.3 Additional U-f Plots 



In this section we present additional plots of the identification variables for /i/x7 and ee7 signatures. 



A.3.1 Additional ^i^i^y Plots 
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Figure A. 5: Multi-muon + photon distributions: AR^fi^y), Ht, vif^)^ do{fi) 



139 



CDF Run II Preliminary | L dt = 305 pb 
9^ 



'I' 



Entries 


24 


Overflow 







2 3 4 
Muon COT Chi2 



CDF Run II Preliminary J L dt = 305 pb"' 



Entries 


24 


Overflow 






20 



-15 
■^10 



• Data 
□zy 

■ Zjet 
Zyy 

■ Non-WZ 



12 3 4 
Muon Axial COT Layers 

(c) Muon N Axial Segments 



CDF Run II Preliminary | L dt = 305 pb 

9^ 



Entries 


24 


Overflow 







4 6 8 10 
Muon Si IHits 



CDF Run II Preliminary J L dt = 305 pb"' 

25|_i i i 1 1 i I i i I i i i i 1 1 i i 1 1 i 1 1 i I i i I i 1 1 i 1 1 1 1 1 i 1 1 i 



14 










03/25/06 


Entries 


24 


Overflow 






20 



i2l5 
c 

V 

■^lO 



• Data 
□Zy 

■ Zjet 
»Zyy 

■ Non-WZ 



1 2 3 4! 
Muon Stereo COT Layers 

(d) Muon N Stereo Segments 



Figure A. 6: Multi-muon + photon distributions: Number of si hits, Numbers of Axial and 

Stereo Segments in COT 
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Figure A. 7: Multi-muon + photon distributions: detector type (/i), AX (CMU), AX (CMP), 
AX (CMX) 



141 



'I' 



«10 
I 8 

LU 

6 
4 
2 
Q 



CDF Run II Preliminary | L dt = 305 pb 

16- 
14 
12 



Entries 


24 


Overflow 






0.02 



0.04 0.06 
Muon Iso Et 



CDF Run II Preliminary L dt = 305 pb' 



Entries 


24 


Overflow 






• Data 

□zy 

■ Zjet 
Zyy 

■ Non-WZ 



0.08 0.1 




0.2 0.3 
Muon Iso Pt 



(a) Muon 



CDF Run II Preliminary L dt = 305 pb 



Entries 


12 


Overflow 







12l Muon Iso Pt 



0.5 1 1.5 2 2.5 
3D Angle(Acollinearity) 



(c) Muon AcoUinearity (3D Angle(/i/i)) 



Figure A. 8: Multi-muon + photon distributions: muon relative calorimeter isolation (Iso Et), 
relative track isolation (Iso px), 3D angle(/x/x) 
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A. 3. 2 Additional 667 Plots 
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Figure A. 9: Multi-electron + photon distributions: AR{e'j), Ht, electron fiduciality (see Chap- 
ter [6]), 7]{e). 
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Figure A. 10: Multi-electron + photon distributions: CES X, CES Z, Acot{6) and Axy (see 
Section [A. 61) . There is one conversion electron and two tridents in the ee'j sample. 
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A.4 Stability Plots for and Wj 



To check electron and muon identification in events with an extra object (such as a photon in 
the signal channel) we plot the rate for Zj and W j in the 8 bins of luminosity (Table I4.ip in 
Figure lA.llI 
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Figure A.ll: Stability plots for Zj (a), Wj (b) for muons; Zj (c), Wj (d) for electrons. The bins 
are those of Table 14. 1[ The uncertainties are statistical only. The luminosity systematic error of 
6% (Chapter [TTjl is not included 



145 



A. 5 CMX vs CMUP muons: Comparison of Isolation Vari- 
ables 



We have checked that muons that go into the CMX system have the same isolation properties as 
CMUP muons. Figure lA. 121 shows the distributions in calorimeter isolation, relative track isolation 
(total pt of tracks in a cone in rj — (f space of radius R = 0.4 around the muon track divided by 
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Figure A. 12: The distributions for CMX (black histogram) and CMUP(red points) muons of 
calorimeter isolation, relative track isolation (see text), absolute track isolation, for Z° P-^fJ'^ 
events (top row) and — > /i^i^ events (bottom row). For the Z° f^^fJ'~^ plots both muons are 
required to be either CMUP or CMX. 
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A. 6 Fake Electrons from Photon Conversions 



There are three dominant sources of fake electrons: a) photons from 7r°, t] and other mesons, that 
convert into asymmetric e^e~ pairs in the material before the COT volume, b) charged hadrons 
in jets that either interact in the electromagnetic volume of the calorimeter or overlap with a 7r° 
or secondary photon in the jet, and c) electrons from the decay of heavy flavor (6, c, and maybe 
even s) . We estimated these backgrounds in Section 110.21 by studying the total pt of tracks in a 
cone in 1] — if space of radius R = OA around the lepton track. 

We consider fake electrons from photon conversions below. Electrons coming from photon 
conversion are identified by conversion algorithm, which looks for couple of opposite sign tracks 
with \Axy\ <0.2 cm and \ Acot{e) \ <0.04. 

For each electron a conversion flag is tested. We define if the electron is flagged as coming from 
a conversion (7 e^e~) or from trident events where a conversion is caused by a bremsstrahlung 
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Figure A. 13: The distributions for same-sign (red points) and opposite-sign (black histogram) 
e^e~ pairs in invariant mass, Pt, and A0; each distribution is shown twice, in linear plots(a, b, 
c), and in log plots(d, e, f). 



Figure IA.13I shows the distributions for same-sign and opposite-sign e^e pairs. The invariant 
mass of the same-sign electrons is shifted with respect to the invariant mass of Z° — >■ e^e~ (Fig- 
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Figure A. 14: The distributions for same-sign (red) and opposite-sign (black) e~^e~ pairs in in- 
variant mass for Z" — e~^e~ MC sample. The observed behavior is similar to what we have in 
data. 

ure IA.131 a). We observe similar behavior (shift in invariant mass distribution) in the Z° e^e~ 
MC as in data (Figure lATH]). 

To develop an understanding of the conversion cut for the candidate events, we summarize the 
same-sign events in Z° e^e~ sample. Out of 199 events in this sample only 21 are not tagged 
as either a conversion electron or a trident. We find that we have 5 same sign muon events, which 
is comparable to 21 non-conversion/trident same-sign electron events. 
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Table A. 5: A breakdown of the source of same-sign electrons in the Z° — > e^e sample. Most 
same-sign events are tagged by the conversion filter as conversions or tridents. 
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A. 7 Checking the /i/i7 and /i7pT Additional Backgrounds 



We have used a number of techniques, described below, to look for additional backgrounds in the 
/i/i7 and /i7^T samples. 

Same-Sign Leptons to Estimate Jets Faking one or More Leptons 

We used results from Ref. [HI] on the numbers of same-sign (SS) and other-sign (OS) muon pairs 
in the dimuon sample. The expected background from SS muon pairs to fifi'-f is calculated as 
follows. 

First we obtain the ratio of W+1 jet events to OS events, x, expected in the dimuon sample: 
X = SS/OS X 1.51±0.05. For the ^fi sample ratio of SS/OS is of order of 0.05% and therefore 
x=0.1%. In the fifi'-f sample we have 12 OS events, and therefore the expected background from 
W+1 jet is negligible, 0.1% x 12 = 0.012±0.001 events. 

Decays in flight of Low Momentum K"^ Faking a High-Momentum Muon 

A low-momentum hadron, not in an energetic jet, can decay to a muon forming a "kink" between 
the hadron and muon trajectories (Figure [A.lSp . In this case a high- momentum track may be 
reconstructed from the initial track segment due to the hadron and the secondary track segment 
from the muon. A kaon that decays before the COT volume results in a muon whose momentum 
is correctly measured; a kaon that decays after the COT is itself correctly measured. These 
contributions are included in the total background estimate (see Section [10.2p . 

The contribution from this background is estimated by identifying tracks consistent with a 
"kink" in the COT. We count the number of times that, proceeding radially along a COT track, 
a "hit" in the n+1 layer of sense- wires is on the other side of the fitted track from the hit in the 
nth layer. Real tracks will have hits distributed on both sides of the fit, and will therefore have 
many "transitions". A mis- measured track from a 5-GeV (for example), on the other hand, 
will consist of two intersecting low-momentum arcs fit by a high momentum track, and will have 
a small number of transitions [85j . 

Figure IA.16I shows the number of transitions in muons in the Z° f^^fJ'^ control sample, and 
in a sample enriched in hadron decays by selecting events with a large > 25 GeV, at least 
one jet and muon that have large impact parameter do > 0.2 cm. Figure IA.17I shows the number 
of transitions for muon tracks with and without silicon hits. The red curve is the distribution for 
the muons from sample. Decays-in-flight have a distribution that peaks much lower, 

with few events above 30 transitions. We see no evidence that any of these tracks are DIF muons. 
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Reconstructed "track" 



Figure A. 15: Decays-In-Flight: 
schematic figure, K fiu. Two 
track segments from K and /i misre- 
constructed as one track. 
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Figure A. 17: Number of transitions for fifi'j (a) fi tracks with SVX hits, (b)/i tracks without SVX 
hits; /i7pT (c) A* tracks with SVX hits, (d) fi tracks without SVX hits. Muons from Z° fJ'^fJ'~ 
are shown as red histogram. 



/i/i7, fJ^'y^T Cosmics Background 

For this we invert cosmic cut, and require the event to be tagged as cosmic [59]. We processed 
the unstripped high-px muon sample with this inverted requirement and found no /i/i7 or /i7^x 
candidate events. 

In addition we scanned our /i7^x ^^"^ f^f^l candidate events with CDF Run II Event Display 
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(Section [3]) and made sure that none of them look hke beam halo or cosmics events. 
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